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MEASUREMENTS IN LIQUID FUEL SPRAYS 
Norman Chigier 

Department of Mechanical Engineering 
Carnegie-Mel Ion University 


The goal of our research on the measurement of liquid fuel sprays is to study 
the effects of atomizer design and to characterize the atomization and vaporization 
processes in sprays under well controlled conditions. The scope of this study is 
limited to the events directly preceding combustion in the liquid fuel sprays. 
Measurement techniques are being used to provide information as a function of space 
and time on droplet size, shape, number density, position, angle of flight and 
velocity. Three spray chambers have been designed and constructed for: (i) air- 
assist liquid fuel "research" sprays; (ii) high pressure and temperature chamber for 
pulsed diesel fuel sprays; and (iii) coal-water slurry sprays. Recent results utiliz- 
ing photography, cinematography, and calibration of the Malvern particle sizer are 
reported. Systems for simultaneous measurement of velocity and particle size dis- 
tributions using Laser Doppler Anemometry/ Interferometry and the application of holo- 
graphy in liquid fuel sprays are being calibrated. 

Parker Hannifin has designed an air assist swirl atomizer for use in basic re- 
search studies at several universities (Carnegie-Mel Ion, Purdue, UC/lrvine). This 
atomizer has been fitted into the basic research spray chamber in which uniform low 
turbulence intensity air flows vertically downward. The liquid fuel spray is in- 
jected downward along with the uniform air flow. Detailed measurements of the dis- 
tributions of particle size, particle velocity, and gas velocity will be made through- 
out the sprays under conditions of well controlled initial and boundary conditions. 

Design and construction of the high pressure and high temperature swirl chamber 
has been completed to withstand maximum pressures of 60 atm. and maximum temperatures 
of 370°C. Two quartz windows, 89mm in diameter, are fitted to the chamber for opti- 
cal access. An electromagneti cal ly controlled pulsed diesel injector has been 
operated with the chamber and preliminary photographic studies have been made. The 
air-assist swirl atomizer will, at a later stage, be installed in this high pressure 
chamber . 

The coal-water slurry spray test facility was operated with a Parker Hannifin 
counter swirl atomizer to investigate fundamental processes leading to combustion 
of coal-water slurries. This research is focuseJ on two areas: the vaporization of 
water and the devolatilization of the coal particles in the spray. Preliminary 
studies on drop size, shape, and distribution were made using photography, cine- 
matography, and the Malvern particle sizer. High quality images of commercial coal- 
water (70* to 30* concentration) slurry spray at various atomizing air pressure and 
slurry feed rate have been made. The magnification of these pictures is 8.4 times 
lifesize and is about l inch distance from the nozzle tip. The spray appeared very 
dense in the photographs. The smallest particle/droplet size that can be distin- 
guished on the edges of the spray is in the order of 40-50 microns. There is evi- 
dence of low-frequency pulsations in the spray in both photographic and high speed 
cinematographic results. 

Samples of diluted coal-water slurry were measured by the Malvern particle 
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sizer for size distribution. The results showed that the coal particle SMD fell 
within the size range of 2 microns to 1.5 microns for an obscuration ranging from 
0 to 90 percent. 


The behavior ot the Malvern 22l\l particle sizer was investigated by analyzing 
a non-varying calibration reticle produced by Laser-K 1 ect ro optics. Some of the 
important parameters varied during the experiment were: photomask tilt angle, 
distance between photomask and receiver lens, and receiver lens focal length. All 
of these variables proved to be very crucial in the determination of size distribu- 
tion by the Malvern. 
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Thre*-Olmens»onal View of the Experimental System 
to Study Coal Water Slurry Sprays 




Fuel: 50% Ethylene Glycol, 50% Water by volume. 
Fuel pressure: 21 psi, Nozzle: 60 degree solid cone. 


C" 

GF I 


5 



Fuel: 50% Ethylene Glycol, 50% Water by volume. 
Fuel pressure: 40 psi, Nozzle: 45 degree solid cone. 
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Fuel: 50% Ethylene (Jlycol, 50% Water by volume. 
Fuel pressure: SO psi. Nozzle: 45 decree solid cone. 
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Slurry: 210 lb/hr, 22 psi; Ain 44 psi 


9 





ORIGINAL PAGE 59 
OF POOR QUALITY 



Slurry: 360 lb /hr, 42 psi; Ain 73 psi 
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Slurry: 300 lb /hr, 50 psi; 


Ain 132 psi 
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FULL SPRAY DIAGNOSI1CS 

Marla A. Bosque 
NASA Lewis Research Center 


lo understand the combustion process of liquid reactants. It Is necessary 
to obtain Information about the distribution and atomization of the Injected 
liquid phase. Turbulence and relative velocities between the fuel's gas phase 
and liquid droplets affect the evaporation, burning rate and pollutant forma 
tlon. This will affect other combustor performance parameters such as combus 
tlon efficiency, pollutant emissions and combustion stability. The measurement 
of size and velocity distribution of the liquid Injected Into the combustor Is 
of vital necessity. The experimental data will be used by numerical modelers 
In deriving mathematical relationships that will describe the combustion 
process . 

Several laser measurement methods are being studied to provide the capa- 
bility to make droplet size and velocity measurements under a variety of spray 
conditions. The Droplet Sizing Interferometer (DSI) has promised to be a suc- 
cessful technique because of Its capability for rapid data acquisition, com- 
pilation and analysis. Its main advantage Is the ability to obtain size and 
velocity measurements In air fuel mixing studies and hot flows. 

Ihe existing DSI at NASA Lewis Is a two-color, two component system. Two 
Independent orthogonal measurements of size and velocity components can be made 
simultaneously. It also uses an off axis large angle light scatter detection. 
The fundamental features of the system are optics, signal processing and data 
management system. The major component Include a transmitter unit, two 
receiver units, two signal processors, two data management systems, two Bragg 
cell systems, two printer/plotters, a laser, power supply and color monitor. 

Ihe laser diagnostic Investigation will consist of various experiment 
configurations for fuel nozzle character 1 zat Ion tests. Water will be used as 
a fuel substitute. General probe volume positions will be obtained for mapping 
the fuel nozzle flow distribution and variations with test conditions. A fuel 
nozzle/swl rler combination will simulate fuel/alr mixing characteristics; small 
particles will be used for seeding the flow field. 

lwo contracts have been awarded for further development and Improvements 
to the system. These contract efforts promise to extend the size range capa- 
bility, reduce beam alignment difficulties and reduce the system sensitivity 
to laser beam quality and differences In the relative Intensity of the beams. 
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fuel spray diagnostics program 

PURPOSE •CONDUCT FUNDAMENTAL COMBUSTION PROCESS RESEARCH 

USING DIAGNOSTIC MEASUREMENTS FROM A UNIQUE LASER 
INTERFEROMETER APPARATUS. 

METHOD • N0N-INTRUS1VE MEASUREMENT PROBE FORMED BY 

COHERENT LASER BEAMS THAT PRODUCE INTERFERENCE 
FRINGE RANDS. 

•DROPLETS/PARTICLES GENERATE LIGHT SCATTER PATTERNS 
WHICH ARE DETECTED WITH PMT . 

YIELD • FAST REAL TIME DATA ANALYSIS OF FLOW CHARACTERISTICS. 
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DROPLET SIZING INTERFEROMETER SYSTEM LINE DRAWING (ONE-COMPONENT) 
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NASA LeRC FUEL SPRAY DIAGNOSTICS PROGRAM 

EXPERIMENT DESCRIPTION 


• SPRAY CHARACTERIZATION - HjO 

-VARIOUS NOZZLE TYPES 

- PRESSURE DIFFERENTIAL EFFECTS SIZE AND VELOCITY DISTRIBUTION 

- SEVERAL SAMPLING PLANES 

• TURBULENCE STUDIES - SEEDED 

-VARIOUS AIR SWIRLER CONFIGURATIONS 

- VANE ANGLE EFFECTS 

- AIRFLOW SPLITS 

• FUEL-AIR REACTION STUDIES 

- VARIOUS EQUIVALENCE RATIOS 

- FUEL PROPERTY EFFECTS 
-MEASURE TURBULENCE INTENSITY 
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RESOURCES 

• FUEL SPRAY DIAGNOSTICS FACILITY 
(SW-5, ERB , NASA LEWIS ) 

• EEVQOftENT AND IMP! Eli NT AT I ON OF AIM ANTED DIAGNOSTIC TECHNIQUES 

CSPECTRON DEVELOPt-ENT LABORATORIES, INC., CALIFORNIA) 

• mASE DETECTION TECHNIQUE STUDY 
(AERUtTRlCSJNC. , CALIFORNIA) 

• FEASUREFENT OF SPRAY CC1BUSTI0N PROCESSES 

(UTS! GRANT, TENNESSEE) 

• CfimriCJN CHARACIERIST1CS IN THE TRANSITION REGION OF LIQUID FLUID SPRAYS 

(DREXEL UNIVERSITY GRANT, PHILADELPHIA) 
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DEVELOPMENT AND IMPLEMENTATION OF ADVANCED 
DIAGNOSTIC TECHNIQUES* 

Cecil F. Hess 

Spectron Development Laboratories, Inc. 


The objectives of this work are to establish the performance and 
limitations of two optical techniques for spray characterization, and to design 
and build an instrument based on these techniques. 

Two techniques have been identified which offer great potential in the 
measurement of sprays. The first is referred to as "1MAX", and it consists of 
a nonintrusive pulse height analyzer. The second is referred to as "Visibil- 
ity/Intensity (V/I)” and it performs a size measurement by examining the visi- 
bility and the pedestal intensity of a Doppler burst. The research conducted 
over this past year indicates that the IMAX technique provides a larger dynamic 
range and higher accuracy than V/I. It also shows that the two-color IMAX 
concept provides a higher S/N primarily because of the high efficiency in 
spectrally separating the two signals. 

The two-color IMAX concept is described in Figure la. Two small beams of a 
given wavelength (4880A) are crossed in the middle of a larger beam of differ- 
ent wavelength (5145A) thus identifying a region of almost uniform intensity 
within the large beam. The two small beams will interfere where they cross and 
a fringe pattern will be formed in the middle of the large beam. Droplets 
crossing the fringes also cross the middle of the large beam. Since the peak 
intensity of the large Gaussian beam is known, a unique relationship between 
droplet size and scattered light is established. The velocity of the droplets 
is also measured using the classical Doppler approach. Figure lb shows a 
schematic representation of the breadboard system used to acquire the reported 
data. 


The size distribution of two kinds of sprays are reported here. The first 
is produced by a Berglur.d-Liu droplet generator with dispersion air. Monodis- 
perse, bimodal, and trimodal sprays with an angle of about 10° were thus pro- 
duced, and the results are shown on Figure 2a. The theoretically predicted 
sizes are 49 pm, 62 pm, and 70 pm, respectively. Good accuracy and resolution 
can, therefore, be observed. The second spray was produced by a pressure 
nozzle. The results are shown on Figure 2b. In order to test the resolution of 
the system, data were obtained using three different size ranges: 5 to 50 pm; 

10 to 100 pm; and 20 to 200 pm. This is one of the most difficult self- 
consistency tests imposed on any technique, and most available techniques will 
show a shift in the predicted data. IMAX shows excellent matching of the data. 

The Visibility/Intensity Technique is described in Figure 3. This techni- 
que makes use not only of the visibility of a Doppler signal but also the peak 
intensity of the pedestal. Both parameters are available in the signal and 
their cross-correlation can be used to eliminate faulty signals produced in 
many practical environments. This technique will especially prevent small 


*Contract No. NAS3-23538 
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particles (high visibility) from appearing as large. Figure 3 shows a Doppler 
trace where both the visibility and the peak intensity of the pedestal are 
indicated. There is a relationship between the size of the droplet and the 
amount of scattered light given by Mle theory. This relationship can be used 
to eliminate signals with an apparently different size by using the following 
logic: Droplets that produce a certain visibility are associated with a given 

size; hence, they should scatter light with a given intensity (characterized by 
1 ). Two exceptions are contemplated: First, droplet with the correct visi- 

bility will scatter different amounts of light due to the Gaussian nature of 
the probe volume's intensity; second, droplets with an erroneous visibility 
will not scatter light with an intensity corresponding to their apparent size. 

Results for both Visibility/Intensity and visibility only are shown in 
Figure 4. Figure 4a shows the data of the spray formed of primary droplets and 
doublets. Figure 4b shows results similar to Figure 4a but obtained with 
visibility only. Notice how much broader the later distributions are. Figure 
4c shows a spray containing primary droplets, doublets and triplets. Figure 4d 
shows results similar to Figure 4c but obtained with visibility only. Notice 
that the later distribution is broader and the resolution is not as fine. 

Both the accuracy and resolution of these measurements are very good. The 
theoretically predicted sizes are S3 pm, 6b pm, and 76 pm. The corresponding 
measured diameters are (32 to 54), (63 to 65), and 73 pm. 

Larger errors (about 20%) can be expected when measuring the size at a 
higher visibility (80%). 

Based on the above concepts, an Advanced Droplet Sizing System (ADSS) was 
developed. Photographs of the optics and electronics are shown in Figure 5. 
Very exhaustive tests are presently being conducted to establish and ensure the 
performance of the ADSS. 
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Figure 5. Electronic Components and Optical Heads 
of the Advanced Droplet Sizing System 
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ANALYSIS AND TESTING OF A NEW METHOD FOR DROP SIZE AND VELOCITY 

MEASUREMENTS USING LASER LIGHT SCATTER INTERFEROMETRY 

W. D. Bachalo 
Aerometrics, Inc. 


Research was conducted on a laser light scatter detection method for measuring 
the size and velocity distributions of spray drops. The work was motivated by the 
need for accurate fuel spray data in a variety of research and development activities 
related to combustion applications. Presently, there is considerable interest in 
obtaining measurements of the fuel drop and air motions within the spray combustor. 

It is known that the turbulent motions influence the reactions by increasing the oxy- 
gen supply to the burning fuel. Thus, the droplet size and velocity distributions 
are important data to be obtained when characterizing the combustion process. The 
need for simultaneous size and velocity measurements has been recognized for even 
general characterizations of sprays. This is especially true when the nozzles are 
operating in turbulent flow environments or when the nozzles generate drops that are 
moving at relatively high velocities and subsequently relax to the ambient flow 
speed in accord with their initial momentum. 

The requirement for simultaneous measurements of drop size and velocity has led 
to the combination of the laser Doppler velocimeter (LDV) with various particle siz- 
ing methods 1 »•’ » 3 . However, experience has shown that the available methods are 
adversely affected by the measurement environment. The conditions include high drop 
number densities, beam attenuation by the surrounding drops and optical access ports, 
and the presence of small combustion generated particulate. Many of these measure- 
ment limitations have been eliminated by deriving a method that performs the measure- 
ments independent of the beam intensity. 

The method consists of an optical system which is the same as a LDV except that 
three detectors are located at selected spacings behind the receiver aperture. Drops 
passing through the intersection of the two beams scatter light which produces an 
interference fringe pattern. The temporal frequency of the fringe pattern is the 
Doppler difference frequency which is linearly proportional to the drop velocity 
while the spatial frequency is linearly related to the drop diameter. The spacing of 
the fringes is also dependent upon the light wavelength, beam intersection angle, 
drop refractive index (unless reflected light is measured), and the location of the 
receiver. Measurement of the spacing of the fringe pattern produced by the scattered 
light may be achieved by placing pairs of detectors at selected spacings in the 
fringe pattern or its image. As the fringes move past the detector at the Doppler 
difference frequency, the detectors produce identical signals but with a phase shift 
proportional to the fringe spacing. The utilization of three detectors ensures that 
phase ambiguity does not occur, provides redundant measurements for signal valida- 
tion and allows an expanded operating range while maintaining good sensitivity. 

An experimental effort was conducted to verify the theoretical analysis, eval- 
uate signal phase processing methods, and investigate the effects of signal quality 
upon the measurements. The optical system was arranged to provide the flexibility 
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needed in evaluating the parameters that influence the measurements. Only two detec- 
tors were used in the present investigations. After investigating several signal 
processing concepts, a breadboard processor was developed and interfaced to an IBM 
microcomputer. Simultaneous phase and Doppler period measurements were made and 
transmitted to the computer to form size and velocity distributions. 

Several drop sources were used to establish the test conditions. A Berglund-Liu 
generator was used to produce monodisperse drop streams. The monodisperse stream 
could be directed through various parts of the measurement volume. Several spray 
nozzles were also used including pressure atomizers of moderate flow rates, high 
volume nozzles that produced large drops and spinning disk atomizers that produced 
small drops in a narrow distribution. 

Initially, tests were performed using the monodisperse drops to verify the theo- 
retical analysis. The results obtained were in complete agreement with the theory 
for all of the optical parameters tested. The parameters included the laser beam 
intersection angle, the detector spacing, and light scatter detection angle. Off- 
axis backscatter was also evaluated. These tests showed that the light scattered 
internally was measured when the drops were transparent. Dye was used to produce 
opaque drops which then resulted in the measurement of reflected light. Measure- 
ments based on reflected light can be made independent of the drop index of refrac- 
tion. 


Spray interference effects were evaluated by passing the laser beams through the 
spray and measuring the monodisperse drop stream. This had an insignificant effect 
on the measurements. The small differences observed were in most part a result of 
the spray interference with the monodisperse droplet stream. The effect of the 
attenuation of the individual laser beams was tested by attenuating one of the beams 
by 90"'.. Such attenuations can occur, for example, when large drops pass through the 
beams outside of the measurement volume. As expected, this had no effect on the 
measurements . 

Sprays were measured to assess the performance of the basic system in real 
environments. In the first case, the spinning disk atomizer was used to generate 
drops in the small size range. Measurements were obtained at three size range set- 
tings and the results compared to determine the consistency of the data. The agree- 
ment was good and the mean size agreed with the results obtained with other methods. 
Direct comparisons were also made of the measurements obtained in a spray generated 
by a pressure atomizer. The data obtained by Delavan Inc. and our method showed 
good agreement of both the mean size and the distribution. 

In summary, the recognized characteristics of the method are: 

• linear relationship between the measured phase angle and drop size 

• size range of approximately 100 at a single optical setting 

• simultaneous size and velocity measurement 

• relative insensitivity to beam or light scatter attenuation 

• high spatial resolution 

• operation is similar to an LDV 
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• adaptable to existing LDV systems 

• can perform measurements independent of refractive index 

An instrument is currently being developed for delivery to NASA Lewis. Develop- 
ment work is continuing on the refinement of the analysis, sample volume characteri- 
zation, and the mass flow measurement capability. Further evaluations of the method 
will be made in combustion environments. The development of the technology for per- 
forming two-component velocity measurements in two-phase flows is required. 
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DEVELOPMENT GOALS 


PRODUCE AN INSTRUMENT TO ACCURATELY MEASURE DROP SIZE AND 
VELOCITY DISTRIBUTIONS TO ADDRESS THE FOLLOWING: 

• SPRAY CHARACTERIZATIONS 

• SPRAY CHARACTERIZATIONS IN COMPLEX TURBULENT 
FLOWS 

• TWO-PHASE TURBULENT FLOW MEASUREMENTS 

• SPRAY DROP MEASUREMENTS IN COMBUSTION 


PHASE/DOPPLER SPRAY ANALYZER 

OBJECTIVE: INVESTIGATE THE CHARACTERISTICS OF A NEW METHOD 

FOR MEASUPING THE DROP SIZE AND VELOCITY 
DISTRIBUTIONS IN SPRAYS 

TASKS: 

. PRODUCE THEORETICAL DESCRIPTIONS OF THE LIGHT 
SCATTERED BY SPHERES AT THE INTERSECTION OF 
TWO LASER BEAMS 

• GENERATE THE THEORETICAL RELATIONSHIP BETWEEN 
THE MEASURED INTERFERENCE FRINGE PATTERN AND 
THE DROP SIZE 

. EXPERIMENTALLY VERIFY THE THEORETICAL ANALYSIS 
USING MONODISPERSE DROPS 

• EVALUATE SIGNAL PROCESSING METHODS 

t TEST THE SELECTED PROCESSING METHOD 
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MEASUREMENT 


DET 1 
DET2 
DET 3 

SCHEMATIC OF THE OPTICAL SYSTEM FOR AN LDV AND THE DROPLET SIZING SYSTEM. 

35 









ORIGINAL PAGE ;g 
OF POOR QUALITY 


/ \ \ 


Wmmxr 

ymwMS 

AM/WVW/Y 


rAlWfWAVAVAVfVfVAWH 


L'OMPUTIP INI I mi Rl Nil I R 1 Niil PATTI RN PRODUCT P RY A PROPUT 


HIGH PASS MITIREP POP PI I R RURST SIGNALS ILLUSTRATING THI PHASE SH1TT PETWLEN 

DETECTORS . 




ORIGINAL PAGE IS 
OF POOR QUALITY 



RtuATIONSHIP FOR THREE DETECTORS. 



COMPARISONS WITH EXPERIMENT. 

THEORETICAL PREDICTION SHOWING THE PHASE VARIATION WITH THE DIMENSIONLESS DROP SIZE. 

3H 


k 



WIT UK : 221 m OT UK : 415 * 

KM IT? : fl - COU MCLt 31 iff 

FWC S?*C : 44.1 .urns SLIMEPI : 58 1 m 
MITT lit : 21T news SLOPE : .751 

HU. CYCLE : 4 UKPLB : 1W 


OJ 

vO 


IIH* 

13 25 45 


l-L III 

123 news 



5 IIINETEKwcws' 177 


WIT LENS : 221 m KVI LENS : 415 » 

KM SEP : f I m COU. MILT: 38 in 

FT*C SMC : 14.2 «i crons SLIT SEP! : 58 I * 

WIST 111 : 213 mcpons SLOPE . : 751 

KIN. CYCLE : 4 SMELTS : 1888 



WIT LEW : 228 « 

IEM SEP : 5 2 ■ 

FB»: SPIC : 24 1 news 
WIST III : 213 «i crons 
KIN CYCLE : i 


KOI LEW : 415 w 
COU MILE: 38 iff 
SLIT SEPI : 58 I « 
SLOPE : .751 

SMI? LC : 1888 


W> CCVJH' 
354 


III 

123 ncron-. 



WIT LEW : 228 m 
K«N SEP : 13 » 

EMC SPIC : 18 7 Microns 
WIST DU : 213 Microns 
NIN CYCLE : 4 


KVI LENS : 415 

COU MCE: 38 4tj 
SLIT SLPI : 5§ I mm 
SLOPE : 751 

SIMPLES : 1888 


FRINGE SPACING VARIATIONS 


OF FOOi< QJALITY 







ORIGINAL PAGE ‘S 
OF POOR QUALITY 


TIME: 

15:26:25 



XMIT LENS : 90 mm 
BEAM SEP : 8 mm 
FRNG SPAC : 7.1 Microns 
UAIST DIA : 87 Microns 
MIN. CYCLE : 3 


i-NEH PARAMETERS 2-REPEAT RUN 


RCUR LENS : 495 mm 
COLL. ANGLE: 30 deg 
SLIT SEPR.: 50.8 mm 
SLOPE : .758 
SAMPLES 


3-LIST RESULTS 


MAX. 

504 


Die: 

23 Microns 

24 Microns 


4-END ? | 


MONODISPERSE DROP SIZE VARIATIONS. 


40 



*!UNS 2» » KW UT : rt • 

R »» U* : 111 - COu MU )* *< 

K tn k : 12 2 IICN’I iu: H*i : 25 < » 

Will 111 : 21) iiiw' (Iff! : TM 

«» CYC'J : I tiffin : m 






OKiCiJiNiAL PASS S3 
OF POOR QUALITY 



B1ANEIER r.icrors) 


VKI LENS : 220 rr 
BEAH SEP : 14 rr 
FRN-i SFfiL : 9.9 r: crons 
WA!$: P1A : 213 morons 
HIN. CYCLE : $ 


m LENS : 495 rr 
COLL ANGLE: 30 de« 
SLII SEEP..: 50 8 rr 
SLOPE : .758 
ShHPL.ES : 500 


BEAM ATTENUATION TESTS - EQUAL BEAM INTENSITIES. 



HAY. COUN 


DlAHETDi(Hicrons) 


m LENS : 220 nr 
BEAH SEP : 14 RR 
FRNG SPkC : 9.9 mcruns 
HAISI DIA : 213 mcrons 
HIN. CYCLE : 8 


RCUP LENS : 495 rr 
COLL. AO: 30 de y 
SLIT SET'S.: 50.8 rh 
SLO l I : .758 
SAKE*. PS : 501- 


BEAM ATTENUATION TESTS - ONE BEAM ATTENUATED. 


(90:10 RATIO) 





HAX.COUNT 

6714 


DIAHEIER(nicrons) 


XMIT LENS : 220 ah 
BEAN SEP : 18,6 hh 
FRNG SPAC : 7,4 Hi crons 
WAIST DIA : 213 Microns 
MIN, CYCLE : 10 

BACKSCATTER TEST 


RCUF LENS : 495 hh 
COLL, ANGLE: 150 deg 
SLIT SEPR,: 50,8 hh 
MAX, DIA. : 280 Microns 
SAMPLES : 823? 

■ OPAQUE DROPS. 


FIRST SURFACE REFLECTION. 



MAX', COUNT 
700 


DIAMETER(nicrons) 


XMIT LENS 
BEAM SEP 
FRNG SPAC 
WAIST DIA 
MIN.CVCLE 


220 HH 
5 HH 

27.8 Microns 
213 Hicrons 
3 


ROIR LENS : 495 hh 
COLL, ANGLE: 150 deg 
SLIT SEPR,: 50,8 hh 
MAX, DIA. : 357 nicrons 
SAMPLES : 1005 


BACKSCATTER TEST - TRANSPARENT DROPS, 
PREDICTION FOR INTERNAL REFLECTION. 


43 




original FACE 13 

OF POOR QUALITY 



0 40 80 120 160 200 240 280 


DROPLET DIAMETER, pm 

COMPARISON OF DROP SIZE DISTRIBUTIONS MEASURED WITH OTHER MEANS. 



0 20 40 60 80 100 120 

DROP DIAMETER, pm 


COMPARISONS OF DROP SIZE DISTRIBUTIONS MEASURED AT THREE SIZE RANGE SELECTIONS. 


44 


OWG'NM- PRSE <3 

or POOR QUALITY 


POTENTIAL CHARACTERISTICS OF THE METHOD 

• LINEAR RELATIONSHIP BETWEEN THE MEASURED PHASE ANGLE AND DROP 
SIZE 

• SIZE RANGE OF 30 OR GREATER AT A SINGLE OPTICAL SETTING 

• OVERALL SIZF RANGE OF 3 TO 2000 MICRONS 

• SIMULTANEOUS SIZE AND VELOCITY MEASUREMENTS 

• RELATIVE INSENSITIVITY TO BEAM OR LIGHT SCATTER ATTENUATION 

• HIGH SPATIAL RESOLUTION 

• OPERATION IS SIMILAR TO AN LDV 

• ADAPTABLE TO EXISTING LDV SYSTEMS 

• CAN DISTINGUISH BETWEEN GAS PHASE AND DROPLETS 

• REDUCED SENSITIVITY TO MISALIGNMENT 

• CAN PERFORM MEASUREMENTS INDEPENDENT OF REFRACTIVE INDEX 


SUMMARY 

EVALUATIONS WERE CONDUCTED ON A NEW SPRAY DROP SIZE AND 
VELOCITY MEASUREMENT METHOD 

• THEORETICAL ANALYSES OF THE LIGHT SCATTERING PHENOMENA 
WERE DEVELOPED 

• EXPERIMENTAL TESTING VERIFIED THE ANALYSES 

t MEASUREMENTS UNDER SIMULATED SPRAY CONDITIONS WERE 
CONDUCTED 

• COMPARISONS OF SPRAY MEASUREMENTS WITH OTHER METHODS 
SHOWlD GOOD AGREEMENT 

• SIGNAL PROCESSING METHODS WERE DEFINED AND TESTED 
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MEASUREMENT OF SPRAY COMBUSTION PROCESSES* 

C. E. Peters, E. F. Arman, J. 0. Hornkohl and W. M. Farmer 
Applied Physics Research Group 
The University of Tennessee Space Institute 


The objectives of this planned three-year investigation of spray combustion 
processes include the following: 

(1) Making measurements of noncombusting spray fields to aid in 
developing adequate computational models of the dynamic interaction 
between the droplets and the turbulent gas phase. 

(2) Making measurements of hydrocarbon-air spray flames to provide 
benchmark results that can be used in the verification and refinement 
of numerical models of the entire spray combustion process. 

(3) Developing techniques for making laser measurements in typical 
spray combustion environments, along with techniques for acquiring, 
processing, displaying, and interpreting the dat 

The freejet configuration selected for the experimental study is shown in 
figure 1. The central droplet-air jet is the exhaust from a plain-jet atomizer 
of the type investigated by Lorenzetto and Lefebvre (ref. 1). At a nominal air 
velocity of 100 m/s, the atomizer provides a Sauter mean droplet diameter of 
60-80 pm. The diameter of the centerbody is 5.4 cm, the diameter of the atomizer 
nozzle exit is 1.27 cm, and the diameter of the outer nozzle exit is 45.7 cm. 

For the combustion experiments, the recirculation zone formed downstream of the 
centerbody acts as a flameholder. 

Conventional pressure, temperature and flowrate instrumentation is used to 
monitor and control the experimental apparatus. The gas velocity field in the 
turbulent spray zone is measured with a laser velocimeter (LV) system, which yields 
two components of the instantaneous velocity vector of a small seed particle in 
the optical probe volume. A schematic diagram of the UTSl-developed LV system 
(ref. 2) is shown in figure 2. The relatively large droplets from the atomizer 
cannot be assumed to be in dynamic equilibrium with the gas phase; consequently, 
amplitude discrimination is used to eliminate those light-scattering events that 
result from large particles in the probe volume. (The diameter of the seed 
particles is only about 2 pm.) Using this LV technique to measure many light- 
scattering events for a given position of the probe volume, one can determine 
two components of the mean gas velocity vector as wej.1 as the variances and the 
covariance. 

The large-droplet size and number density, along with one component of the 
velocity, are measured with the Particle Sizing Interferometer (PSI) technique 
initially proposed by W. M. Farmer (ref. 3). A schematic diagram of the PSI system 
developed at UTSI is shown in figure 3. 


*NASA Grant NAG3-370 
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In the hydrocarbon spray combustion experiments, space and time-resolved 
measurements of temperature and species number density will be made with either 
laser-Raman or laser-fluorescence techniques. The equipment required for these 
measurements is already available at UTSI. These measurements will be made with 
the collaboration of Prof. J. W. L. Lewis. 

The first year of this research program was devoted to preparing the experi- 
mental equipment. A pictorial view of the new spray combustion facility is shown 
in figure 4. This facility was designed, the required components were procured 
and/or fabricated, and the facility was installed in one of the UTSI laboratories. 

The outer air flow is drawn from the atmosphere by a fan located in the downleg 
of the facility; a variable-speed DC fan-drive motor provides well controlled 
bellmouth exit speeds up to approximately 11 m/s. As shown in figure 4, turning 
vanes and screens are used to provide a nearly uniform flow at the bellmouth 
entrance. The atomizer air is supplied by the UTSI low-pressure compressor system, 
which provides up to 0.2 kg/s of air flow at pressures up to 10 atm. After passing 
through a pressure regulator, the atomizer air passes through a choked venturi, 
which is used to provide an accurate measurement of the flow rate. The atomizer 
liquid (either water or hydrocarbon fuel) is supplied from a nitrogen-pressurized 
tank located in a heated enclosure outside the laboratory. The pressure drop 
across a parallel array of capillary tubes is used to determine the liquid flow rate. 

Not shown in figure 4 is the optical table on which the various laser measure- 
ment systems are mounted. This table, which is located beneath the freejet test 
section, provides translation in three directions; a total axial travel of 
approximately 75 cm is available. 

At present (February 1984), all subsystems of the spray combustion facility 
are operational. All the conventional instrumentation systems have been calibrated 
and the microcomputer software for processing the signals has been developed. The 
LV system has been installed and aligned, and all other preparations have been made 
for commencing the water-droplet experiments. We plan to complete the LV and PSI 
measurements on the water-droplet system in mid-1984. At that time, the hydrocarbon 
combustion experiments will be initiated. 
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AUTOMATIC HOLOGRAPHIC DROPLET ANALYSIS 
FOR LIQUID FUEL SPRAYS* 


Alan C. Stanton, Gerald W. Stewart and H. John Caulfield 
Center for Chemical and Environmental Physics 
Aerodyne Research, Inc. 


Substantial attention has been directed toward both modeling and experi- 
mental characterization of liquid fuel sprays in order to develop a fundamental 
understanding of the fluid mechanical and chemical Interactions which govern the 
combustion process in these two-phase flows (ref. 1). Two critical parameters 
needed to characterize spray combustion are the size and velocity distributions 
of the fuel spray droplets. Several l&eer-based optical techniques, Including 
scattered intensity, laser Interferometry, and laser velocimetry, have been 
applied in recent years to measurements in fuel sprays (ref. 2). An alternative 
approach is holographic analysis of sprays. In which the three dimensional 
spatial distribution of droplets is recorded in a hologram for subsequent 
analysis (ref. 3). For such an approach to be practical, however, an automated 
system for analysis of the holograms is required. The goal of the present re- 
search program is to demonstrate an automated approach for droplet hologram anal- 
ysis and to identify critical requirements for adaptation of this technique to 
measurements in spray combustion systems. 

The basic scheme for automated holographic analysis Involves an optical sys- 
tem for reconstruction of the three-dimensional real image of the droplet field, 
a spatial scanning system to transport a digitizing x-y image sensor through the 
real Image, and processing algorithms for droplet recognition which establish the 
droplet sizes and positions. The hardware for this system is straightforward. 

In the present demonstration experiment, we are utilizing the expanded and colli- 
mated beam from a 5 mW helium-neon laser for hologram reconstruction, an Imaging 
lens for magnification of the real image field, and a video camera and digitizer 
providing 512-by-5l2 pixel resolution with R-bit digitization. A mechanical 
stage is used to scan the hologram in three-dimensional space, maintaining con- 
stant image magnification. A test droplet hologram provided by Dr. Michael 
Farmer of tne University of Tennessee Space Institute Is used for development and 
testing of the Image processing algorithms. 

For the case of fuel sprays, where the liquid droplets may be assumed spher- 
ical, the image processing problem is simplified in that there is no rotational 
degree of freedom, and special edge- tracking algorithns for tracing the boundary 
of complicated shapes are not required. The problem reduces to one of recogniz- 
ing droplets by comparing pixel Intensities against a threshold level which is 
selected on the basis of the average droplet and background intensities. The 
droplet location In three dimensional space is then established by finding the 
position of best focus in the image field. We are studying two approaches for 
finding the best focus. In the simpler approach, the droplet image areas are 


*Work performed under Contract No. NAS3-24094 with G. Neeley of NASA Lewis 
Research Center as program monitor. 
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determined at successive image planes by counting the number of pixels above 
threshold. The plane with the smallest image area is the best focal plane. The 
second approach, which is expected to be more sensitive, analyzes Intensity gra- 
dients by examining intensity differences between adjacent pixels at the image 
boundary. The image plane with the sharpest intensity gradients is chosen as the 
plane of best focus. 

An extension of this technique, which we anticipate will be explored in a 
future phase of the program, is the determination of droplet velocities and tra- 
jectories by the use of multiple exposure holograms (ref. 4). In this approach, 
a pulsed laser is used to create a hologram by two or more exposures at well- 
known intervals. The droplet image locations in the reconstructed hologram are 
determined using the boundary-defining algorithms described above. Finally, sta- 
tistical pattern recognition techniques would be used to associate droplet image 
locations at later times with their locations at earlier times. 
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AUTOMATED droplet h olog ram anai ysis 
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AND LOCATION 
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HOI OGRAM 



53 








ORIGINAL PAGE 19 
OF POOR QUALITY 



INTENSJ] Y St' ANN I NG 

COMPARE PI XI I INTENSITIES 
AGAINSI THRE SIKH 11 1 1 VII 
FOR DROPILT Rl COGNITION 

OBTAIN DROPI I I ARI AS BY 
COIIN I ING nilMHI R 01 PI XI I s 
AROVI IIIRISHOin 

REST FOCAI I’lANI IS PI AMF 
01 MINIMUM ARI A 


F DGF GRADIENTS 

DEFINE DROPIET EDGES BY 
COMPARISON OF INTENSITIES 
OF AD.IACEN1 PIXEIS 

BISI FOCAl PIANE IS PLANE 
Willi SHARPEST LPGE GRADIENTS 


Mill 1 I PI I IXPOSIIRI HOI OGItAPHY 


. Mill T I PI I IXPOSIIRI 1101 OGRAMS MAP DROPI 1 1 DISPI ACEMENTS IN VD SPACE 

. APPIY DROPI I I Rl C0GN1 1 ION AIGORITHMS TO FIND DROPLET 10CATI0NS 

. IISI SI AT ISI I CAI PA Ilf RN RECOGNITION TECHNIQUES TO ASSOCIATE DROPIET 
IMAGES Al I ARI IE R AND EATER TIMES 

. OBI AIN DROPI I I VI I OCI I II S FROM Ml ASIIRI D DISPI ACEMENTS AND KNOWN 
T I Ml INIIRVAIS 





N8 4 20532 


PREDICTIONS OF SPRAY COMBUSTION INTERACTIONS* 
J-S. Shuen, A.S.P. Solomon and G. M. Faeth 
The Pennsylvania State University 


Measurements were completed in dilute particle-laden jets, nonevaporating 
sprays and evaporating sprays injected into a still air environment. The flows 
are stationary, turbulent, axisymmetric and conform to the boundary layer 
approximations while having well-defined initial and boundary conditions — to 
facilitate use of the data for evaluation of analysis of the processes. Mean 
particle (drop) sizes were in the range 2-210 microns over the entire data base. 
The following measurements were made (as appropriate for the flow): mean and 

fluctuating phase velocities; mean particle mass flux; particle size; and mean 
gas-phase Reynolds stress, composition and temperature. Three models of the 
processes, typical of current practice, were evaluated using both existing and 
the new data: (1) a locally homogeneous flow (LHF) model, where interphase 

transport rates are assumed to be infinitely fast; (2) a deterministic 
separated flow ( DSF ) model where finite interphase transport rates were 
considered but effects of turbulence/drop interactions were ignored; and (3) a 
stochastic separated flow (SSF) model where effects of finite interphase 
transport rates and turbulent fluctuations were treated using random sampling 
for turbulence properties in conjunction with random-walk computations for 
particle (drop) motion. All three models used a k-e-g model for the continuous 
phase — which performed well in earlier studies of single-phase jets. The LHF 
and DSF models did not provide very satisfactory predictions over the present 
data base. In contrast, the SSF model generally provided good predictions and 
appears to be an attractive approach for treating nonlinear interphase 
transport processes in turbulent flows containing particles (drops). Current 
work is considering measurements and analysis of combusting sprays. 


INTRODUCTION 


The main objective of this investigation was to complete measurements of 
the structure of dilute particle-laden jets and sprays in order to support 
development of analysis of these processes. Test configurations involved 
injection into still air, yielding stationary, turbulent, axisymmetric flows, 
having well-defined initial and boundary conditions, which conform to the 
boundary layer approximations. The complexity of interphase transport was 
increased systematically, considering particle-laden jets, nonevaporating 
sprays, evaporating sprays and combusting sprays, in turn. In order to insure 
that appropriate measurements were made, model development was also 


*NASA Grant No. NAG3-190 with R. Tacina of Lewis Research Center serving as 
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undertaken — considering methods typical of current practice. The results of 
the study have application to the development of rational design methods for 
aircraft combustion chambers as well as other devices involving spray 
evaporation and combustion. 

The following description of the study is brief. Full details and a 
complete tabulation of data can be found in papers, reports and theses 
completed under the study (refs. 1-11). Studies of sprays and particle-laden 
jets have recently been reviewed by one of us (refs. 12 and 13 ) and this 
discussion will not be repeated here. Based on these reviews, and others cited 
therein, there is general agreement that a most pressing need for gaining a 
better understanding of spray processes is the creation of a well-defined set 
of measurements of the structure of these flows — motivating the present 
investigation. 

In the following, experimental and theoretical methods are described first 
of all. This is followed by a discussion of results already completed for 
particle-laden jets, nonevaporating sprays and evaporating sprays. The paper 
concludes with a brief discussion of current worx being undertaken in 
combusting sprays. 


EXPERIMENTAL METHODS 


Test Apparatus 

The same test configuration was used for the particle-laden jets and 
nonevaporating and evaporating sprays, cf . , figure 1. The injector was 
directed vertically downward within a screened enclosure. The injector and 
enclosure were traversed, since optical instrumentation was mounted rigidly. 

An exhausc system was used, to prevent recirculation of fine particles, 
however, its operation had negligible influence on flow properties at the 
measuring station. 

Test conditions are summarized in tables I and II for the particle-laden 
jets and the sprays, respectively. In each case, an air jet from the same 
injector was also tested, in order to establish test procedures by comparison 
with existing measurements. The particle-laden jets were essentially 
monodisperse with particle diameters in the range 79-207 microns while the 
sprays had Sauter mean diameters (SMD) in the range 30-87 microns. In general, 
measurements were confined to dilute regions of the flow where void fractions 
exceeded 99 %. All the flows were turbulent, with initial jet Reynolds numbers 
exceeding 10000. 


Instrumentation 

A single-channel laser-Doppler anemometer (LDA) was used to measure mean 
and fluctuating velocities of the continuous phase. High concentrations of 
seeding particles were used to avoid biasing due to flow particles (in the 
following, "particle" is used to represent "drop" unless the distinction is 
important). Data densities were high so that measured quantities could be 
time-averaged. The LDA was also used to obtain mean and fluctuating particle 
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velocities in the particle-laden jets — after terminating the flow of seeding 
particles. 

Double-flash photography was used to check LDA measurements of particle 
properties in particle-laden jets and to measure drop size and velocity 
correlations in the sprays. Drop size distributions obtained in this manner 
were corrected for depth-of-focus bias. Drop size distributions were also 
measured by slide impaction — corrected for Reynolds number bias. Finally, a 
Malvern particle sizer, which operates by Fraunhoffer diffraction of a laser 
beam, was used to monitor drop size distributions in the nonevaporating sprays, 
however, this instrument was ineffective in the evaporating sprays due to beam 
steering by gas-phase density gradients. 

Mean particle mass flux and liquid flux in the nonevaporating sprays were 
measured by isokinetic sampling at the mean gas velocity using a diverging 
probe. Isokinetic sampling with a heated probe, followed by analysis with a 
gas chromatograph, was used to measure mean total Freon-11 concentrations in 
the evaporating sprays. A shielded fine wire (25 micron diameter) thermocouple 
was used to measure mean gas-phase temperatures in the evaporating sprays. An 
impact plate was used to measure injector thrust. 

Measurements in the particle-laden jets were conducted for x/d = 1-50, 
where x is distance from the injector and d is the injector diameter. Adequate 
spatial resolution for the sprays could only be achieved for x/d = 50 with 
measurements extending to x/d = 500-600, due to the small exit diameter of the 
injector. For all test conditions, baseline calibration tests were conducted 
to establish predictions of transport rates to individual particles (drops). 


THEORETICAL METHODS 


General Description 

Three models were considered: (1) a locally homogeneous flow (LHF) model 

where interphase transport rates are assumed to be infinitely fast; (2) a 
deterministic separated flow (DSF) model where finite interphase transport 
rates were considered but effects of turbul ence/drop interactions were ignored; 
and (3) a stochastic separated flow (SSF) model where finite interphase 
transport rates and interactions between drops and turbulent fluctuations were 
treated using random sampling for turbulence properties in conjunction with 
random-walk computations for particle (drop) motion. A k-e-g model was used to 
find properties of the continuous phase for all three models, since this 
approach provided good structure predictions for single-phase jets (ref. 12). 
The test conditions correspond to steady, turbulent, axisymmetric boundary 
layer flows having low Mach numbers where effects of viscous dissipation of the 
mean flow and radiation are small. Other assumptions vary for the LHF, DSF and 
SSF models and will be treated separately in the following. 


LHF Model 

The LHF approximation implies that both phases have the same velocity and 
are in local thermodynamic equilibrium at each point in the flow — which is only 


57 


exact for im'i nitely-smal l particles. Therefore, the flow corresponds to a 
single-phase fluid with an unusual equation of state due to the presence of 
particles. The analysis employed Favre (mass)-averaged governing equations in 
conjunction with the conserved-scalar formalism of Bilger (ref. Hi). This 
procedure eliminates ad hoc neglect of density velocity correlations and 
effects of buoyancy in the governing equations for turbulence quantities. The 
relationship between scalar properties and mixture fraction, needed by this 
approach, followed past practice in this laboratory (ref. 13). The LHF model 
does not require detailed information concerning initial drop sizes and 
velocities; therefore, LHF calculations were begun at the injector exit. 


DSF Model 

Both separated flow models adopt the features of the LHF model for the gas 
phase. The dispersed phase was treated by solving Lagrangian equations of 
motion for the particles and then computing source terms for interphase 
transport which appear in the governing equations for the gas phase. This 
involves dividing the particles into a number of groups at the initial 
condition and then computing their subsequent motion. 

Void fractions in the region of computation always exceeded 99% : 
therefore, the dispersed phase volume, particle collisions and effects of 
adjacent particles on interphase transport rates could be ignored with little 
error. Ambient conditions for particles were taken to be local mean gas 
properties; therefore, effects of turbulent fluctuations on interphase 
transport, turbulent dispersion and turbulence modulation were ignored — typical 
of most current spray models (ref. 13). Effects of varying local ambient 
conditions, however, were considered. 


SSF Model 

The SSF model treats turbulence /drop interactions by computing drop 
trajectories as they move away from the injector and encounter a succession of 
turbulent edd ies--using Monte Carlo methods. Properties within each eddy are 
assumed to be uniform, but to change in a random fashion from one eddy to the 
next. Trajectory calculations are the same as the DSF model, except that 
instantaneous eddy properties replace mean gas properties. Eddy properties are 
found by making a random selection from the probability density functions (PDF) 
of velocity and mixture fraction. A drop is assumed to interact with an eddy 
as long as its relative displacement with respect to the eddy is less than a 
characteristic eddy size and its time of interaction is less than a 
characteristic eddy lifetime. All these parameters are directly found from the 
k-r-g computations. 


RESULTS AND DISCUSSION 


Particle-Laden Jets 


Initial 

measurements 


evaluation of the models was undertaken using existing 

in the literature (refs. 1,2). It was found that the LHF model 
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was only effective for flows containing tracer-sized particles, that the DSF 
model was not effective for any flow, and that the SSF model yielded 
encouraging predictions of flow structure and particle dispersion. The 
evaluation was not definitive, however, due to uncertainties in initial 
conditions for existing measurements. 

Initial conditions were fully defined for present tests with 
particle-laden jets --removing earlier limitations. Typical structure 
measurements and predictions are illustrated in figures 2-6 (u and v are axial 
and radial velocities, G is particle mass flux and subscripts o and c denote 
initial and centerline conditions). In figure 2, the agreement between 
predicted and measured velocities in the gas jet is good — establishing a 
baseline for the work. The LHF and SSF predictions are nearly identical =<rd 
are also in good agreement with measurements. Particle velocities, illustrated 
in figure 3, exhibit more significant effects of the model. The LHF approach 
underestimates particle velocities since effects of slip are ignored while the 
SSF model yields good results. Results for centerline particle mass flux, 
illustrated in figure h , are similar. 

Radial variations o r ga3 and particle properties are illustrated in 
figures 5 and 6. All gas-phase predictions are similar and are in reasonably 
good agreement with measurements, however, only the SSF model yields good 
predictions of mean and fluctuating particle properties. The DSF model 
underestimates particle spread rates, since effects of turbulent dispersion are 
ignored while the LHF model overestimates spread rates since slip is ignored. 

The flows were too dilute to test predictions of turbulence modulation by 
the SSF model. Sensitivity studies showed that predictions were most sensitive 
to initial particle properties. 


Nonevaporating Sprays 

A portion of the results for nonevaporating sprays is illustrated in 
figures 7-12. In this case, the LHF predictions are initiated at the injector 
exit while the separated flow model predictions are initiated at x/d = 50, 
where adequate initial conditions could be measured. Results for mean gas 
velocities and liquid fluxes along the centerline (figures 7 and 8) are 
similar’: the SSF model yields good predictions while the LHF model 

overestimates the rate of development of the flow. SMD variation along the 
axis is illustrated in figure 9. The SMD increases gradually along the axis, 
due to turbulent drop dispersion. The SSF model correctly predicts this trend 
while the DSF model yields an opposite trend since turbulent dispersion is 
ignored . 

Radial profiles of mean liquid flux are illustrated in figure 10. 

Turbulent dispersion of drops causes these flows to extend beyond r/x = 0.2, 
which is the usual boundary for a single-phase jet. The SSF model correctly 
predicts this trend — which is an encouraging finding — while the other models 
are unsatisfactory. 

Radial profiles of mean drop velocity for the case 1 spray and the SMD of 
both sprays are illustrated in figures 11 and 12. The correlation between drop 
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size and velocity is predicted reasonably well by the SSF model — as is the SMD 
variation. 

Turbulence modulation was again not an important factor in these sprays 
and could not be adequately evaluated. Specification of initial drop 
properties was found to have the greatest influence or predictions. 


Evaporating Sprays 

A portion of the findings for evaporating sprays is illustrated in figures 
13-16. Tends are qualitatively similar to results for nonevaporating sprays, 
however, differences between the models are somewhat reduced due to effects of 
drop evaporation. 


Combusting Sprays 

Current tests and analysis are extending results to combusting sprays. 
Initial tests are limited to dilute conditions where the flame is primarily 
fueled by gas (methane) with a monodi3perse stream of drops in coflow with the 
gas fuel. This system simplifies measurements and the presentation of results 
and allows greater concentration on the interesting fundamental problem of 
drop/combusting turbulent flow interactions. Subsequent work will consider 
flames fully fueled with polydisperse sprays. 


CONCLUSIONS 


Major conclusions are: (1) the SSF model provides a useful approach for 

treating turbulence/drop interactions in particle-laden flows and sprays with 
minimal empiricism and a capability to treat nonlinear effects; (2) the LHF 
model is a useful simplification for very well-atomized sprays (particle/drop 
size less than 10 urn) but effects of slip are important for most practical 
sprays; and (3) the DSF model was not effective for any of the flows examined 
here, since effects of turbulent particle dispersion and effects of turbulent 
fluctuations on interphase transport rates are important in most practical 
sprays. 

Specification of initial conditions is the most critical aspect of spray 
structure predictions. Effects of turbulence modulation were small for present 
flows, and current SSF model treatment of this phenomenon could not be 
evaluated decisively. Current work is extending the data base and the model 
evaluation to dilute combusting sprays. 
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ORIGINAL PAGE IS 
OF POOR QUALITY 

TABLE I. -SUMMARY OF PARTICLE-LADEN JET TEST CONDITIONS 


Flow 

Air Jet 


Dilute Particle 

-Laden 

Jets 

Case 

— 

1 

2 

3 

U 

SMD (microns) 

— 

79 

119 

119 

207 

Loading Ratio 

— 

0.20 

0.20 

0.66 

0.66 

Injector Exit Conditions: 6 






Velocity (m/a) 






Air 

32.1 

26.1 

29.9 

25.2 

25.3 

Particle 

— 

2k. 1 

29.2 

21 .9 

18.5 

Particle Mas3 Flux (kg/m-s) 

— 

6.1 

6.5 

18.9 

19.0 

a Injector exit diameter of 10.9 
b At flow centerline. 

mm; particle density 2620 kg/m’. 



TABLE 

II. -SUMMARY OF SPRAY TEST CONDITIONS 3 




Nonevaporating 

Evaporating 

Flow j 

Air Jet 

Sprays 

Sprays 

Case 

— 

1 

2 

1 

2 

Injected Fluid 

air 

air/oll b 

alr/oll b 

air/Freon-1 1 

alr/Freon-1 1 

Flow Rates, mg/s 






Gas 

338 

338 

216 

188 

120 

Liquid 

0 

600 

1900 

1950 

1899 

Loading Ratio 

0 

1.8 

6.5 

7.7 

15.8 

Jet Momentum, mN 

120 

137 

70 

106 

60 

Initial Velocity, m/s° 

175 

196 

**3 

65 

30 

Reynolds number 0 

26000 

30000 

29000 

91000 

36000 

SMD, microns 

— 

30 d 

87 d 

31 e 

58 e 

Spray Angle, deg. 

— 

30 

33 

27 

29 


a Spraylng systems air-atomizing Injector, 1.2 mm exit diameter. 
b Vacuum pump oil, Sargent-Weleh, Cat. No. 1907K25. 
c A3suming LHF with air viscosity '"or Reynolds number. 
d Measured with Malvern, Model 2200 Particle Sizer 0 x/d « 12.6. 
e Mea3ured by 3llde impaction at x/d - 50. 
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Figure 1. 

SKETCH OF THE EXPERIMENTAL 
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Figure 3. 

CENTERLINE MEAN PARTICLE VELOCITY 
IN PARTICLE-LADEN JETS 



Figure 2. 

CENTERLINE MEAN GAS-PHASE VELOCITY 
IN PARTICLE-LADEN JETS 



Figure U. 

CENTERLINE MEAN PARTICLE MASS FLUX 
TN PARTICLE-LADEN JETS 
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ORIGINAL FACE JS 



Figure 5. 

GAS-PHASE PROPERTIES AT x/d 40 
FOR CASE 1 PARTICLE-LADEN .1 El- 
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Figure 7 

CENTERLINE MEAN GAS-PHASE VELOCITIES 
FOR NONEVAPORATING SPRAYS 
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Figure 6 

PARTICLE PROPERTIES AT x/d - 40 
FOR CASE 1 PARTICLE-LADEN JET 
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Figure 8 

CENTERLINE MEAN LIQUID FLUX 
FOR NONEVAPORATING SPRAYS 
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Figure 9 
CENTERLINE SMD 
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Figure 11 

MEAN DROP VELOCITIES AT x/d = 250 
FOR CASE 1 NONEVAPORATING SPRAY 
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Figure 10 

RADIAL PROFILES OF MEAN LIQUID FLUX 
IN NONEVAPORATING SPRAYS 




OATA 

CASE 

THEORY (SSF) 

4 

O 

A 

1 

2 


ol 1 — 

1 

— 1 

1 l_ 

000 004 

008 

0 12 

016 020 


r/« 

Figure 

12 


SMD AT 

x/d 

- 250 

FOR 

NONEVAPORATING SPRAYS 


65 





ORIGINAL PAGE 19 
OF POOR QUALITY 



T T ITT? — *f 

DATA CASE theory 

O I DSF 

A l SSF 


I T 1 T T 


" «oJ 





_J I 1 L 1 _L_lJ I 

40 60 *0 IOO 200 

I/O 




400 600 600 


Figure 15 

CENTERLINE SMD FOR EVAPORATING SPRAYS 


Figure 13 
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Figure 1 

CENTERLINE MEAN FREON-11 
CONCENTRATION FOR EVAPORATING SPRAYS 
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Figure 16 

RADIAL PROFILES OF MEAN FREON-11 
CONCENTRATION FOR EVAPORATING SPRAYS 
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EFFECT OE LIQUID DROPLETS ON TURBULENCE 
STRUCTURE IN A ROUND GASEOUS JET 

J 

A.M. Mostafa and S.E. Eighobashi 
Mechanical Engineering Department 
University of California, Irvine 


Accurate prediction of spray combustion is extremely difficult due to the 
complex physical and chemical phenomena encountered in this two-phase process. The 
interaction between droplets and the turbulent fluid, turbulence effects on chemical 
reaction and heat transfer (and hence on droplet vaporization) are just a few 
examples of the complexity. In order to understand the nature of these 
interactions, coordinated experimental and theoretical studies need to be performed 
in a stepwise manner thus isolating the phenomenon to be investigated. A turbulent 
non-reacting gaseous jet laden with solid particles or evaporating droplets is a 
flow which allows the study of the interactions between the two-phases. 

The recent experiment of Modarress, et al. (1982, 1983), which was performed 
in parallel with the present work, provided a much needed data to help understand 
the behavior of two-phase turbulent jets and validate the theoretical models. 
Eighobashi and Abou-Arab ( 1983) reviewed existing turbulence models for two-phase 
flows and indicated that these models are based on ad hoc modifications of single- 
phase turbulence models. They developed a two-equation turbulence model for 
incompressible dilute two-phase flows which undergo no phase changes. 

In order to validate the proposed model, a turbulent axisymmetric gaseous jet 
laden with spherical uniform-size solid particles is studied by Eighobashi, et al. 
(1984). The predictions of the mean flow properties of the two-phases and the 
turbulence kinetic energy and shear stress of the carrier phase show good agreement 
wicn the experimental data. 
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Mostafa and Elghobashi (1981) extended the proposed model for steady 
incompressible two-phase flow including phase change. This model is tested for the 
flow of a turbulent axisvmmetric gaseous jet laden with multisize evaporating liquid 
droplets by Mostafa and Elghobashi (1983). To avoid the problem of density 
fluctuations of the carrier phase at this stage, only isothermal flow is considered 
and vaporization is assumed to be due to the vapor concentration gradient. The 
droplets are classified into finite size-groups. Each group is considered as a 
continuous phase interpenetrating and interacting with the carrier phase. Predicted 
results include distributions of the mean velocity, volume fractions of the 
different phases concentration of the evaporated material in the carrier phase, 
turbulence intensity and shear stress of the carrier phase, droplet diameter 
distribution and the jet spreading rate. The results are analyzed based on a 
qualitative comparison with the corresponding single phase jet flow. 

More validation testing of the model is needed via well-defined experiments. 
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OBJECTIVES 


ORIGINAL PAGE \3 
OF POOR QUALITY 


To predict the f1o« of i turbulent gaseous jet 
laden with multisize liquid droplets undergoing 
vaporization. 


To compare the effects, on the flow field, of 
vaporizing droplets uith those of single phase 
jet. 


ASSUMPTIONS 


Droplets in a given size-range constitute a 
dispersed phase - thus for k sizes there mill be 
(k ♦ 1) phases in the flow. 

Each phase behaves as a continuum 
(macroscop i ca 1 1 y ) . 

No collisions between droplets (dilute spray). 

Droplets remain spherical as they decrease in 
size. 

Velocity of vapor leaving droplet surface is 
equal to that of the droplet. 

Concentration gradient is the only driving force 
for evaporation. 

Drag relations of a solid sphere apply to liquid 
droplets (internal circulation effects on 
friction drag are counterbalanced by evaporation 
effects on pressure drag). 
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' Fig 1 Normalized mean velocities and volume fractions 

profiles at z/D * 20 and at X ■ 0.5 
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Fig 2 The turbulence intens ity distribution under 
different nass loading ratios at z/D * 20 


Fig 3 


The shear stress distribution under different mass 
loading ratios at z/D * 20 
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Fig 8 The axial distribution of the maximum turbulence 

intensity under different mass loading ratios Fig 9 The axial distribution of the maximum shear stress 

under dif ferent mass loading ratios 



Fig 10 Radial disc r lbur ioi.u of the local droplets diameters at different 
axial locations at it X >0.5 








Fig 11 The spreading rate under different T.ass loading ratios 


SUN* ARY 


• The proposed turbulence model (validated 
eipen mentally for solid particles) Is applied to 
predict the turbulent round jet laden with 
multisize vaporizing droplets. 

• Predicted changes In the flow field (lean and 
turbulent quantities) are significant and should 
be considered In liquid spray combustion 
calculations. 

A 

a Experiment to validate the proposed model mill be 

performed mhen funds become available. 
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COMBUSTION CHARACTERISTICS IN THE TRANSITION REGION 
OF LIQUID FUEL SPRAYS 

Nicholas P. Cemansky, Izak Namer, Robert J. Tidona and Hamid Sarv 

Drexel University 


There have been a number of studies involving the behavior of fuel sprays in 
the transition region which encompasses droplets in the 10-80 ym size range. In 
particular, the effect of droplet size on minimum ignition energy of a spray has 
been examined. Using swirl atomizers to generate fuel sprays, Ballal and Lefebvre 
(ref. 1) showed that ignition energies increase with increasing Sauter Mean Diame- 
ter (SMD). On the other hand, the experimental work of Chan and Polymeropoulos 
(ref. 2), in an 8-35 ym monodisperse droplet size range, showed an optimum drop- 
let diameter around 15 ym for which the ignition energy was minimum. However, 
these systems had limitations induced by droplet size distribution or variations 
in the flow velocities due to the requirements for monosized droplet formation. 
Thus, in November 1982, an investigation of the droplet ignition requirement was 
initiated with the following specific objectives: 

1) To map the minimum ignition energies for monodisperse fuel sprays 
in the transition region as a function of droplet size } number 
density, equivalence ratio, etc.; 

2) To isolate the optimum droplet diameter at which the ignition 
energy requirement is minimum; 

3) To study the effects of changing vaporization environment on the 
minimum ignition energy; 

4) To characterize the ignition requirements of polydisperse fuel 
sprays in terms of the droplet size distribution; and 

5) To determine the effects of droplet size and size distribution on 
flame speed. 

The experimental facility developed for this study consists of a Berglund-Liu 
Vibrating Orifice Monodisperse Aerosol Generator which produces monosized droplets 
with diameters having a standard deviation of less than 1% of the mean droplet 
diameter (ref. 3). Based on the vibrating orifice principle, the droplet size can 
be calculated from the operating conditions of the instrument (fuel feed rate and 
vibration frequency). After generation, the droplets are subjected to a flow of 
dispersion air to prevent coagulation. The air/fuel aerosol passes through a flow 
reducing section where dilution air is added to achieve the desired stoichiometry. 
The monodisperse fuel spray then flows through the test section, a square pyrex 
tube, which houses the ignition electrodes. Both electrodes are mounted on micro- 
meter traversing mechanisms to vary the spark gap. In its current configuration 
(Fig. 1), the experimental facility is operated with downward flow in order to 
prevent any accumulation of fuel within the combustor. 
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An ignition system has been developed to supply a capitance spark with indepen 
dently variable energy, energy density, and duration. The spark energies are deter 
mined by integrating (over time) the product of the voltage and current, recorded 
by a two-channel digital oscilloscope. 

Efforts in the first year have concentrated on the development of the experi- 
mental facility and performing exploratory measurements. Preliminary measurements 
have been made in the 10-74 pm monodisperse droplet size range using n-heptane as 
fuel. Minimum ignition energies are measured as a function of droplet size, equi- 
valence ratio, and flow velocity. Prior to these measurements, the optimum spark 
gap and duration corresponding to minimum ignition energies are determined for each 
operating condition. Further measurements using different fuels including alcohols 
will be made soon. 


REFERENCES 


1. Ballal, D.R., and A.H. Lefebvre, "Ignition and Flame Quenching of Flowing 
Heterogeneous Fuel-Air Mixtures. " Combustion and Flame 35 , 155, 1979. 

2. Chan, K.K., and C.E. Polymeropoulos , "An Experimental Investigation of the 
Minimum Ignition Energy of Monodisperse Sprays," Paper No. ESSCI 81-21, 

Eastern States Section, The Combustion Institute Meeting, November, 1981. 

3. Berglund, R.N., and B.Y.H. Liu, "Generation of Monodisperse Aerosol Standards, 
Environmental Science and Technology 7, 2, 147, 1973. 
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overall objectives 


BACKGROUND 


Study the effects of droplet size 
and sire distribution on 

• Mini mum Ign i t i on Ener g i es 

• Flammability Limits 
e F I ame Speeds 


e Spray Burning Widespread 

- 30-50% of Total Energy Consumption 

• Modes of Spray Combustion 

- Pr em i xed 

- D i f f us i ve 

- Tr ans • t i on 
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• Transition Region Effects 

- Increasea Flam® Speeds 

- Reduced N0 V Emissions 

- Br oadenea F I ammoD i lity Limits 

- Lower Ignition Requirements 


OBJECTIVES OF IGNITION STUDIES 


• Map minimum ignition energies for mono- 
disperse fuel sprays in transition region 

e Determine optimum droplet size corresponding 
to minimum ignition energy 

e Study pr evapor i za t i on effects on minimum 
i gn i t i on ener gy 

• Examine size distribution effects in po I y- 
(Jisperse sprays on ignition requirements 
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Figure I Schematic of the Spray Burner Facility 


EXPERIMENTAL VARIABLES 


SUMMARY 


• Droplet Size (10-74 jjm) 

• Number Density 


• A monodisperse aerosol generator has been modi- 
fied to study i gn i t i on r equ i remen ts, f I ammab i I i ty 
limits, and flame speeds in the transition region 


e Equivalence Ratio 


• An ignition system has been developed and tested 


e Extent of Prevapor izat ion 

• F I ou Ve I oc i ty 

• Fue I Type 

• Oxidizer Composition 


• An optical drop sizing system has been designed 
and fabrication is near I y complete 

• Prei minary measurements of droplet size effects 
on the minimum ignition energy for n-heptane 
sprays have been performed 

e Parametric studies of droplet size effects on 
minimum igni t ion energies of various fuels 
in.luaing a I choho I s are in progress 
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IURBULENT SWIRLING COMBUSTION* 

Holger T. Sommer, Raymond J. Mosula and Ellice Seiden 
Department of Mechanical Engineering 
Carnegie-Mel Ion University 


EXPER I MEN 1 AL EFFORT 
Turbine Combustion Facility 

An experimental facility constructed to simulate swirling gas turbine combustion 
has been subjected to extensive redesign and modification. The program undertaken 
over the course of the past year consisted essentially of two phases: (1) an assess- 

ment of the existing test facility, and (2) a major redesign and reconstruction of 
those components deemed in need of alteration. As a result of this extensive effort, 
the test facility at Carnegie-Mellon University has achieved a higher degree of pre- 
cision and reproducibility. Furthermore, it is believed that those undesirable 
extraneous effects which plague any fluid mechanics experiment have been substan- 
tially reduced. Hence, it is at this time that accurate and reproducible flow 
measurements can begin. 


Measurements and Results 

Measurements of the non-reacting turbulent flow produced by two confined, co- 
and/or counter-swirling jets have been obtained by means of a two-color laser Doppler 
velocimeter. These results are compared with earlier experiments by Vu and Gouldin 
(ref. 1), and the numerical predictions of Ramos (reii. 2 and 3). It is shown that 
under both swirl conditions, a closed recirculation zone is created at the combustor 
center line. This zone is characterized by the presence of a one-cell toroidal vor- 
tex, low tangential velocities, high turbulent intensity, and large dissipation 
rates of the kinetic energy of turbulence. 


THEORETICAL EFFORT 

Numerical Modeling of Turbulent Swirling Flows 

The major thrust of the theoretical effort has been the adaptation of previously 
developed computer codes by Lilley (ref. 4) and Ramos (ref. 5) to the above described 
experiment. Both codes are based on the k-e turbulence model, while Ramos addi- 
tionally utilizes a k-S. turbulence model. 


★ 

This project was supported by the NASA Lewis Research Center under contract number 
NAG3-401 . 
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in his numerical solution, Ramos (ref. 2) includes streamline curvature effects. 
The streamline curvature effect is being improved and a correction factor for the 
k-e model will be developed during the remaining funding period. Because the codes 
utilize different numerical solution techniques to solve the finite difference equa- 
tions, comparison between numerical results and measurements will be used to evaluate 
the quality of both techniques. 

Special care will be taken to model the recirculation zone casued by the pre- 
sence of co- and counter-swirling flows. It will be possible, due to clearly defined 
experimental boundary and initial conditions, to investigate the influence of tur- 
bulence models on complex reacting flows. These conditions have been found to be 
critical in numerical studies. 


REFERENCES 
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ERROR REDUCTION PROGRAM* 

A PROGRESS REPORT 

Saadat A. Syed 

Engineering Division; Pratt & Whitney Aircraft Group 


The combustion chamber of the aircraft gas turbine has to be designed to satisfy 
the stringent dual engine requirements of reduced fuel consumption and increased 
durability. Combustor development based on conventional methods to meet these 
requirements is proving to be increasingly costly and time consuming. A mathe- 
matical model of the combustor would be very helpful in reducing cost and time of 
the design cycle as it would permit most of the work to be done on a computer. Such 
models are already being used in a limited manner in the industry. More widespread 
use of these models will be accelerated by the removal of some of the known 
deficiencies in the present codes. One of these deficiencies is the numerical 
error, or numerical diffusion, that can be introduced into the calculations under 
certain conditions. This numerical error is due to the finite differencing scheme, 
usually upwind or hybrid, used in these codes. 

The objective of the Error Reduction Program is to evaluate available finite 
difference schemes for minimum numerical diffusion, and to incorporate the best 
available scheme into a three-dimensional combustor performance code. 

As a first step toward this end, five schemes (refs. 1 through 5) were selected 
for initial evaluation. Some of these schemes were combined with bounding schemes 
to eliminate physically unrealistic undershoots and overshoots. Two basic criteria 
were used to judge these schemes: they should be conservative, and should produce 
solutions that exhibit no extraneous maxima or minima (boundedness property). The 
accuracy of the schemes was evaluated by performing the truncation error analysis, 
and running one- and two-dimensional test cases and comparing the calculated 
solutions against the exact solutions. Based on this evaluation, two schemes were 
selected: QUDS, Quadratic Upstream Differencing Scheme, and BSUDS2, Bounded Skew 
Upstream Differencing Scheme Two. The first scheme was proposed by Leonard (ref. 
1), and the second scheme is the one proposed by Raithby (ref. 2), which is bounded 
by a new bounding scheme. 

The selected two schemes were coded into a two-dimensional computer code, 2D-TEACH, 
and their accuracy and stability were evaluated by running several test cases (refs. 
6 through 8). It was found that BSUDS2 was more stable than QUDS. It was also 
found that the accuracy of both schemes is dependent on the angle that the 
streamlines make with the mesh, QUDS being more accurate at smaller angles and 
BSUDS2 being more accurate at larger angles. However, both schemes, at all angles, 
were more accurate than the existing hybrid scheme. BSUDS2 was selected to be 
extended into three dimensions, primarily because it was more stable. This scheme 
is currently being incorporated into a three- dimensional code, 3D-TEACH. 


* This work was conducted under NASA Contract NAS3-23686. 
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NEED FOR ERROR REDUCTION 


OBJECTIVE 



Convection and diffusion of a scalar step in constant property 
uniform inclined flow. Pe B Pe y 60 


To identify and inco'porate best available 
finite difference error reduction scheme 
into a 3D combustor performance code. 


PROGRAM SCHEDULE 
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DISCRETIZATION SCHEMES 
EVALUATED 


• Agarwal Differencing Scheme (ADS) 

• Quadratic Upwind Differencing Scheme (QUDS) 

• Skew Upwind Differencing Scheme (SUDS) 

• Spline/Hermetian Schemes 

• Cubic Spline Scheme (CSS) 

• Glass and Rodi Hermetian Scheme (GRHS) 

• Flux Blending Schemes 

• Bounded One 

• Bounded Two 


SCHEMES SELECTED 


• QUDS — Quadratic Upwind Differencing Scheme 

• BSUDS2 — Bounded Skew Upwind Differencing 

Scheme Two 


WHAT IS BOUNDEDNESS 

In the absence of sources the value at node P should 
be bounded by the values at surrounding nodes. 

N 



S 
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QUADRATIC UPWIND 
DIFFERENCING SCHEME 
(QUDS) 

• Finite volume method 

• Upwind biased quadratic interpolation for convection terms 

• Central differencing for diffusion terms 

• Employs extended nine point modecule 

$w '^^WW * + .u >o 

= - + Vi^p + Vs^yy.u < o 

• Scheme is conservative 

• Coefficients can become negative 

• Solution can become unbounded 



BOUNDED SKEW UPWIND 
DIFFERENCING SHCEME TWO 
(BSUDS2) 

• Finite volume method 

• Central differencing for diffusion terms at all Peclet 

numbers 

• Central differencing for convection terms for absolute 

Peclet numbers less than two 

• Blending of upwind and skewed 

upwind differencing for Peclet 
numbers greater than two 

• Employs compact nine point 

molecule 

• Scheme is conservative 

• Produces bounded solutions 
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FIRST LAMINAR FLOW TEST CASE 

QUDS more accurate over all 


XR 

h 


7r 


Experimental (F. DURST) 



I 1 1 1 

10 100 1.000 10.000 

Number of nodes 


SECOND LAMINAR FLOW TEST CASE 



Flow in a sudden expansion with swirl; Re = 450. swirl no. = 0.66. 
at 45-degree vane angle, contraction ratio = 9.0 
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STAGNATION STREAMLINES FOR SECOND 
LAMINAR FLOW TEST CASE ON 40 x 20 MESH 



STAGNATION STREAMLINES FOR 
SECOND LAMINAR FLOW TEST CASE 
ON 78 x 40 MESH 
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FIRST TURBULENT TEST CASE 
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1275 

Axial distance, mm 


2540 


Inlet boundary conditions 


Mean velocity 

Temperature 

Pressure 

Boundary layer thickness 
Turbulence intensity 


60 1 1 ft sec 
59 F 
14 7 psia 
0 33 inch 
0 003 


(18 32 m/sl 
(41 2°C) 

1101 35 kN/m 2 l 
(8 4 mm) 


CALCULATED REATTACHMENT POINTS 
(MEASURED 7 STEP HEIGHTS) 


Case 

HYBRID 

BSUDS2 

QUDS 

Dumb (15x17) 

2.4 

2.25 

2.28 

Coarse (26x29) 

5.2 

5.4 

5.5 

Fine (50x56) 

5.7 

5.9 

Unstable 

Fine adjusted 
(74x53) 

5.8 

L 

5.8 

1 

Unstable 
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SECOND TURBULENT TEST CASE 



SECOND TURBULENT TEST CASE 
WITH SWIRL 

Inlet boundary conditions; 

• Measurements at x = 5mm were used 

• This plane was considered as inlet plane 


SECOND TURBULENT TEST CASE 
WITH SWIRL 
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SECOND TURBULENT TEST CASE WITH SWIRL 
COMPARISON OF 
AXIAL VELOCITY PROFILES 


Coarse grid (40x35) 

HYBRID 

QUDS 

BSUDS2 


x * 51mm oat 

x * 1 52mm » a 


y 

1 

i l 

V* 3 

3 


7 -" , v‘ " - 

■ ■ ■ ‘ » ‘ * *■ » ••• 

• * « 1 . 00 


V O 610 0,# ° °* VO §10 


x = 305mm 



V 0 610 


QUDS seems more accurate 


No difference in accuracy for fine grid (47x43) 


SECOND TURBULENT TEST CASE WITH SWIRL 
COMPARISON OF 
TANGENTIAL VELOCITY PROFILES 

Coarse grid (40x35) 

HYBRID 

QUDS 

BSUDS2 



x - 1 52mm 



ooe 


VO §10 




Y 0 §10 


QUDS seems more accurate 
No difference in accuracy for fine grid (47x43) 
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SINGLE CELL CALCULATIONS 
OF CONVECTION OF NORMAL DISTRIBUTION 


<£(x) exp I '/» CIA* I * I 
[V x/A.A* A/<t) 



\ 



UDS 


0 75 i i 

0 0 7.5 15 0 22 5 30 0 37 5 45 0 

6 . degree 

Influence of the flow angle on the error 
for various difference schemes 


OBSERVATIONS 

• QUDS unstable — needs a better solver 

• No significant differences between HYBRID, QUDS 

and BSUDS2 for intelligent fine grids 
(1000 1500 nodes) 

• QUDS more accurate than BSUDS2 most of the time 

• BSUDS2 more accurate for second laminar test case — 
very strong streamline curvature 

• Flows with swirl more sensitive to difference schemes 
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There have been a variety of methods used to solve the Havier-Stokes 
equations. Most schemes which solve the steady state form of there equations 

employ the SIMPLE (Semi-Implicit Method for Pressure-Linked Equations) method 
developed by Patankar and Spalding in 1972 (ref.l). Although this scheme has 
proven to be quite effective, its convergence rate can be improved. This 
investigation focuses on two approaches which accelerate the solution of the 
steady state Navier-Stokes equations. 

The SIMPLER algorithm (ref. 2), a revised version of SIMPLE, provides a 
more accurate pressure field for each iteration through the momentum 
equations, thereby speeding convergence. PISO (Pressure Implicit Split Operator, 
ref. 3) performs a secondary correction of the velocity and pressure fields 
(after the typical pressure correction) which enhances convergence. Both 
schemes account for terras neglected in the SIMPLE approach, but do so in 
slightly different ways. 

A series of calculations of two-dimensional driven cavity flow and flow 
over a step were made to examine the effect of geometry on the performance of 
these schemes. Computations were carried out on a series of progressively finer 
grids. The effect of relaxation number on convergence rate was analyzed, using 
results from SIMPLE as criteria for performance correlation. 

Preliminary results indicate the following: 

1. The improved schemes promoted convergence by up to sixty percent for the 

driven cavity and forty percent for flow over a step. 

2. For the driven cavity problem, the efficiency of PISO and SIMPLER 

increased as the number of nodes increased. 

3. To ensure faster convergence, higher relaxation numbers must be applied. 
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NUMERICAL MODELING OF TURBULENT FLOW 


R. W. Claus 

NASA Lewis Research Center 


Three dimensional combustor calculations are currently 
stretching the computer hardware capabilities and the 
computing budgets of gas turbine manufacturers. One of the 
main reasons for this relates to the large number of complex 
physical processes occuring in the combustor. Airflow, fuel 
spray, reaction kinetics, flame radiation, and , not the 
least of which, turbulence must be modeled and the related 
differential equations solved. Obtaining accuracte solutions 
to these modeled equations entails another difficulty. 
Current combustor codes are, generally, based on the SIMPLE 
algorithem developed by Patankar and Spaulding, ref. 1. One 
of the key features of this algorithem is the use of hybrid 
differencing to approximate The convective terms in the 
governing equations. In most practical calculations this 
results in the use of first order accurate upwind 
differencing for most of the flow field. The overly 
dissipative solutions obtained using this scheme can mask 
important flow field features, ref. 2. Ideally, this can be 
overcome through the use of additional grid points. In three 
dimensional calculations , however, this becomes impractical 
since the computational work involved increases greatly as 
grid points are added. 


Other discussions in this conference will address methods 
to improve the accuracy of combustor flow field calculations 
and methods to speed the convergence of the modeled 
equations. This talk will focus on aspects of merging these 
two new technologies. The improved accuracy discretization 
schemes have a negative impact on the speed of convergence 
of the modeled equations that the improved solution 
algorithems may not overcome. A description of the causes of 
this problem and potential solutions will be examined. 


REFERENCES 
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ADVANCED NUMERICS FOR MULTI-DIMENSIONAL FLUID FLOW CALCULATIONS 

S. P. Vanka 

Argonne National Laboratory 


In recent years, there has been a growing interest in the development and use of 
mathematical models for the simulation of fluid flow, heat transfer and combustion 
processes in engineering equipment (ref. 1). The equations representing the multi- 
dimensional transport of mass, momenta and species are numerically solved by finite- 
difference or finite-element techniques. However despite the multitude of differ- 
encing schemes and solution algorithms, and the advancement of computing power, the 
calculation of multi-dimensional flows, especially three-dimensional flows, remains a 
mammoth task. The following discussion is concerned with the author's recent work on 
the construction of accurate discretization schemes for the partial derivatives, and 
the efficient solution of the set of nonlinear algebraic equations resulting after 
discretization. The present work has been jointly supported by the Ramjet Engine 
Division of the Wright Patterson Air Force Base, Ohio, and the NASA Lewis Research 
Center. 


FINITE-DIFFERENCING 


An efficient finite-difference scheme must represent the differential equations 
to a high accuracy and must he stable. In spite of the many formally higher order 
schemes that have been proposed, practical calculations have bee limited to the use 
of simple first order schemes. The reasons for this are that either the proposed 
schemes are unstable in general circumstances or they have limited range of formal 
accuracy. The complexities of the construction of accurate schemes stem from the 
following flow characteristics. 


a) 

High grid Peclet number 



b) 

Stream lines skewed with 

grid 

lines 

c) 

Pressure gradients 



d) 

Nonlinear source terms. 

also 

containing derivatives. 


In designing any generally useful finite-difference calculation procedure, it is 
necessary to consider all of the above factors carefully. 

The present work of the author is concerned with the Integration of several 
individual ideas into a general purpose algorithm. Two promising schemes that are 
being considered are the influence scheme approach of Chen et al. (ref. 2) and of 
Raithby et al. (ref. 3) and the compact differencing approach of Kriess (ref. 4). In 
the influence scheme approach, the differential operator is analytically Integrated 
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over two (or three) dimensional subdomains after linearizing It with previous Iterate 
values. A profile is assumed for variations along the boundaries. This results in 
an algebraic linear relation of the form 
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Vi + Vi 
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where $ -j are neighbor values and Oj represent the neighbor source terras. Cj and Dj 
represent coefficients. The advantages of the above approach are that It fully 
considers the mult ^-dimensional effects and is stable at high Peclet numbers. 

In a compact differencing scheme the derivatives are represented to a fourth 
order accuracy but with a tridiagonal structure. The tridlagonal structure Is 
advantage for the solution of equations and for the treatment of boundary 

conditions. The compact differencing scheme however has been found In our 
preliminary work to be unstable (and inaccurate) at the high Peclet numbers of 
current Interest. However the concept can still be used to represent the various 
other derivatives that occur In the source terms and the pressure gradients. The 
above two approaches are currentl, being pursued. 



SOLUTION ALGORITHM 


The solution of the nonlinear algebraic equations obtained after discretization 
Is a mammoth task. The required computer times are large and convergence is not 
always assured. Current work in this direction has led to a fully-coupled solution 
of the nonlinear equations of momentum and continuity. The coupled set of nonlinear 
equations are solved with a Newton-Raphson method and efficient sparse matrix 
techniques. The turbulence equations are solved decoupled from the momentum and 
continuity equations. The equations are also arranged as blocks which are then 
preordered to reduce the computer storage and time. This algorithm has been extended 
to reacting flows and successful calculations have been made of Isothermal turbulent 
flows and turbulent confined diffusion flames. The algorithm has been observed to be 
rapidly convergent and Insensitive to variations In grid aspect ratio, flow Reynolds 
number and the number of finite-difference nodes. 


Results of Some Calculations 

Several calculations of laminar and turbulent recirculating flows have been made 
with different finite-difference grids. These have been documented In detail in 
Vanka (refs. 5, 6 and 7). In all the cases tested, rapid diminution of the residuals 
in the equations has been observed. Calculations have been made with finite- 
difference grids up to (80 x 95) size for turbulent flow in a sudden expansion. In 
the present report two tables are presented which demonstrate generally the observed 
rates of convergence. The first one Is for a turbulent flow In a sudden expansion 
for the experiments of Craig et al. (ref. 8), with a (80 x 95) grid. The second one 
shows the rate of convergence for a confined diffusion flame for a configuration 
studied by Lockwood et al. (ref. 9). The latter calculations employed a (32 x 27) 
grid. 
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Table 1. Convergence Rate for Craig's Sudden-expansion Flow, 80 x 95 Grid 


Iter . 

6u/ui n 

6v/ui n 

6 P /u in 

6 k 

6e 

1 

1.391E-1 

2.904E-2 

2.766E-1 

5.846E-5 

8.858E-5 

5 

2.879E-2 

1.059E-2 

1.155E-2 

4.355E-5 

3.344E-5 

10 

7.300E-2 

5.500E-2 

1.171E-1 

4.584E-5 

3.336E-4 

15 

2.717E-3 

3.224E-4 

8.086E-4 

3.682E-5 

5.972E-5 

20 

1.508E-5 

1.360E-5 

2.597E-5 

4.246E-6 

1.970E-6 


Table 

2 . Convergence Rate 
(Experiments of 

for Confined Turbulent Diffusion 
Lockwood et al.) with 32 x 27 grid 

Flame Calculations 

Iter. 
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3.65E-5 

1.54E-5 
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GRID FLEXIBILITY AND PATCHING TECHNIQUES 

T. G. Keith, L. W. Smith, C. N. Yung, 

S. H. Barthelson, and K. J. DeWItt 

The University of Toledo 

The numerical determination of combustor flowflelds Is of great value to 
the combustor designer. An a priori knowledge of the flow behavior can speed 
the combustor design process and reduce the number of experimental test rigs 
required to arrive at an optimal design. Even 2-D steady Incompressible Iso- 
thermal flow predictions are of use; many codes of this kind are available 
[1], each employing different techniques to surmount the difficulties arising 
from the nonlinearity of the governing equations and from typically Irregular 
combustor geometries. We will look at mapping techniques (algebraic and ellip- 
tic PDE), and at adaptive grid methods (both multi-grid and grid embedding) as 
applied to axlsymmetrlc combustors. 


ALGEBRAIC GRID GENERATION IN THE STARPIC CODE 

The solution of the Turbulent Reynolds Equations for arbitrary geometries 
can be handled In several ways [1]. One popular technique Is to represent a 
boundary with discrete 90° steps. At lernat 1 vely , the geometry can be mapped 
Into a square solution plane via a variety of available transformations . The 
simplest transformation Is an algebraic one as shown In Figure 2. 

On transforming the governing equations, new chain rule terms are Intro- 
duced, each of which contains a dh/dx term which comes Into play only at axial 
locations where there Is a sloping or curved wall segment, thereby mathemati- 
cally affecting expansion and recirculation In the flow. The complexity of 
extra terms Is offset by the ability of the technique to handle general bound- 
ary shapes, as well as the simplification of boundary conditions In the mapped 
plane. For these reasons algebraic mapping shows great promise. Figures 1 
through 6 show the effects of algebraic mapping on the STARPIC code. 


ELLIPTIC PDE GRID GENERATION FOR A BIFURCATED 
COMBUSTOR INLET DIFFUSER 

The purpose of a combustor Inlet diffuser Is to convert kinetic energy to 
pressure; a reliable prediction method for diffuser flows will lead to the 
efficacious design of those diffusers. In this work body fitted coordinate 
transformations are employed In the solution of turbulent flows in a bifur- 
cated combustor Inlet diffuser. The work Is aimed at comparing the numerical 
solution for different kinds of transformations with experimental results. 
Furthermore, different numerical methods will also be employed for comparison. 
Figures 7 through 12 Illustrate this work. 


I 


MULTIGRID AND GRID EMBEDDING TECHNIQUE FOR 
TURBULENT FLOWS IN COMBUSTOR 

A multigrid method [2] applied to axlsymmetrlc turbulent flows Is useful 
In speeding convergence of numerical schemes. This method reduces global 
errors by economical relaxation of errors on coarse meshes Instead of by labor 
Intensive relaxation on one fine mesh. The multigrid equations, a summary of 
the multigrid algorithm, and relationships between coarse and fine meshes for 
a staggered mesh are depicted In Figures 13 through 17. 

Grid embedding [3,4] similar to the multigrid method, uses a fine mesh In 
steep gradient regions to refine the solution locally, where needed. It elim- 
inates the wasteful process of using a fine grid globally. Including In regions 
where It Is not needed. Grid embedding Is depicted tor a staggered mesh. As 
Indicated In the last figure, this method Is still under development. Like the 
multigrid method, It will steadily be extended to more arbitrary combustor 
geometries . 
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ADAPTED FROM 
STARPIC 


CO NTRASTING STARPIC AND STARMAP 

ADVANTAGES OF MAPPING: 

ANY SINGLE VALUED FUNCTION VALID FOR h(»). 

SOLUTION OCCURS ON A SQUARE MESH IN THE 5 - n PLANE. 
BOUNDARY CONDITIONS ARE SIMPLIFIED. 


DISADVANTAGES: 

ADDITIONAL TERMS ARISE VIA TRANSFORMATION. 
HANDLING OF "WALL FUNCTIONS" NOT OBVIOUS. 
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STATUS 

NEARING OPERATION OF COMPLETELY MAPPED CODE. 

FOLLOWING FIRST SUCCESSFUL RUN. CODE WILL BE 
TESTED WITH A VARIETY OF GEOMETRIES IN TWO 
PHASES: 

PHASE (I) VALIDATION OF STARMAP IN 

DIRECT COMPARISON WITH STARPIC 
RESULTS: 

PHASE (ID GENERALIZATION TO OTHER GEOMETRIES 
NOT PRESENTLY HANDLED BY STARPIC. 
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STATUS 

TRANSFORMATION OF THE GEOMETRY HAS BEEN ACCOMPLISHED. 
TRANSFORMED GOVERNING EQUATIONS HAVE BEEN DISCRETIZED. 
COMPUTER CODES TO SOLVE TRANSFORMED EQUATIONS ARE 
BEING PREPARED. 
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STATUS 

MULTIGRID METHOD IS BEING APPLIED TO THE STARPIC 
CODE 


FUTURE DEVELOPMENTS 

REFINE EXISTING MULT I GRID METHOD 

DEVELOP EMBEDDED GRID METHOD 

EXTEND BOTH METHODS TO MORE ARBITRARY GEOMETRIES 
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MASS AND MOMENTUM TRANSPORT EXPERIMENTS WITH SWIRLING FLOW 

Bruce V. Johnson and Richard Roback 
United Technologies Research Center 


An experimental study of mixing downstream of axial and swirling coaxial jets is 
being conducted to obtain data for the evaluation and improvement of turbulent trans- 
port models currently employed in a variety of computational procedures used through- 
out the propulsion community. The axial coaxial jet study was completed under Phase 
I and is reported in Ref. 1. The swirling coaxial jet study, which is the subject of 
this paper, was conducted under Phase II of the contract and is reported in Ref. 2. 

A TEACH code was acquired, checked out for several test cases, and is reported in 
Ref. 3. A study to measure length scales and to obtain a limited number of measure- 
ments with a blunt trailing edge inlet is being conducted under Phase III of the con- 
tract . 

Qualitative and quantitative studies were conducted of the flow downstream of 
swirling coaxial jets discharging into an expanded duct. The ratio of annular jet 
diameter and duct diameter to the inner jet diameter were approximately 2 and 4, re- 
spectively. The inner jet peak axial velocity was approximately one-half the annular 
jet peak axial velocity and the mean swirl angle in the annular stream was approxi- 
mately 30 degrees. Results from the studies were related to the five shear regions 
within the duct: (1) the wake region downstream of the inlet, (2) the shear layer 
between the jets, (3) the annular recirculation region, (4) the reattachment region, 
and (5) the centerline recirculation region (Fig. 1). 

A flow visualization study was conducted using fluorescence dye as a trace 
material and high-speed motion pictures to record the dye patterns in selected r-z 
and r-6 planes. The results of the flow visualization study are summarized with the 
following observations (Figs. 2 and 3): 

1. The high intensity turbulent eddies in the shear layers were not axisymmetric or 
periodic. The large scale waves and eddies apeared to have a range of wave- 
lengths . 

2. Downstream of the centerline recirculation region, the flow initially had large 
lobes which evolved into a vortex swirl pattern. There was little radial mixing 
downstream of the recirculation cell. 

3. Two major mixing regions were observed: (1) at the interface between the inner 
jet and the centerline recirculation zone, and (2) at the interface between the 
inner jet and annular jet streams. 

4. Mixing at the interface of the inner stream and the recirculation cell diluted 
the inner jet concentration and the resulting mixture of fluid from the 
recirculation cell and the inner jet was entrained by the swirling annular 
stream . 

A detailed map of the velocity, concentration , mass turbulent transport rate 
and momentum turbulent transport rate distributions within the test section was ob- 
tained to provide data for the evaluation and improvement of turbulent transport 
models. Data sets of velocity component pairs were obtained simultaneously to deter- 
mine mass turbulent transport rate, concentration and velocity. Probability density 
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functions (p.d.f.s) of all the forementioned parameters were obtained from the data 
sets. Mean quantitites (Figs. 4, 6 and 7), second central moments (Fig. 5), correla- 
tion coefficients, skewness and kurtosis were calculated to characterize each data 
set. Following are the principal results from this study: 

1. Mixing for swirling flow was completed in one-third the length required for non- 
swirling flows (Figs. 8 and 9). 

2. The principal momentum turbulent transport was in the r-z plane, i.e. uv, and 
was attributed to the axial velocity gradients. Peak momentum turbulent trans- 
port rates were approximately the same as for the nonswirling flow condition 
(Fig. 6). 

3. The axial mass turbulent transport is gradient rather than countergradient as 
occurred for nonswirling flow. The peak axial mass transport rates were greater 
than the peak radial mass transport rates even though the axial concentration 
gradients were approximately one-seventh the radial gradients (Figs. 4 and 7). 

4. The mass turbulent transport process for swirling coaxial flow is very complica- 
ted. Mixing appears to occur in several steps of axial and radial mass trans- 
port coupled with a large radial mean convective flux. The transport process 
appears to begin with high concentration fluid from the inner jet mixing with 
fluid from the centerline recirculation zone. The diluted inner jet fluid is 
then convected by the negative axial velocities into the large eddy shear region 
between inner and annular streams. 

5. Axial mass turbulent transport correlation coefficients as high as 0.5 were 
measured. These correlation coefficients were less than the peak mass transport 
correlation coefficients obtained for nonswirling flow although the axial mass 
transport rate for swirling flow was greater than that for nonswirling flow. 
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VISUALIZATION OF FLOW 
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MEAN VELOCITY AND INNER JET 
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FIG. 5 

FLUCTUATING VELOCITY AND INNER JET CONCENTRATION PROFILES AT z = 25 mm 
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FIG. 8 
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THE INFLUENCE OF LARGE-SCALE MOTION ON TURBULENT 
TRANSPORT FOR CONFINED COAXIAL JETS* 

David C. Brondum and John C. Bennett 

Mechanical Engineering Department 
University of Connecticut 


The existance of large-scale coherent structures in turbulent 
shear flows has been well documented in the literature (ref.l). The 
importance of these structures in flow entrainment, momentum 
transport and mass transport in the shear layer has been suggested 
by several researchers (ref. 2, 3). Comparisons between existing 
models and experimental data for shear flow in confined coaxial jets 
(ref. 4, 5) reinforce the necessity of further investigation of the 
large scale structures. These comparisons show the greatest 
discrepency between prediction and actual results in the developing 
flow region where the large scales exist. It was also observed that 
the momentum transport rate comparisons were very bad. Finally, 
Schetz (ref. 6) has reviewed mixing flows and concluded that 
large-scale structures were essential aspects of future modeling 
efforts. 

The present analytical effort is a follow-on to a previous 
contract (ref. 7) under which the experimental data was collected. 

The data consists of all three velocity components as well as 
concentration data for a turbulent shear flow in a confined coaxial 
jet with a sudden expansion. The analysis undertaken is similar to 
that of Blackwelder and Kaplan (ref. 3) which was an analysis of a 
turbulent boundary layer. As they discovered, the bursts within the 
boundary layer were found to contribute greatly to the momentum 
transport. As part of the previous contract effort, flow 
visualization was used extensively and identified large-scale 
structures within the shear layer. The appearance of these 
large-scale structures in the mixing layer has led to the search for 
a burst detector along with detection criteria. Preliminary 
conditional sampling of the data has been performed and has revealed 
much information about the transport processes. All analysis to 
this point has been limited to data acquired at one axial location 
(102 mm from the origin). 

Various conditionally sampled parameters of the data have been 
plotted and trends observed. Initial efforts were directed toward 
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the use of these results in the identification of a suitable burst 
detector for the flow and also the development of better 
understanding of the effects on the flow (i.e., mass transport). 

One such plot is the scatter plot where the velocity fluctuation is 
plotted as a function of the concentration fluctuation or where 
fluctuation of one velocity component is plotted against the 
fluctuation of another velocity component. These types of plots 
help establish the existence of the large scale in the flow and 
assist in separating these structures from the general low level 
turbulence (Fig. 1,2). They are also beneficial in determining if 
the large scale correlation is significant (non-zero, Fig.l or zero 
Fig. 2) . 

Other interesting parameters can be established from 
conditonally sampled data pairs that are contained in the large 
structure. Based on the scatter plots, the initial choice for a 
burst detector was large excursions (positive or negative) in 
concentration fluctuations. The first parameter investigated was 
the large scale fraction (or ratio of samples with large 
concentration excursions compared to the total sample size). 

Plotting this ratio as a function of radius (at a given axial 
location) yields some trends in the data (Fig.l). The percentage 
appears to be maximized in the shear layer where the large-scales 
are known to exist. This adds credibility to the choice of detector 
of course. 

Having established some credibitity for the detector, the 
efforts were concentrated on the effects of large-scale structure on 
mass transport. Specifically, the ratio of conditionally-sampled 
mass transport to overall mass transport was plotted (Fig. 4). As 
with boundary layers, the large-scale structures in shear layers 
appear to exert considerable influence on mass transport. 

Having established procedures for investigating large-scale 
structures, future efforts will be directed toward several goals: 

(1) establishing a model of a typical 3-D structure, (2) 
investigating the regional extent of the influence of large-scale 
structures, (3) investigating the influence of large-scale 
structures on counter-gradient transport identified previously 
(Ref. 7), and (4) investigating the recently completed swirl data for 
similar large-scale effects. 
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VELOCITY VISUALIZATION IN GASEOUS FLOWS 


R. K. Hanson, B. Hiller, C. Hassa and R. A. Booman 
Department of Mechanical Engineering 
Stanford University 


Techniques yielding simultaneous, multiple-point measurements of velocity 
in reacting or nonreacting flowfields have the potential to significantly 
impact basic and applied studies of fluid mechanics. This research program 
is aimed at investigating several candidate schemes which could provide such 
measurement capability. The concepts under study have in common the use of 
a laser source (to illuminate a column, a grid, a plane or a volume in the 
flow) and the collection of light at right angles (from Mie scattering, 
fluorescence, phosphorescence or chemiluminescence) using a multi-element 
solid-state camera (100 x 100 array of photodiodes). 

The talk will include an overview and a status report of work in 
progress with particular emphasis on the method of Doppler-modulated 
absorption. This technique involves monitoring fluorescence from a plane 
in the flow which is illuminated with a sheet of light from a tunable, 
narrow-1 inewi dth laser source. The laser wavelength is set to coincide 
with an absorption transition of a seeded molecular species, usually 
iodine, which is Doppler-shifted according to the local velocity of the 
gas (relative to the direction of the laser beam). The resulting 
velocity-induced variations in absorption appear directly as variations 
in the fluorescence intensity, and hence quantitative, simultaneous 
measurements of fluorescence intensity at a large number of flowfield 
points can be used to infer the velocity at these points. Separate 
techniques and results are reported for supersonic and subsonic flows 
near room temperature. 
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OUTLINE OF PRESENTATION 

• MOTIVATION (NEED FOR MULTIPLE-POINT MEASUREMENTS) 

• SUMMARY: TECHNIQUES UNDER INVESTIGATION 

STREAK RECORDING 
MIE SCATTERING 

LASER-INDUCED PHOSPHORESCENCE 
•LASER MARKING 

•DOPPLER-MODULATED ABSORPTION/FLUORESCENCE ( DMA/F) 
SUPERSONIC FLOWS 
SUBSONIC FLOWS 

• EXPERIMENTAL RESULTS 

• RELATED RESULTS ON SPECIES/TEMPERATURE VISUALIZATION (AFOSR) 


STREAK RECORDING 

• SEEDED PARTICLES (DROPS) ARE ILLUMINATED BY CW LASER SHEET 

• MIE SCATTERING GIVES STREAKS ON FILM (RETICON) 

• GRADED INTENSITY BEAM (TRANSVERSE TO SHEET) GIVES POTENTIAL 
FOR 3-D RECORDING, IF PARTICLES ARE SPHERICAL 


TRANSFER OPTICS 
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LASER MARKING: CONCEPTS 

• A PULSED TUNABLE LASER IS USEC TO *NMK* A LINE OR A ARID 
IN THE FLO'VFIELD 

• LASER- INDUCED PHOSPHORESCENCE ' OR CHEMILUMINESCENCE } PROVIDES 
LONG-LIVED RADIATION FROM EXCITED FLUID ELEMENTS. OR 

• LASER- INDUCED PARTICULATE FORMATION ENABLES TRACKING BY 
MIE SCATTERING 

• MULTIPLE EXPOSURE OF GRID PA'Tr-’\s ARE RECORDED ON A SINGLE 
FRAME OF RETICON CAMERA (OR FILM) 

• CANDIDATE MATERIALS: 3IACETVL (DROPLET OR VAPOR) FOR PHOSPHORESCENCE 

NO, /CO FOR CHE” I LUMINESCENCE 
c 

SF 6 FOR PARTICLE FORMATION 



Flow Mark 1 ng by 
Laser-induced Phosphorescence 

Nitrojen ♦ 3.3 Ft acetal 
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DOPPLER-MODULATED ABSORPTION: SUPERSONIC METHOD 
Scans of Fluorescence Intensity vs Laser Frequency 
for an Off-Axis Pixel 



Relative Frequency CGHz] 

i 

I 

I 

1 

COMPARISON OF MEASURED AND CALCULATED VELOCITY FIELDS 
IN NEAR FIELD OF UNDEREXPANDED JET 
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EXTENSION TQ SUBSONIC FLOWS 


• FIX WAVELENGTH 


• MONITOR VARIATION IN INTENSITY BETWEEN FORWARD AND 
COUNTERPRCPAG flT ING BEAMS TO INFER VELOCITY 


• RECORD AT HIGH REPETITION RATE TO DETERMINE U(x,y.t) 


$(<•> 



u = C (iv/v ) 


Jet Axis 
I 



Counterpropaga ting 
Beams 


1'. = jS/(5S/?v). 

o 


= (AS/S q ) Cg(v o )/Og/0v). ] 
'o 


g = linesHape function 


VELOCITY VISUALIZATION 

SIMULTANEOUS MULTIPLE-POINT VELOCITY MEASUREMENTS BY 
SENSING DOPPLER-MODULATED LASER ABSORPTION WITH A DETECTOR SPRAY 



LO»«o 


Sketch «»t detection setup sho\% ini; the flow Held 
prohed hv a thin sheet ol laser light in the tenter plane »•! the 
jet and the camera imaging the fluorescence distribution 
perpendicular to the incoming sheet For < laritv, only one ol 
the lour sheets is drawn 



DIRECTION OF FLOW ; | 
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Flow Field of a subsonic round jet near nozzle exit. 
The important features are the potential core in which the 
center velocity remains constant and the growing annular 
turbulent shear layer. Three measured velocity profiles are 
shown. 


Velocity Visualization 
by Laser- I ndtieed F 1 uor escence 
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Characteristic radii of the subsonic round jet versus 
axial position. R., is the nozzle radius. R, is the radius of the 
hall -velocity point, and Rj is the radius of the shear-layer edge. 


SUMMARY 

OOPPLER-MOOULA TED ABS ORPTIO N/ FLUORESCENCE 

SIGNIFICANCE: MULTIPLE-POINT MEASUREMENTS IN UNSEEDEO FLOW 

STATUS: FIRST GENERATION EXPERIMENTS COMPLETED SUCCESSFULLY 

SUBSONIC FLOW (FIXED WAVELENGTH) 

SUPERSONIC FLOW (SCANNING CONCEPT) 

LIMITATIONS: LENGTHY RECORDING TIMES 

LIMITED VELOCITY RESOLUTION (1-5 M/SEC) 

FUTURE: IMPROVED TEMPORAL RESOLUTION 

STREAMLINED DATA PROCESSING 
COMBUSTION MEASUREMENTS 
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INVESTIGATIONS OF FLOWFIELDS FOUND IN TYPICAL 
COMBUSTOR GEOMETRIES 


David G. Lilley 
Oklahoma State University 


Studies are concerned with experimental and theoretical research on 2-D 
axi symmetric geometries under low speed, nonreacting, turbulent, swirling flow 
conditions. The flow enters the test section and proceeds into a larger chamber 

(the linear expansion ratio D/d = 2, 1.5 and 1) via a sudden or gradual expansion 
(side-wall angle a = 90 and 45 degrees). A weak or strong nozzle (of area ratio A/;. 
= 2 and 4) may be positioned downstream at x/D = 2 to form a contraction exit to the 
test section. Inlet swirl vanes are adjustable to a variety of vane angles with 

values of 4> = 0, 38, 45, 60 and 70 degrees being emphasized. The objective is to 
determine the effect of these parameters on isothermal flowfield patterns, time-riean 
velocities and turbulence quantities, and to establish an improved simulation in the 
form of a computer prediction code equipped with a suitable turbulence model. The 
goal of the on-going research is to perform experiments and complementary 
computations with the idea of doing the necessary type of research that will yield 
improved calculation capability. This involves performing experiments where time- 
mean turbulence quantities are measured and taking input conditions and running an 
existing prediction code for a variety of test cases so as to compare predictions 
against experiment. Hence the validity of turbulence model modifications can be 
assessed. If fact, they are also being deduced directly from the measured stresses 
and velocity gradients. 

Three Ph.D. Theses (Rhode, Jackson and Abujelala) and five M.S. Theses (Janjua, 
Yoon, McKillop, Sander, and Scharrer) have evolved in connection with the 
investigation over the last three years. New features of the present year's study 
include the investigation of: a more complex range of swirl strengths; swirler 

performance; exit nozzle sizes and locations; expansion ratios; and side-wall 
angles. Their individual and combined effects on the test section flowfield are 
being observed, measured and characterized. Recent studies on the test facility 
include: 

1. Swirler Performance - with five-hole pitot probe 1 and single-wire hot- 
wire 2 for time-mean and turbulence characteristics. 

2. Gross Flowfield Characteri zation - flow visualization has been achieved 
via photography of neutral ly-buoyant helium-filled soap bubbles and smoke 
produced by an injector and smoke wire. 3-5 

3. Time-Mean Flowfield Characterization - with five-hole pitot probe for 
time-mean velocity field and velocity gradients for a full range of swirl 
strengths. 3 -6 

4. Turbulence Measurements - with single wire 7 ’ 8 , cross-wire 9 and triple- 
wire 1 ^ hot-wires, including an assessment of the accuracy and directional 
sensitivity of the single-wire technique 11 , enabling all Reynolds stress 
components to be deduced. 
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5. Turbulence Modeling - limitations and empirical extensions of the k-e 
model have been analyzed 12 , but relating Reynolds stresses to appropriate 
time-mean velocity gradients are yielding effective model improvements 13 . 

6. Computer Predictions - using an advanced computer code 14 , tentative 
predictions^* have now been supplemented by predictions made using 
recently-acquired realistic inlet conditions for a complete range of swirl 
strengths 15 , including expansion ratios and downstream nozzle effects 16 . 
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INVESTIGATIONS OF FLOWF I ELDS FOUND IN TYPICAL 
COMBUSTOR GEOMETRIES 
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1 . Introduction 
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Research 

Object i ves 

1 . 

To 

determine the effect of: 


* 

1 

swi rl vane an,gl e 4> 


* 

si de-wa 11 angl e a 


* 

downstream blockage area ratio 



and location L/D 


on : 



★ 

isothermal flowfield patterns 


* 

time-mean velocities 


★ 

turbulence quantities 

2. 

To 

establish an improved simulation 


the 

form of: 


★ 

computer prediction code 


★ 

suitable turbulence model 


1 . 3 Research Approach 

1 . T ime-mean fl owf i el d characterization by five- 
hole pitot probe measurements and by flow 
visual ization. 

2. Turbulence measurements by a variety of single- 
and multi-wire hot-wire probe techniques. 

3. Flowfield computations using the computer code 
developed during the previous year's research 
program . 
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2 . Progress 

2.1 Progress During First Year (1981) 

1 . Test facility - design and construction, including 

variable angle swirler and expansion blocks. 

2. Experimental techniques - flow visualization, five-hole 

pitot probe and reduction computer code. 

3. Computational code development - advanced version of a 

pr i mi t i ve - va r i abl e solution procedure - STARPIC. 

4. Flowfield characterization - emphasis on time-mean flow 

character, using flow visualization, five-hole 
pitot probe and tentative predictions. 


2.2 Progress During Second Year (1982) 

1. Higher swirl strengths - flow visualization and five- 
hole pitot probe measurements. 

2 . Swirler ef f ectiveness - exit profiles of u, v, w, and p 
accurately established. 

3. Flowfield predictions - the need for specifying the 

inlet conditions very precisely. 

4. Downstream blockage - effects of area ratio and axial 

location studied by flow visualization, five-hole pitot 
probe and corresponding computer predictions. 

5 . Turbulence measurements - in nonswirling and swirling 

flow via one-wire, two-wire and three-wire hot-wire 

methods . 
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2.3 Progress During Third Year (1983) 

1. Sma 1 ler expansion ratios - D/d = 1, 1.5 and 2 

2. Turbulence measurements - at higher swirl strengths, 
with and without downstream blockage 

3* Turbul ence model i ng - developments from measured 
Reynolds stresses and time-mean velocity gradients 
4 * Flowfield predictions - for corresponding situations 


3. Results C4> = 45 degrees] 


3.1 Swirler Performance 






3 . 2 Time-Mean Velocity Profiles 
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3.4 Turbulent Shear Stress Profiles 





3 . 5 P redicted Prof i 1 es 

Using optimum k-c parameters of Ref. 12 
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Closure 

1. Measurements of effect of swirl, expansion angle and 
contraction blockage _or^ flowfield patterns. 

2. Computer prediction studies of associated phenomena. 
Turbulence modeling developments for the simulation of 
swirling recirculating flow. 
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DILUTION JET EXPERIMENTS IN 
COMPACT COMBUSTOR CONFIGURATIONS 

Isaac Greber 
James Zizelman 

Department of Mechanical and Aerospace Engineering 
Case Western Reserve University 

This project concerns the effects of cooling jets on the 
velocity and temperature fields in a compact reverse flow 
combustor. The work is motivated by the need to limit the 
temperatures of post-combustion gases in jet engines to values 
within the endurance capabilities of turbine blades. The 
application requires not only that the temperature be kept 

sufficiently low but also that a suitably tailored temperature 

profile be provided at the combustor exit, with higher 

temperatures generally permissible at the blade tip than at the 
blade root because of higher centrifugal loads at the root. 

Flows in reverse flow combustor accelerate both 
longitudinally because of area changes and transversely because 
of flow turning. The current project started with flow 
visualization experiments in water, using aqueous solutions of 
zinc bromide to model the relatively higher density of cooling 
jets. These flow visualization experiments were conducted in 
simple two-dimensional configurations designed to examine 
separately the effects of longitudinal and transverse 

acceleration. The next phase of the work consisted of 
temperature measurements in a model reverse flow combustor, using 
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a rake of thermocouples which could be moved both longitudinally 
and transversely to provide transverse temperature profiles at a 
number of longitudinal stations. Single jets, and rows of 
cooling jets of different spacing, were injected at several 
locations with varying flow rates and temperatures to produce a 
useful range of momentum ratios and density ratios. In addition, 
a semi-empirical calculational model was developed to predict the 
behavior of a single cooling jet in the reverse flow combustor 
conf iguration . 

Results of these experiments show that single jet 
temperature trajectories are swept toward the inner wall of the 
turn, whether injection is from the inner or outer wall. A 
widely spaced row of jets produces a trajectory similar to that 
of a single jet. As spacing is reduced, jet penetration is also 
reduced, and the cooling jets tend to remain close to the wall 
from which they are injected. These results suggest that 
suitable cooling and temperature distribution tailoring can be 
accomplished without injecting cooling jets upstream of the turn, 
and thus it appears that combustors can be made significantly 
smaller than current designs. 

The current phase of this project is directed toward the 
acquisition of velocity measurements in the turn section of the 
combustor. To this end a combined pitot-static tube and 
thermocouple probe was constructed. Following successful tests 
of this probe, a five-probe rake using the same design was 
constructed. Velocity and temperature distributions are being 
obtained using this rake. Because of the large quantity of data 
required for adequate representation of the velocity and 
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temperature fields, the data is being directly computer acquired. 

An attempt to develop a semi-empirical calculational model 
for a row of jets, similar to that for a single jet, was limited 
by the absence of fundamental mixing information for rows of 
jets. Consequently, experiments were performed to obtain 
entrainment rates for rows of jets. In these experiments volume 
flow rates were deduced by integrating velocity distributions 
obtained from hot wire anemometer transverses. A Master's thesis 
describing these results is in preparation. In a companion 
effort, vortex models of a jet in cross-flows have been developed 
in an attempt to understand the relative importance of viscous 
and inviscid phenomena in determining the jet trajectories. 

Throughout this work we have been fortunate to have the 
support and advice of Steve Riddlebaugh of the NASA Lewis 
Research Center, and we are pleased here to thank him. 
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FIGURE 3. CUTAWAY SKETCH OF TEST COMBUSTOR 



FLOW 


EXHAUST 



MEASURING 

STATIONS 


6 ROTATING- 


+ v i / -Tl-I TRAVERSING 

1 5 THERMOCOUPLE 
RAKE 


S 1 

INJECTION 
POINT B 


INJECTION . 

POINT C 
START OF 
TURN SECTION 


FIGURE <i. COMBUSTOR TURN SECTION GEOMETRY 
AND TEMPERATURE MEASURING STATIONS 
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PROFILES 

FIGURE 5. TEMPERATURE DISTRIBUTION IN THE COMBUSTOR 
TURN WITHOUT DILUTION JETS 



FIGURE 6 . JET TRAJECTORIES SHOWING EFFECT 
OF MOMENTUM RATIO ON A JET 
INJECTED FROM THE INNER WALL 
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DENSITY RATIO 2.75 
MOMENTUM RATIO 
5.67 
n 7.65 
o 9.01 


FIGURE 7. JET TRAJECTORIES SHOWING EFFECT 
OF MOMENTUM RATIO ON JET 
INJECTED FROM THE OUTER WALL 
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FIGURE 8. COMPARISON WITH MODEL FOR A SINGLE JET 
INJECTED FROM THE INNER WALL 
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FIGURE 9. COMPARISON WITH MODEL FOR A SINGLE JET 
INJECTED FROM THE OUTER WALL 
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FIGURE to. BEHAVIOR OF A WIDELY SPACED ROW OF JETS 
(SPACING RATIO OF 7.41) INJECTED FROM 
THE INSIDE WALL 
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FIGURE 11. BEHAVIOR OF A CLOSELY SPACED ROW OF JETS 
(SPACING RATIO 2.47) INJECTED FROM 
THE INSIDE WALL 
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FREE STREAM TURBULENCE AND DENSITY RATIO EFFECTS ON THE 
INTERACTION REGION OF A JET IN A CROSS FLOW 1 


9 

C.E. Wark“ and J.F. Foss 
Department of Mechanical Engineering 
Michigan State University 


Jets of low temperature air are introduced into the aft sections of gas 
turbine combustors for the purpose of cooling the high temperature gases and 
quenching the combustion reactions. Research studies, motivated by this 
complex flow field, have been executed by introducing a heated jet into th| 
cross stream of a wind tunnel. The investigation by Kamotani and Greber 
stands as a prime example of such investigations and it serves as the prin- 
cipal reference for the present study. 

The low disturbance level of the cross stream, in their study and in similar 
research investigations, is compatible with an interest in identifying the 
basic features, of this flow field. The influence of the prototypes' strongly 
disturbed cross flow is not, however, made apparent in these prior 
investigations. The present study provides a direct comparison of the thermal 
field properties for a low (u/U^ = 0.67%) and a high (u/V^ = 34%) 

disturbance level condition. 

A novel technique was used for the data acquisition (Figure 1). Sixty-four 
fast response (\ - 4msec) thermocouples were simultaneously sampled (and 
corrected for the time constant effect) at a downstream plane close to the jet 
exit (x/d*3-4). Various measures were used to characterize the thermal fields 
for the disturbed and undisturbed conditions, for two different momentum flux 
2 2 

ratios ( J * p . V 7 / P V~=16,64), and for three overheat conditions 

j j 00 00 

(1*22.2, 47.1 and 61.1). 

Two forms of data acquisition were used for this study. Stochastic values 

were obtained from triplet, values: T (t-6t), T ( t ) and T (t + 6t), where 

K k k 

K=l....fa4 is the thermocouple designation and the tit ( 6 1=0 . 64msec ) values 
were used to form the central difference time derivative from which the 
corrected temperature value: { T ( t > } , was determined. These values were 

K C 

retained for t, t +At , t + 2 At . . . . t +NAt where N* 5,000 for the undisturbed and 
10,000 for the disturbed condition and At=100ms which ensured the statistical 
independence from its neighboring value. Instantaneous samples: At = 0.64ms 
and N«1 ,225 were stored for J*16, 64 and T. * 61.1 °c. 


0 Supported by NASA Grant 
M.S. 1984, presently at 


NAG 3-245 

Illinois Institute of Technology 
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RESULTS 

Histograms, formed from the independent samples, were sufficiently smooth to 
approximate a probability density function; examples are shown in Figure 2 for 
a thermocouple from the central region of the jet. A striking result, from 
all such histograms, is that the peak (non-dimensional) temperature did not 
exceed 0.25. Hence, even at the relatively close x/d locations of the present 
study, molecular diffusivity has played a dramatic role in the reduction of 
the temperature for the fluid elements of the jet. 

The individual histograms, and the isothermal patterns that are fit to the 
mean values of the corrected temperatures, both reveal that the magnitude of 
the overheat exerts a significant dynamic effect on the jet in a cross flow 
problem; see Figures 3 and 4. It is pertinent to note that the increasing 
distance of the thermal center of gravity, from the source plate with in- 
creasing jet temperature, is opposite to that which would be expected from 
buoyancy and enhanced entrainment effects. The former agrees with the 
Kamotani and Greber results, the latter is in disagreement. The physical 
reasons for this unexpected behavior are not apparent; the consistency of the 
trend: an increased penetration of the thermal field with an increased jet 
temperature, is consistently observed and serves as a primary focus for 
further study. A further comparison with the results from Kamotani and Greber 
(see Figure 5) reveals a substantial difference between the jet penetration 
for the two studies. This difference is (at least in part) attributed to the 
present use of a sharp-edge nozzle and to the large (6/d - 1.3) boundary layer 
for the present study. 

One clearly apparent effect of the large disturbance condition is the mi- 
gration of the centroid of the thermal field; see figure 6. Similar data for 
all of the instantaneous scans (6 disturbed, 3 undisturbed) are summarized in 

Table I; note: K is the normalized correlation coefficient between the 
s t 

vertical (s) and transverse (t) displacements of the thermal center of 
gravity. 


TABLE I. Cross stream disturbance effect 
on centroid of temperature field 


J 


und isturbed 

c.f . 

16 

undisturbed 

c . f . 

64 

disturbed c. 

, f. 

16 

disturbed c. 

. f. 

64 


s 

c.g. 

t 

c-g. 

St 

.341 

.253 

.067 

.233 

.233 

-0.13 

1.84 

1.45 

-0.33 

1 .16 

1 .12 

-0.45 


The lower mean values, of the disturbed cross stream mean isotherms, are 
apparently caused by these large centerline migrations; figure 7 shows an 
instantaneous isotherm oattern that reveals maximum temperature values of the 
same order as those for .he undisturbed case. 
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Figure la. Schematic of Jet Supply and Exit Conditions 

Note: Lateral (i.e.,z) position of the jet corresponds to the 

undisturbed position. The jet exit was moved into the 
low speed portion of the shear layer: u/U =0.25, for the 
disturbed flow condition. 0 


u. 



Figure lb. Schematic Representation of Jet Trajectory, Showing Lab 
coordinates (x, y, z) and Jet Coordinates (t, n, ;) 
Thermocouples were located in the y-z plane at x=3.5(fcr J=641 
and 3.85 (J*16). 
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Figure lc. Schematic Representation of the Temperature 
measurement system 
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Figure le. Thermocouple Array Assembly 


ORIGINAL PAGE !9 
OF POOR QUALITY 





















ORIGINAL PAGE IS 
OF POOR QUALITY 


.J HOI MMM VAluat 



Figure 4. Overheat Influence on s Component Centroid 

(calculated using laothera contours) 



Figure 5d. Location of Temperature Centerline (R"4) 
from Kamotanl and Greher (1971) 
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Figure 5b. Legation ot Temperature Centerline (R»6> 
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Figure 6. Variation of s Component of Centroid with Time 

J ‘ 16 ’ T jet ' T - = 61 ' 10C 
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MODELING OF DILUTION JET FLOWF I ELDS 

by J.D. Hoi deman 
NASA Lewis Research Center 


and 

R. Srinivasan 

The Garrett Turbine Engine Company 


Consi der at 1 ons of dilution rone mixing in gas turbine 
combustion chambers have motivated several studies of the 
mixing char acter i st l cs of a row of jets injected normally into 
a flow of a different temperature in a constant area 
duct."«- l ~ 4 Recently, experiments have been performed 
to extend these investigations to include geometric and flow 
variations char acter i st i c of most gas turbine combustion 
chambers, namely a variable temperature mainstream, flow area 
convergence, and opposed rows of jets, either in-line or 
staggered see figures 1 It 2. 

The present paper will compare temperature field measurements 
from selected cases in these investigations with distributions 
calculated with an empirical model based on assumed vertical 
profile similarity and superposition'' and with a 3-D 
elliptic code using a standard K-E turbulence model. 7 
The results will show the capability (or lack thereof) of the 
models to predict the effects of the principle flow and 
geometric var i abl es. 


Variations with Orifice Size and Spacing. At constant 

orifice area, changes in orifice size and spacing can result in 

jets which vary from under -penetrat i on to over-penetr at i on . 

This is shown in figure 3a) and b) , for jets from closely 
spaced small orifices and widely spaced larger orifices 
respectively. The empirical model reproduces the data very well 
in the small orifice case, since the data are consistent with 
the major assumption in the empirical model, that all vertical 
temperature distributions can be reduced to similar Gaussian 
profiles. The empirical model does not do as well in the larger 
orifice case however, as the jets have impinged on the opposite 
wall and the vertical profiles are not similar. 


The analytical model calculations made with approx i matel y 
20,000 nodes, although in qualitative agreement with the data, 
under — predict the mixing. That is, the temperature gradients, 
especially in the transverse direction, are too steep. This 
result is typical of the analytical model calculations to be 
shown in this paper . For the smal 1 -or i f l ce case a coarse-grid 
calculation using less than 6000 nodes was also performed. This 
solution illustrates the diffusive nature of the calculation 
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and the significant influence of qrid selection on the solution 
obtai ned . 

Coupled Spacing and Momentum Flux Ratio . Examination of 
the experimental data revealed that similar profiles can be 
obtained over a range of momentum flux ratios, independent of 
orifice diameter, if orifice spacing and momentum flux ratio 
are correctly coupl ed. v This is shown in figure 
4a) to c ) . 

In all of the combinations shown here, the empirical model 
results are in very good agreement with the data, as the data 
are consistent with the Gaussian profile assumption. The 
analytical model calculations using appr ox 1 mal el y 20,000 nodes 
for these cases aqree qua 1 l t l t l vel y with the data, as in the 
previous figure. In the medium momentum flux ratio case, a 
second calculation was performed with the same total number of 
qrid points, but with the nodes slightly more concentrated in 
the vicinity of the jet exit. As can be seen in figure 4b), 
these two results are not substant i al 1 y different. 

Variable Temperatur e Mainstream. The influence of a 
non— l sother mal mainstream flow on the profiles for medium 
momentum flux ratios with S/Ho^.G and Ho/D=4 can be seen by 
comparing figures 5 & 6. The shape of the experimental profiles 
in figure 6 suggests modeling them as a superposition of the 
upstream profile and the corresponding j et s—i n-an—i sother mal 
mainstream distribution. v This gives only a crude 
approximation however, as seen in the empirical model results, 
because of the cross-stream transport of mainstream fluid due 
to the blockage, which is not accounted for in superimposing 
the distributions. 

In the variable temper ature mainstream case the analytical 
model results aqree well with the experimental data, especially 
on the jet centerplane, but the transverse mixing is 
underpredicted, as in the corresponding isothermal mainstream 
case in figure 5. 

Opposed Rows of In-line Jets. For opposed rows of jets, 

with the orifice centerlines m-line. the optimum ratio of 

orifice spacing to duct height is one-half of the optimum value 

for single— side injection at the same momentum flux 

ratio.'" As as example consider the single-side case 

with S/Ho=.5 and Ho/D=4 in figure 5 and the opposed row of 

in-line jets with S/Ho=.25 and ho/D— 8 in figure 7. 

Iho empirical model predicts the opposed-jet case very well, 
verifying the primary assumption that the effect of a plane of 
symmetry is similar to that of an opposite wall. 1 ’-"" 

Mote that the experimental profiles on both sides of the plane 
of symmetry support the Gaussian profile assumption. The 
analytical mode) results show the steep transverse and lateral 
gradients seen in almost all of the previous calculations also, 


but are in otherwise good agreement with the data. 

Opposed Rows of Staggered J ets. For opposed rows of jets, 

with the orifice centerlines staggered, the optimum ratio of 
orifice spacing to duct height is double the optimum value for 
single-side injection at the same momentum flux 
ratio."” As an example consider the single-side case 
with S/Ho=.5 in figure 5, and the opposed row of staggered- jets 
with S/Ho=l in figure 8. 

The empirical model does not handle this complex case well, as 
the fluid dynamic interactions here are not amenable to a 
direct extension of the simple Gaussian profile and 
superposition type modeling appropriate <or most of the 
single-side and opposed-jet cases of interest. The analytical 
model calculations give slightly better agreement with the data 
than does the empirical model, and would be expected to improve 
with overall improvements in the capability of the 3-D codes. 
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TURBULENCE CHARACTERISTICS OF AN 

AXI SYMMETRIC REACTING FLOW 

Richard D. Gould, Warren H. Stevenson 
and H. Doyle Thompson 

Purdue University School of 
Mechanical Engineering 


Turbulent sudden expansion flows are of significant theoretical and practical 
importance. Such flows have been the subject of extensive analytical and experi- 
mental study for decades, but many issues are still unresolved. Detailed informa- 
tion on reacting sudden expansion flows is very limited, since suitable measurement 
techniques have only been available in recent years. The present study of reacting 
flow in an axisymmetric sudden expansion was initiated under NASA support in 
December 1983. It is an extension of a reacting flow program which has been carried 
out with Air Force support under contract F33615-81-K-2003. Since the present 
effort has just begun, results are not yet available. Therefore a brief overview 
of results from the Air Force program will be presented to indicate the basis for 
the work to be carried out. Details may be found in reference 1. 

Laser velocimetei measurements of mean streamwise velocity and turbulence in- 
tensity were made in the highly turbulent flow field following a sudden pipe ex- 
pansion. Both isothermal and reacting flows were studied. A fused quartz test 
section was used to permit laser velocimeter measurements throughout the flow field 
for x/H values from 0.33 to 15. A lean partially premixed propane-air mixture 
(<p = 0.28) was used to keep the wall temperature low enough so that steady state 
operation was possible. The inlet flow condition was that of fully developed 
turbulent pipe flow with a centerline velocity of 22 m/s. The corresponding 
Reynold's number based on step height H and inlet centerline velocity U. was 
5.5 x io4. 1 

Three complete sets of LDV measurements were made. Two were in isothermal 
flow (biased and unbiased) and one was made in the reacting flow (biased) . The 
biased and unbiased cold flow measurements were compared to determine the effects 
of velocity bias on the measurements. The hot and cold flow measurements were 
compared to determine the effect of combustion on the structure of the flow field. 
(Unbiased hot flow data could not be obtained because of particle seeder limita- 
tions. ) 

The scope of the investigation is outlined in figure 1. The experimental 
apparatus is shown in figures 2-4 and test parameters are presented in figure 5. 
Examples of some of the data obtained are given in figures 6-15 and conclusions 
which can be drawn are given in figure 16. A summary of the further research 
planned under the present NASA sponsored program is presented in figure 17. 
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SCOPE OF THE STUDY « 

l* Measured streamwise neon velocities and turbulence 
Intensities. 

2. 3 different mecsurements were mode at eacn measurement 
location to form 3 complete data sets. Tney are 
identified as: 

♦ biased 1 

♦ unbiased J 

♦ biased 

3. Reattacnment point was located. 

4. Integrated mass flux at each measurement Diane was 
calculated ana was used as a continuity cnecx. 


5. Measured temperature at a Diane located 17 steD heignts 
aownstream of the sudden expansion, 

fa. Comparison of cold flow aata with 2/E/FiX, 

*AF Contract ^lS-Sl-K-TOOS 


Figure 1 


cold flow 
hot flow 
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Telescope 


Polor izotion 


^ 0 ~§ 


• 2*68 mm 


1 1 * 44 mm 



(movable) (fixed) 


fiijure i HIV Optics Package 
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All Meas temenls in mm 
Not to Scale 




Injure 4 Geometry of Axi&yiHmetr ic lest Section 


TEST PARAMETERS 


FLTW rAFA v E'E-.S 

U t * 22.07 m/s ( FITOT TUBE MEASUREMENT) 
Se 3 * 1.1 « 10 5 
Se M * 5.5 * 10 4 

OVERALL F/A * 0.3183 (GASEOUS PROPANE) 
OVERALL EQUIVALENCE RATIO, > * 0.23 


QLQHETRY 


0j N * 76.2 in ( 3 in) 

0 0uT * '52.4 nm (6 in) 

AREA RA'IO, A r . A 

STEP HEIOrlT, ri = 38.1 nvn (1.5 in) 


PARTICLE ARRIVAL RATE: 
SAMPLE RATE: 

SEED: 


UNBIASED COLO 
> 20,000 

50/sec 

OOP 


BIASED COLD 
500-1500/ sec 

FREE" (4700/sec) 


BIASED HOT 
500-1500/sec 

FREE" (4700/sec) 


Figure 5 
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Figure 8. Measured Mean Streaniwise Velocity Profiles in Cold Flow 
(o - unbiased cold flow. A - biased cold flow, --- 2/E/FIX prediction) 


NOkflALIZLU* 1 NTLGRATLL) MASS FLUX 
(COLD F10W) 


*/« UNHIASIU I1ASIO 

1/3 1.00 

1 1.044 1.065 

3 1.056 1.200 

5 1 . 002 1.129 

1 1.012 1.206 

9 1.032 1.238 

•I 1.019 1.212 

15 1031 1.166 


•VALUES NORMAL 1 2 E WITH ASSUMED 1/7 POWER LAW 
INLET VELOCITY PROFILE 

Figure 9 
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ORIGINAL P2GE tS 
OF POOR QUALITY 



u/U, 

Figure 10. Measured Mean Streaiuwise Velocity Profiles in Cold Mow and in Reacting Flow 
(o - unbiased cold flow, U • biased hot flow) 
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CONCLUSIONS 


FUTURE WORK* 


1. Velocity measurements were made in a sudden exDansion 
turbulent reacting flow at atmosDheric pressure. Few 
difficulties were encountered in making these measure- 
ments. 

2. The reacting flow exhibited several structural differences 
relative to the cola flow: 

a) The recirculation zone was smaller ana more intense 
in the reacting flow. Reattachment occurred at X/H = 
7.4 for the hot flow as compered to 8.6 for the cold 
flow. 

b) The position of peak turbulence intensity was ai <f erent 
in the reacting flow, although the maximum normalized 
turbulence intensity was approximately the same in 
both cases. 

c) Local turbulence intensity was reduced by up to 60 
percent in the reacting flow due to the higher local 
velocities resulting from the heat release. 

3. velocity bias had the expected effect on tne mean velocity 
data. A significant shift in the normalized turbulence 
distribution in the shear layer was found to exist due to 
bias. 


Figure 16 


1. Develop 2 color - 2 component LDV system 

• Measure mean and fluctuating velocity 
simultaneously for two components 

• Measure Reynolds stress correlation 

• Obtain power spectrum and length scales. 

2. utilize correction lens system to make radial velocity 
measurements 

• Allows uv R Reynolds stress correlation 
measurement 

• Important if studying e'fects of swirl 

3. Utilize high freauency response thermocouple to make 
dynamic temperature measurements in reacting flow 

• Balance terms in the turbulent kinetic energy 
conservation eauation in both hot and cold flow 

• Identify regions where the various turbulence 
mechanisms dominate in the hot and cold flow 

Production 

Advection 

Diffusion 

Dissipation 

•NASA Grant NAG 3-502 


Figure 17 
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DIRECT NUMERICAL SIMULATION OF REACTING FLOWS 

James J. Riley and Ralph W. Metcalfe 
Flow Research Company 


The objectives of thi^ work are (i) to to extend the technique of direct numeri- 
cal simulations to turbulent, chemically reacting flows, (ii) to test the validity of 

the method by comparing computational results with laboratory data, and (iii) to use 
the simulations to gain a better understanding of the effects of turbulence on chemi- 
cal reactions. In particular, we address the effects of both the large-scale struc- 

ture and the smaller-scale turbulence on the overall reaction rates, examine the rela- 
tionship between infinite reaction rate and finite reaction rate chemistry, and com- 
pare some of the results of our calculations with existing theories and laboratory 
data. The direct numerical simulation method involves the numerical solution of the 
detailed evolution of the complex turbulent velocity and concentration fields. Using 
very efficient numerical methods (e.g., pseudospectral methods), the fully nonlinear 
(possibly low-pass filtered) equations of motion are solved and no closure assumptions 
or turbulence models are used. Statistical data are obtained by performing spatial, 
temporal, and/or ensemble averages over the computed flow fields. 

The scope of work to do this involved the following. First, existing computer 
codes were modified to treat the present problem. Next, extensive numerical testing 
of the computer codes was performed. Finally, in order to examine the effects of the 
mixing layer turbulence, both the large-scale structure and the smaller-scale turbu- 
lence, on the overall reaction rates, a sequence of three problems was computed: (i) 

reactions on a unidirectional (one-dimensional) mixing layer, (ii) reactions on a 
mixing layer experiencing large-scale, two-dimensional vortex rollup, and (iii) reac- 
tions on a tti ree-d imensional turbulent mixing layer. The simulations of the two- 
dimensional mixing layer with vortex rollup are intended to model the large-scale 
structure in the mixing layer, whereas the three-dimensional simulations contain both 
the large-scale structure and the smaller-scale turbulence. 

The numerical testing involved the comparison of computed results with exact solu- 
tions for a number of different cases. Both rigid body rotation and vertex rollup 
flow fields were used. The work greatly extended the results of Orszag (*ef. 1) for 
the advection of a passive scalar on a rigidly rotating flow field (the color problem) 
to include also diffusion, chemical reaction, and more complex flows. We have found 
that high accuracy of the spectral methods observed in the advection case is also 
obtained when these further complications are present. Our results indicate that 
spectral numerical methods may prove to be useful in the future both for solving the 
model equations for combustor processes as well as for future studies of chemically 
reacting turbulent flows employing direct numerical simulations. 

The approach of direct numerical simulations allowed extensive examination and 
interpretation of the reaction process. From the one-dimensional simulations, we 
found that we could easily compute finite reaction rates near the fast reaction limit, 
that the results were fairly insensitive to the initial conditions (for the class of 
initial conditions computed), and that the results were in reasonable agreement with 
theoretical predictions. 
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For the two-dimensional simulations, we found that in all of the cases computed, 
the vorticity field and the product field approximately coincided. From the computa- 
tion of volume averages, we observed the enhancement of the overall reaction rate due 
to the vortex rollup. It also appeared that the merging of vortex cores was a more 
significant mechanism in increasing the overall reaction rate than the straining of 
the reaction interface. Significant species segregation was apparent, so that, for 
example, the product of the average concentrations was approximately equal to and 
opposite the correlation between the fluctuating concentrations. Some of the higher 
order correlations were also examined. 

In the three-dimensional simulations, the contour plots indicated that the vortex 
rollup takes longer to develop, and the vortices and braids are not as distinct as in 
the two-dimensional case. Also, the vortices that develop are not strongly corre- 
lated laterally. Both the contour plots and the statistical results indicated that 
the spatial segregation was also not as strong as in the two-dimensional case, 
probably due to the weaker vortex rollup as well as the effects of smaller-scaie , 
three-dimensional turbulence. 

Comparisons were made between the simulation results and results using similarity 
theory. Approximately linear growth rates of various computed length scales, in- 
cluding the mean velocity half-width, the mean vorticity thickness, and the mean pro- 
duct thickness, were obtained and were in agreement with the theory. Similarity 
scaling was found to collapse quite well the results for the average reactant concen- 
trations, the rms fluctuating reactant concentrations, the concentration correlations, 
the average product concentrations, and the rms fluctuating product concentrations. 

Some limited compai isons were made with laboratory data. Computed profiles that 
were qualitatively similar to corresponding laboratory profiles were obtained for the 
average reactant concentrations, the rms fluctuating reactant concentrations, the 
average product concentrations, the rms fluctuating product concentrations, and the 
concentration correlations. 

We have made some comparisons with existing theories. Donaldson and Hilst have 
suggested, in addition to using the equations for the average concentrations, in- 
cluding the equations for the concentration fluctuations and correlation (ref. 2). 
These equations can be closed by neglecting certain triple moments when compared to 
certain lower order terms. However, our results indicate that the triple moment 
terms are as important as other terms in the equations, so that an assumption of this 
type will probably lead to poor predictions. Mason and Spaulding have proposed a 
model for the mean reaction term, suggesting that it will be proportional to the 
average concentration of the lean species divided by a turbulent time scale (ref. 3). 
We have found that such an assumption will only be moderately successful if applied 
to our case. Finally, Toor has suggested estimating the concentration correlation of 
the reacting species in terms of that for the nonreacting case, which is much easier 
to model (ref. 4). Although this was proposed mainly for statistically homogeneous 
flows, we find that it is a reasonable approximation for our reacting flow 
s imulations. 
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OBJECTIVES 


• Begin the Extension of Direct Numerical Simulations to Combustor 
Flows 

• Develop computational methodology 

• Use Direct Numerical Simulations in Order to 

• Better understand the effects of turbulence on chemical reactions 

• Suggest model improvements 

• Test Validity of the Method by Comparing Comnutational Results 
with Laboratory Data and Theoretical Models 


PROBLEM CONSIDERED 

• Temporally Growing, Three-Dimensional Mixing Layer 

• Binary, Single-Step, Irreversible Chemical Reaction 

(A + nB — ► products) with negligible heat release 

• Very Small Mach Number 

• Initial Conditions - NonPremixed Species 
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PROBLEM CONSIDERED 



«'(z) 



C R (x.O) 




SCOPE OF WORK 


ORIGINAL PAGE IS 

OF POOR QUALITY 


• Modify Existing Computer Codes to T reat the Present Problem 

• Perform Numerical Testing of the Computer Codes 

• Compute the Following Sequence of Problems 

• Unidirectional (one-dimensional) flow 

• Large-scale, two-dimensional vortex rollup 

• Three-dimensional turbulent mixing layer 


EQUATIONS OF MOTION 

• Conservation of Momentum 

00 i a a 2 

o U. + U.' 5 U. - - — « p + V « 2 U - 

O t I I O X: I Q 0 X. Ox. I 

'I ** I J 

• Conservation of Mass 


-5 — u. = 0 

d X: I 


• Conservation of Species A and B 

a 


o_ 

a 

a 


0 C + u. j* C — — R C C + D 

0 1 1 j 0Xj 1 112 


dj C 

1 dxj 1 


C_ + u. o C - R C C + D « 2 c 


0t 2 T u j 0 Xj 2 


2 1 2 2 0x^ 2 
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• Initial Boundary Conditions 


DIRECT NUMERICAL SIMULATIONS 


• Compute the Solution of the Navier-Stokes Equations Numerically 

• Initial Conditions (or Upstream Conditions) 

Random, turbulent-like - so that the flow field starts out 
as turbulent-like - or 

Small amplitude perturbations -- compute the flow through 
transition to turbulence 

• Compute the Time Development of the Detailed Structure of the 
Flow Field 

• No Reynolds Averaging 

• Statistical Results are Obtained by Performing Spatial, Temporal, 
or Ensemble Averages over the Computed Flow Field 

• Analogous to Laboratory Experiments 


NUMERICAL METHODS 


Pseudo-Spectral Methods 

• Expand the dependent variables in fast-converging series 


E.G., 


ulxl - N<K 2lkle ' kX 


• Resulting Equations 


+ i£ u.(k-k')k: C 1 (k') = - R 1 2 C^k-k'jC^k') -Djk^lk 


Convection 


Reaction 


(k> 
Diffusion 
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NUMERICAL TESTS 


• Color Problem (Rigid Body Rotation of Passive Scalar) 

• Convection Alone (Orszag, Seinfeld et al.) 

• Convection Plus Reaction (Exact Solutions Available) 

• Convection Plus Diffusion (Exact Solutions Available) 

• One-Dimensional, Diffusion-Reaction Problem 

• Test Against Exact Solutions 

• Two-Dimensional Vortex Rollup 

• Convection, Diffusion, and Reaction 

• "Exact" Solution for 6 = C - C is Known 

from Previous Calculation ^ 2 


NUMERICAL TEST PROBLEMS 


• Two Dimensional Rigid Rotation (Color Problem) 

( 2 _ 2 

*0 2 0 

• Analytical Solution for Identical Initial Concentration Fields 
(No Diffusion) 


C(x',z',t) = Qtx' i Z'.C) — 

1 + RtC(x',z',0) 

• Analytical Solution with Diffusion (No Reaction) 


C(x',z',t) 


1 

nvn 

x' 2 /«2 

^4 

1 / 1 + 4Dt \i 

CAP 

'l +4B1\ 

/l + 4Dt U 

, + 4Dt \ 

A 4 )< 

4) 

\ x o A 

z o / 
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Initial Concentration Field for Color Problem Teats 



Computed Concentration Field After One Revolution - No Diffueion 
or Reaction 
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Computed Concentration Field Alter One-Halt Revolution 
Di t f us ion , K " 1 .0 


Computed Concentration Field After One-Half Revolutii 
Di t fusion , R • U .0 
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TWO-DIMENSIONAL SIMULATIONS 


■ 


• Case I: Most Unstable Mode Based Upon Linear Theory 

(Fundamental) 

• Results in a Single Rollup 

• Case II: Subharmonic of the Most Unstable Mode 

• Results in a Single Rollup, of Approximately Twice the Size 
and Time Period 

• Case III: Fundamental and Subharmonic 90° Out-of-Phase 

• Results in a Double Rollup 


EQUATIONS (CONT.) 

• Define Conserved Scalar 0 = C A ~ C B 

Then 0 satisfies 

0 u-V0 = Q (No R dependence) 

dt dxf 

with initial conditions 0(t = O) = 0(0) = C^(0) — Cq(0) 

• Infinite Reaction Rate Limit 

(0 0^0 _\ 0 0^0 
1 - |o 0<o 2 ~|-0 0<o 

Cp can be obtained from conserved scalar 
0 = Cp + C^y 

where 0 also satisfies Equation 0 

with initial conditions 0(t = O) = 0(0) = C^(0) 
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Figure 28a t • 12 


Figure 28c t « 36 


Plots of Concentration Contours for Species A for a Sequence of 
Times for Case 3 (Fundamental and Subharmonic Added Together 
Out-of-Phase) - Two-Dimensional Simulations 


Plots of Concentration Contours for Species A and B at 

t * 12 for Case 1 (Fundamental Mode Alone) - Two-Dimensional 

Simulations 








Total Product Versus Time for Two Different Values of the 
Reaction Rate Coefficient R “ One~D linens lona 1 Simulations. 
Also Shown Is the Infinite Reaction Rate Solution 
[Equation (37)] 


PRODUCT 
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Plot of Product Concentration Contours at t s 24 for Case 3 
(Fundamental and Subharmonic Added Together Out-of-Phase) 

- Two-Dimensional Simulations 


Figure 45c t ■ 18 


Plots of Lateral (y) Vorticity Contours for a Sequence cl Times 
in an x-z Plane Located a Lateral Distance L y /2 from that 
Depicted in Figure UU - Single Rollup Case, Three-Dimensional 
Simulations, Run JN20B 
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MfcAN PRODUCT TH1CKNFSS 


r 



Mean Produce Thickness, zp, Versus Time - Three-Dimens iona 
Sinulat ions 


NON DI HENS 10NAL PRODUCT TIIICRNF.SS 


i 




Nondimensional Produce Thickness Versus Time - 
Three-Dimensional Simulations 


a 


Average of the Product of the Concentrations, the Product of 
the Averages, and the Concentration Correlation Versus z at 
t • 12 - Realization 1, Double Rollup Case, Three-Dimensional 
Simulations 










CONCLUSIONS 


• Color Problem Tests Showed That 

• Spectral methods produced high accuracy solutions in the 
presence of advection, reaction, and diffusion 

• Weak spatial filter insured numerical stability 

• One Dimensional Simulations Showed That 

• Calculations close to the infinite reaction rate can be done 

• Results were fairly insensitive to initial conditions 

• Results were in reasonable agreement with theoretical 
predictions 

• Two Dimensional Simulations Showed That 

• Reaction zone was significantly stretched by vortex rollup 

• Vorticity and product fields coincided 

• Significant species segregation and "flame shortening" 
effects were present 


CONCLUSIONS (CONT.) 

• Three Dimensional Simulations Showed That 

• Good agreement with self-similarity theory was obtained 

• Technique can treat entire range of reaction rates - from very 
slow to infinite 

• Computed profiles qualitatively similar to profiles from 
laboratory data 

• Quantitative comparisons will be made soon 

• Nondimensional product thickness was 0.22 - compared to 
values of 0.3 to 0.35 from laboratory data 

• No adjustable parameters are used 

• Some presently used estimates of correlation products 
can be significantly in error 
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CHARACTERISTICS OF INHOMOGENEOUS JETS IN CONFINED SWIRLING AIR FLOWS 

Ronald M. C. So and S. A. Ahmed 
Arizona State University 

Flow fields In the neighborhood of confined jets proved to be of great 
interest for gas turbine combusto’" designers. Although there have been 
numerous investigations on the characteristics of confined (1-4) jets, very 
little quantitative data exist for confined jets In a swirling flow. We have 
recently started an experimental program to study the characteristics of 
Inhomogeneous jets In confined swirling flows to obtain detailed and accurate 
data for the evaluation and Improvement of turbulent transport modeling for 
combustor flows. Our work was also motivated by the need to investigate and 
quantify the Influence of confinement and swirl on the characteristics of 
inhomogeneous jets. 

The flow facility at Arizona State University was constructed In a simple 
way which allows easy Interchange of different swlrlers and the freedom to 
vary the jet Reynolds number. The velocity measurements were taken with a one 
color, one component DISA Model 55L laser-Doppl er anemometer employing the 
forward scatter mode. Standard statistical methods are used to evaluate the 
various moments of the signals to give the flow characteristics. 

The current project started with the measurements of the velocity field. 
The second phase of the work will concentrate on the Investigation of the 
scalar fluxes utilizing concentration probes. Finally, the two different 
sensors will be used simultaneously to determine the turbulent momentum and 
mass fluxes In the whole flow field. 

The present work was directed at the understanding of the velocity field. 
Therefore, only velocity and turbulence data of the axial and circumferential 
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components are reported for Inhomogeneous jets in confined swirling air flows. 
Results to date show that the jet centerline velocity decreases rapidly in a 
short distance for both Helium and air jets. However the similarity between 
Helium and air jets ends here. For air jets, the jet-like behaviour In the 
flow disappears at about 20 diameters downstream of the jet exit. This 
phenomenon Is independent of the Initial jet velocity. When this stage is 
reached for the mean flow# the turbulence field also decays to that of the 
background swirling flow. For Helium jets, the jet-like behaviour Is noticed 
even at 40 diameters downst. earn of the jet exit. The turbulence field also 
reflects the same behaviour. Since the jets are fully turbulent (therefore. 
Independent of jet Reynolds numbers) and the jet momentum fluxes for both air 
and Helium jets are the same, the cause of this difference In behaviour is 
attributed to the combined action of swirl and density difference. This 
behaviour could explain some of the common observations In gas turbine 
combustors. 

The completion of the proposed work will make a substantial contribution 
to the understanding and predictive capability of complex turbulent swirling 
flows with Inhomogeneous jets. 
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OBJECTIVE 


PROGRAM OF WORK 


• MIXING IN COMBUSTOR 

• ISOTHERMAL VARIABLE DENSITY 
MIXING 

• SWIRL EFFECTS 

• SCALE INFORMATION FOR 
COMBUSTOR DESIGN 


• EXPERIMENTAL STUDIES 

I - AIR JETS 

2- HELIUM JETS 

3- CONCENTRATION AND FLUX MEASURE 
MENTS 

4- EFFECTS OF VARYING SWIRL 


• DATA FOR MODEL VERIFICATION 


EXHAUST 

DUCT 



Schematic Oiaqram of the Flow Syston 
(Note that the laser axis is perpendicular to the plane of the paper) 





; 
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HELIUM OR AIR 
SUPPLY 
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Schanatic Diagram Showing the Test Section In Detail 


MOTOR 



Block Diagram of the I DV and Data Acquisition System 
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CP POOR Q JALIlY 


FLOW CONDITIONS 


R-* 5 5* 10 


QUANTITY 

MEASURED 

AIR 

JET 

HELIUM JET 

Uj(mf/*ec) 

25 4 

66 8 

25 4 

66 8 

R «, 

14 x I0 4 

3 8 * I0 4 

18 x |0 S 

4 8 x I0 3 

M ( /M 

0 071 

0 491 

0010 

0.068 


Uj • 770 


Check on the Flow 'ymmetry 


U ( -66 8 
X/Oj»l 3 

. SEEDING UPSTREAM 
FLOW ONLY - 2< 

• SEEDING UPSTREAM 
FLOW AND JET TOO 


Check on the Seed Effects 
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(a) Centerline Velocities at Different Flow Conditions 
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GRAVITATIONAL BUOYANCT - g • 

CENTRIFUGAL ACCELERATION - 

_ (19 m/s) 2 

6.2 x 10 -2 » 


9.8 m/s 2 

i 2 

R 


; 6 x 10^ m/s 2 


RATIO = 600 


DISTANCE FOR He TO RISE ONE R 


x - Uj t 

SINCE ( P H ♦ | * 9 ( P A * P h) 


IF 


Pa 1 7P, 




1/2 


x . Ui IR/.Tg)^ 2 

BJ Bp 


= 800 


EFFECTS OF ROTATION ON MIXING 

FLUID PARTICLES ARE IN RADIAL EQUILIBRIUM DUE TO 
RADIAL PRESS. GRAD ■ CENTRIFUGAL FORCE 


PROMOTE MIXING 



, B 

ft -0 i 

DISPLACED FLUID B 
CONTINUES TO MOVE 
TOWARDS FUEL SURFACE 


IMPEDE MIXING 



- A 

il-0 1 

DISPLACED FLUID A 
WILL BE FORCED BACK 
TOWARDS 1 
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Centerline Velocities at Different Flow Conditions 



Axial Distributions of / u ,2 /V at Different Flow Conditions 
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SATURATED LASER FLUORESCENCE IN 

TURBULENT SOOTING FLAMES AT HIGH PRESSURE 

Galen B. King 
Campbell D. Carter 
Normand M. Laurendeau 

Flame Diagnostics Laboratory 
School of Mechanical Engineering 
Purdue University 


The primary objective of this project is to develop a quantitative, single 
pulse, laser-saturated fluorescence (LSF) technique for measurement of radical 
species concentrations in practical flames. The species of immediate interest is 
the hydroxyl radical. Measurements will be made in both turbulent premixed and dif- 
fusion flames at pressures between 1 and 20 atm. Interferences from Mie scattering 
will be assessed by doping with particles or by controlling soot loading through 
variation of equivalence ratio and fuel type. The efficacy of the LSF method at 
high pressure will be addressed by comparing fluorescence and absorption measure- 
ments in a premixed, laminar flat flame at 1-20 atm. 

Signal-averaging over many laser shots is sufficient to determine the local 
concentration of radical species in laminar flames. However, for turbulent flames, 
single pulse measurements are more appropriate since a statistically significant 
number of laser pulses is needed to determine the probability distribution function 
(PDF) . PDFs can be analyzed to give true average properties and true local kinetics 
in turbulent, chemically reactive flows. 


PROGRESS OF WORK 


The construction phase of the high-pressure facility is nearing completion. 

This includes the fabrication of not only the high-pressure vessel, but also the 
optical table in which the pressure vessel is mounted and other support systems, 
including the high pressure burners and the gas flow control system. 

The pressure vessel is constructed from a 109 cm long by 20 cm ID stainless 
steel tube, to which four 6.4 cm ID tubes are welded. Slip-on flanges are welded at 
the ends of these smaller tubes to hold the optical windows. Central to the design 
of the pressure vessel is the use of an internal burner-translation system rather 
than translation of the entire vessel to access the various points in the flame. 

This arrangement allows two dimensional movement of the burner within the pressure 
vessel. Maximum travel is greater than 11 cm in the vertical and 3.8 cm in the 
horizontal directions. The translation stage is driven by stepper motors which are 
in turn controlled by an Apple computer. The position of the burner is monitored 
with two linear transducers and a pressure transducer monitors the pressure within 
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the vessel. By counting the number of steps sent to the stepper motors, the posi- 
tion of the burner may be controlled to better than 3 pm in both travel directions. 

Because of its size and mass, the pressure vessel will be located between two 
optical bread boards rather than on top of a single optical table. A framework to 
support the pressure vessel has been added to an existing table base and the optical 
bread boards have been mounted on the table framework. Additional modifications to 
the table are now underway, including (1) a cylindrical plexiglas shield to protect 
laboratory personnel in case of pressure failure and (2) steel plates which will 
surround the vessel to form a flat surface between the bread boards. A small winch 
is also being mounted in the ceiling of the laboratory to aid in movement of the 
pressure vessel during experiments. 

Although initial plans for the high-pressure flame facility called for control 
of gases by mass flow transducers, a simpler gas flow control system will be used 
at this time. This strategy will allow greater variation in burner type during the 
exploratory stage of our research. The control system we have constructed uses 
ordinary regulators and rotometers. This system, while simple and easily varied, 
is limited to a useful pressure range of 10 atm. Therefore, we plan to eventually 
add mass flow controllers to the flame facility. 

Two burners have been designed for initial testing of the high pressure flame 
facility. The first is a premixed water-cooled burner constructed with a Hasteloy, 
hexagonal -pattern, flame holder. The main body of this burner has been constructed 
and work is now in progress to shape the flame holder. A second flat flame burner 
is being built by McKenna Products. This burner is a scaled down version of their 
standard water-cooled sintered bronze, flate flame burner. The flame diameter for 
both burners is approximately 3 cm to avoid self-absorption problems at high 
pressure. 

Significant advances have also been made in our ability to measure single-pulse 
OH radical concentrations in turbulent flames. Measurements of superequilibrium 
hydroxyl concentrations in nonpremixed flames at 1 atm using LSF have been made in a 
collaborative interaction between the Flame Diagnostics Laboratory at Purdue Univer- 
sity and Corporate Research and Development at General Electric Company. Tha re- 
sults of this study were presented at the 1983 Fall Meeting of the Western States 
Section of the Combustion Institute and a paper has been submitted to the Twentieth 
International Symposium on Combustion. 

An improvement has also been made in our sampling technique as we can now dis- 
criminate against bad laser shots. With the use of a second sampling oscilloscope, 
the power level of the laser pulse is monitored and fed into a computer along with 
the fluorescence signal. We are currently using this technique for measurements of 
PAHs and will eventually incorporate it into our procedure for generating PDFs. In 
this method, the analogue-to-digital converter and both sampling oscilloscopes are 
triggered by the laser pulse and thus single pulse data for both the laser power 
and the fluorescence signal can be collected. Data points in which the laser pulse 
shape or energy level are not satisfactory are then eliminated from the data set. 
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LEAN LIMIT PHENOMENA 
C. K. Uw 

Northwestern University 


Extensive theoretical and experimental research have been accomplished on the 
various aspects of the structure and extinction of preraixed flames. The 
presentation will focus on only one specific topic of the lean/rich limit phenomena, 
namely the concept of flammability limits in the presence of flame interaction, and 
the existence of negative flame speeds. Attention is called to other topics 
reported in the publications listed; contributions from them are no less significant 
but will not be presented due to the lack of time. 

To appreciate the importance of flame interaction, we first note that although 
predictions of complex flame phenomena are frequently based on understanding of 
isolated flames, it is clear that in most situations the combustion flow field is 
composed of an ensemble of fiamelets, of different intensity and extent, which 
continuously interact with each other and thereby can cause significant 
modifications of the bulk combustion behavior. 

In the present invest Ration downstream interaction between two counterflow 
preraixed flames of different stoichiometries are experimentally studied. Various 
flame configurations are observed and quantified; these include the binary system of 
two lean or rich flames, the triplet system of a lean and a rich flame separated by 
a diffusion flame, and single diffusion flames with some degree of preraixedness . 
Extinction limits are determined for methane/air and butane/air mixtures over the 
entire range of mixture concentrations. 

The results show that the extent of fJame interaction depends on the separation 
distance between the flames which are functions of the mixtures' concentrations, the 
stretch rate, and the effective Lewis numbers (Le). In particular, in a positively- 
stretched flow field Le<l (>l) mixtures tend to interact strongly (weakly), while 
the converse holds for flames in a negatively-stretched flow. 

An important consequence of flame interaction is the modification of the 
flammability limits of a mixture when it is stratified. If the mixture is only 
weakly stratified, then the conventional flammability limits for uniform mixtures 
may still apply. However, with more extensive stratification the stronger portion 
of the mixture may burn Independent of the weaker portion, while the weaker portion 
can also burn by receiving support from the stronger flame. The mixture is made 
more flammable as a whole. 

Our study has also established the existence of negative flames whose 
propagation velocity is in the same general direction as that of the bulk convective 
flow, being supported by diffusion alone. Their existence demonstrates the tendency 
of flames to resist extinction, and futher emphasizes the possibility of very lean 
or rich mixtures to undergo combustion. 
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Flame interaction is also of relevance to the modeling of turbulent flames. 
Within a turbulent flow field stratifications in temperature and reactant 
concentrations invariably exist, caused by mixing inhomogeneities or, more subtly, 
through the coupling between stretch and preferential diffusion even for an 
initially homogeneous mixture. Thus a flame kernel initiated within a turbulent 
eddy will manifest effects of downstream inteaction between flames of different 
intensities. 

Publications 


1. "Effects of Heat Loss, Preferential Diffusion, and Flame Stretch on Flame-Front 

Instability and Extinction of Propane/Air Mixtures," by S. Ishizuka, K. 
Miyasaka, and C. K. Law, Combustion and Flame , Vol. 45, pp. 293-308, (1982). 

2. "On the Opening of Premixed Bunsen Flame Tips," by C. K. Law, S. Ishizuka, and 

P. Cho, Combustion Science and Technology , Vol. 28, pp. 89-96, (1982). 

3. "On Stability of Premixed Flames in Stagnation-Point Flow," by G. I. 

Sivashinsky, C. K. Law, and G. Joulin, Combustion Science and Technology , 
Vol. 28, pp. 155-159, (1982). 

4. "An Experimental Study of Extinction and Stability of Stretched Premixed 

Flames," by S. Ishizuka and C. K. Law, Nineteenth Symposium on Combustion , 
the Combustion Institute, Pittsburgh, Pa., pp. 327-335, (1983). 

5. "Heat and Mass Transfer in Combustion: Fundamental Concepts and Analytical 

Techniques," Proc. of ASME-JSME Joint Thermal Engineering Conference , (Y. 
Mori and W. J. Yang, Ed.), Vol. 2. pp. 535-559, (1983). ** Plenary Paper** 

6. "An Invariant Derivation of Flame Stretch," by S. H. Chung and C. K. Law, 

Combustion and flame , Vol. 55, pp. 123-125 (1984). 

7. "Extinction of Premixed Flames by Stretch and Radiative Loss," by S. H. Sohrab 

and C. K. Law, to appear in International Journal of Heat and Mass Transfer . 

8. " On the Determination of Laminar Flame Speeds from Stretched Flames," by C. K. 

Wu and C. K. Law, submitted. 

9. "An Experimental Investigation on Flame Interaction and the Existence of 

Negative Flame Speeds," by S. H. Sohrab, Z. Y. Ye, and C. K. Law, submitted. 

10. "Effects of Preferential Diffusion on the Burning Intensity of Curved Flames," 

by M. Mizomoto, Y. Asaka, S. Ikai, and C. K. Law, submitted. 
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OVERALL OBJECTIVES 


1 


o TO UNDERSTAND THE FUNDAMENTAL 

MECHANISMS GOVERNING EXTINCTION OF 
PREMIXED FLAMES. 

o TO QUANTIFY FLAMMABILITY LIMITS OF 
MIXTURES AND IDENTIFY DOMINANT 
MECHANISMS IN OPERATION AT THESE 
LIMITS. 

o TO IDENTIFY POSSIBILITIES OF 
EXTENDING THESE EXTINCTION/ 
FLAMMABILITY LIMITS. 


2 

POTENTIAL EXTINCTION MECHANISMS 

o HEAT LOSS. 

o KINETIC TERMINATION. 

o FLAME STRETCH: 

FLOW NONUNIFORMITY 
FLAME CURVATURE 
FLAME ACCELERATION 


2 

PREVIOUS CONTRIBUT IONS 

o FLAME STRETCH, COUPLED WITH 

PREFERENTIAL DIFFUSION (NON-UNITY 
LEWIS NUMBER) EFFECTS, WAS IDENTIFIED 
TO BE THE DOMINANT MECHANISM IN 
OPERATION AT THE FLAMMABILTY LIMITS. 

o FLAMMABILITY LIMITS, DETERMINED BY 
USING THE COUNTERFLOW FLAME, AGREE 
WELL WITH EXISTING DATA. 


2 

PRESENT ACCOMPLISHMENTS 

o IDENTIFIED AND QUANTIFIED THE 
IMPORTANCE OF FLAME INTERACTION. 

o INDENTIFIED THE EXISTENCE OF LAMINAR 
FLAMES WITH NEGATIVE FLAME SPEEDS. 

o STUDIED EFFECTS OF FLAME CURVATURE ON 
BURNING INTENSITY. 

o CRITICALLY RE-EXAMINED EXISTING 
METHODS IN LAMINAR FLAME SPEED 
DETERMINATION; PROPOSED NEW 
METHODOLOGY. 

o THEORETICAL STUDIES ASSOCIATED WITH 
THE ABOVE TOPICS. 


_5 

IMPORTANCE OF F ' AME INTERACTION 

o PRACTICAL COMBUSTION FLOW FIELDS 
CONSIST OF ENSEMBLES OF FLAMELETS. 

O INDIVIDUAL FLAMELETS BURN WITH 
DIFFERENT INTENSITY AND EXTENT. 

o INTERACTION BETWEEN THEM CAN EITHER 
CONTRACT OR WIDEN THE EXTINCTION 
LIMIT. 

o ESPECIALLY RELEVANT TO TURBULENT 
FLAME MODELING AND FLAME PROPAGATION 
IN FLOWS WITH CONCENTRATION 
STRATIFICATION OR POOR MIXING 
CHARACTERISTICS. 
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METHODOLOGY 


6 


8 


o ESTABLISH TWO PREMIXED FLAMES OF 
UNEQUAL STRENGTH BY THE COUNTERFLOW 
OF 'TWO STREAMS WITH DIFFERENT MIXTURE 
CONCENTRATIONS. 

O MAP THE FLAME CONFIGURATIONS AND 
EXTINCTION LIMITS AS FUNCTIONS OF 
STRETCH AND FUEL CONCENTRATIONS OF 
THE UPPER (SI.) AND LOWER ( ft ) 

STREAMS. 

o MEASURE THE TEMPERATURE PROFILE FOR 
THE INTERACTING FLAMES. 


]_ 

EXTINCTION LIMITS (METHANE-AIR 
MIXTURES) 

o FOR SI L < 4%, S^ » 6% AND IS 

INSENSITIVE TO ft^; WEAK INTERACTION; 
STRONG FLAME BURNS INDEPENDENTLY, 

WEAK FLAME PARASITIC. 

o FOR 4% < (ft L , S^) < 6%, 

(S2 L + ft^/2 - 5%; STRONG INTERACTION; 
SYMBIOTIC COMBUSTION. 

o INTERACTION IS STRONG (WEAK) FOR LEAN 
(RICH) MIXTURES. 


PRACTICAL CONCLUSIONS 

o METHANE-AIR MIXTURES LEANER THAN 5% 
REMAINS NON-FLAMMABLE IF 
CONCENTRATION STRATIFICATION (OR 
MIXING INHOMOGENEITY) IS WITHIN 4-6%; 
A HIGHER STRATIFICATION WILL INDUCE 
BURNING. 

o MIXTURES ABOVE 6% WILL ALWAYS BURN, 
HENCE AUGMENTED FLAMMABILITY LIMIT. 

o INTERACTION IS STRONG (WEAK) FOR LEAN 
(RICH) MIXTURES. 

o FLAMMABILTY LIMITS LESS (MORE) 
SUSCEPTIBLE TO CONCENTRATION 
STRATIFICATION FOR LEAN (RICH) 
MIXTURES. 

o PROPANE-AIR AND BUTANE-AIR MIXTURES 
BEHAVE OPPOSITE TO METHANE-AIR 
MIXTURES. 

o MIXING INHOMOGENEITY PROMOTES 
BURNING. 

o EXISTENCE OF NEGATIVE FLAME SPEEDS 
DEMONSTRATES THE RESISTENCE OF FLAMES 
TO EXTINCTION. 

o flame: stretch and preferential 

DIFFUSION ARE CRITICAL FACTORS IN 
DETERMINING FLAM. CHARACTERISTICS, 
ESPECIALLY EXTINCTION. 

o FRACTIONAL DEVIATION OF LEWIS NUMBER 
FROM UNITY CAN CAUSE SAME EXTENT OF 
DEVIATION OF FLAME TEMPERATURE. 


o THEORY SUBSTANTIATES THE ABOVE 
BEHAVIOR. 
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Figure i. Possible modes of downstream interaction between two 
premixed flames. PF : Premixed Flame; DF : Diffusion 
Flame; LPF/RPr: Lean/Rich Tremixed Flame; PPF : Partially 
Premixed Flame; PDF- Partial Diffusion Flame. 




Figure 2. Mapping of the flame 

configurations and extinction 
boundaries tor methane/air 
mixtures . 


Figure 3. Mapping of the flame 

configurations and extinction 
boundaries for butane/air 
mixt ores . 
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Figure 4. Temperature profiles of methane/air flames: (a) Symmetric flames; 

(b) Asymmetric flames; (c) Negative flame; (d) Triple flames; (e) 
l. Diffusion flame; 2. Partial diffusion flame, concentrations are 
fuel + 10% stoich. in upper flow and air + 10% stoich. fuel in 
lower flow. 
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Figure 5. The estinctlon boundary for two lean 
premised methane/alr flames. 


Figure 6. Distance between lover flame and stagnation 
surface, f, as a function of Cl . 
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INTERACTIVE COMPUTER MODELING OF COMBUSTION CHEMISTRY 

and 

COALESCENCE-DISPERSION MODELING OF TURBULENT COMBUSTION 

David T. Pratt 
University of Washington 


The goals of this research project are as follows: 

1. Develop an interactive computer code for simulation of a high-intensity 
turbulent combustor as a "single point" inhomogeneous stirred reactor [l] . This 
will be developed from an existing batch processing computer code CDPSR [2]. 

2. Use the interactive CDPSR code as a guide for interpretation and direction 
of DOE-sponsored companion experiments utilizing Xenon tracer with optical laser 
diagnostic techniques to experimentally determine the appopriate mixing frequency, 
and for validation of CDPSR as a mixing-chemistry model for a laboratory jet- 
stirred reactor. 

3. Incorporate the coalescence-dispersion model for finite rate mixing into an 
existing interactive code AVCO-MARK I, to enable simulation of a combustor as a 
modular array of stirred flow and plug flow elements, each having a prescribed 
finite mixing frequency, or axial distribution of mixing frequency, as appropriate. 

4. Further increase the speed and reliability of the batch kinetx^s integrator 
code CREKID [3] by rewriting in vectorized form for execution on a vector or parallel 
processor, and by incorporating numerical techniques which enhance execution speed 
by permitting specification of a very low accuracy tolerance [4]. 


REFERENCES 


1. Pratt, D. T.: Mixing and Chemical Reaction in Continuous Combustion, in 

Progress. Progress in Energy and Combustion Science, vol. 1, N. A. Chigier,Ed., 
Pergamon Press, 1976. 

2. Pratt, D. T.: Coalescence/Dispersion Modeling of High Intensity Combustion, 

AIAA J. Energy, 1 , 3, 177180, 1979. 

3. Pratt, D. T. : CREKID: A Computer Code for Transient, Gas-phase Combustion 

Kinetics, Paper No. WSCI 8321, Western States Section/The Combustion Institute, 
1983. 

4. Pratt, D. T. : Exponential-Fitted Methods for Integrating Stiff Systems of 

Ordinary Differential Equations: Applications to Homogeneous Gas-Phase 

Chemical Kinetics, paper presented at the 1984 JANNAF Propulsion Meeting, 

New Orleans, LA., February 1984. 
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*** HARK 2 1 **■* 

•MARK2I* IS AN INTERACTIVE VERSION OF THE MARK-II COMBUSTOR MODEL . 
THIS IS A PRELIMINARY DESIGN TOOL, A MEANS OF GAINING INTUITIVE INSIGHT 
INTO EFFECTS OF CHANGES IN FUEL-AIR MIXING OR PARTITIONING ON TURN-DOWN 
RATIO, COMBUSTION EFFICIENCY AND FQLLUTANT FORMATION RATES. 

AN INITIAL DATA SET IS TAKEN FROM DATA FILE * MAKN2 . DA T * BUT CAN BE ALTERED 
INTERACTIVELY, AND USED IN CONSECUTIVE RUNS. 

MARK-II REPRESENTS A SIMPLE DRAGC, COMBUSTOR CONSISTING OF A 
MAXIMUM OF 9 FLOW ELEMENTS WITH THE ADDITION OF A SINGLE RECYCLE 
ELEMENT. FLOW ELEMENT TYPES MAY INCLUDE ! 

1) NON-REACTING MIXERS (’MIX*), IN WHICH THE CHEMICAL REACTIONS ARE 
ASSUMED TO HAVE STOPPED DURING THE MIXING PROCESS} 

2) PERFECTLY STIRRED REACTORS <*PSR*>, WITHIN WHICH INTENSE SELF- OR 
BACK -MIXING IS ASSUMED TO OCCUR, SO THAT THERE ARE NO AXIAL GRADIENTS} 

3) PLUG FLOW REACTORS ( * F’FR * ) . 

THE USER MAY DEFINE THE MODEL AS HAVING UP TO 9 ELEMENTS IN 
SERIES WITH AIR AND FUEL INLET JETS AT EACH ELEMENT. THE RECYCLE ELEMENT 
MAY BE OF ANY OF THE THREE FLOW TYPES, AND MUST RECYCLE FROM A HIGHER 
NUMBERED ELEMENT TO A LOWER. COOLING BOUNDARY LAYER EFFECTS AND CHEMICAL 
REACTIONS WITHIN THE BOUNDARY LAYER ARE NOT CONSIDERED. 

PLEASE WAIT A MOMENT WHILE INITIALIZATION IS COMPLETED. 

INITIALIZED — PRESS ^RETURN> TO BEGIN — 


*** INPUT 

data »** 
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FAST ALGORITHMS FOR COMBUSTION KINETICS CALCULATIONS: A COMPARISON* 

Krishnan Radhakrishnan** 

NASA Lewis Research Center 


Many practical problems arising in chemically react5.ng flows require the 
simultaneous numerical integration of large sets of chemical kinetic rate 
equations of the type shown in figure 1. The initial value problem is that of 
finding the composition and temperature at the end of a prescribed time 
interval, given the initial mixture composition and temperature, the pressure, 
and the reaction mechanism. Multi-dimensional modeling of reactive flows 
requires the integration of the system of ordinary differential equations 
(ode's) given in figure 1 at several thousand grid points. To make such 
calculations practicable, it is necessary to have a very fast batch chemistry 
integrator . 

To identify the fastest algorithm currently available for the numerical 
integration of chemical kinetic rate equations, several algorithms have been 
examined. In the present paper, we summarize our findings to date — details 
are available in references (1) and (2). The algorithms examined in this work 
include two general-purpose codes EPISODE and LSODE (refs. 3 and 4), and three 
special-purpose (for chemical kinetic calculations) codes CHEMEQ (ref. 5), 
CREK1D (refs. 6 and 7), and GCKP84 (refs. 8 and 9). In addition, an explicit 
Runge-Kutta-Merson differential equation solver (ref. 10) (IMSL Routine DASCRU) 
is used to illustrate the problems associated with integrating chemical kinetic 
rate equations by a classical method. These methods are summarized in figure 3. 

The algorithms summarized in figure 3 were applied to two test problems 
drawn from combustion kinetics. These problems, summarized in figure <1, 
included all three combustion regimes: induction, heat release and 
equilibration. Figures and 6 present variations of the temperature and 
species mole fractions with time for test problems 1 and 2, respectively. Both 
test problems were integrated over a time interval of 1 ms in order to obtain 
near-equilibration of all species and temperature. 

Of the codes examined in this study, only CREK1D and GCKP84 were written 
explicitly for integrating exothermic, non- isothermal combustion rate 
equations. These therefore have built-in procedures for calculating the 
temperature (T). For the other codes, two different methods, labeled as Methods 
A and B, oe re used to compute T. The following convention was adopted in naming 
these other cedes: those using temperature method A were given the suffix-A 
(e.g. LSODE-A) and those using temperature method B were given the suffix-B 
(e.g. LSODE-B). In Method A, T was calculated from the mole numbers and the 
initial mixture enthalpy using an algebraic energy consei vation equation (given 
in figure 7) and a Newton-Raphson iteration technique. In this method, the 
temperature is not an explicit independent variable, so the number of ode's is 


* Work partially funded by NA'A Grants NAG3-1'17 and NAG3-294. 

** NRC-NASA Research Associate: on leave from The University of Michigan, Dept, 
of Mechanical Engineering and pplied Mechanics, Ann Arbor, MI <18109. 
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equal to the number, NS, of distinct chemical species in the mixture. The 
integrator therefore tracks only the solutions for the species mole numbers. In 
Method B, the temperature was treated as an additional independent variable and 
evaluated by integrating its time-derivative given in figure 7. In this method, 
the number of ode's is equal to NS+1, and the integrator tracks the solutions 
for both the temperature and the species mole numbers. 

All codes were run on the NASA Lewis Research Center's IBM 370/3033 
computer using single-precision accuracy, except GCKP84 which was in double- 
precision. A typical computational run consisted of initializing the species 
mole numbers, temperature, and CPU time. The integrator was then called with 
values for the necessary input parameters. On return from the integrator, the 
total CPU time required to solve the test problem was calculated. Other 
performance parameters were also recorded — see reference (2) for details. 

Figures 8 and 9 present the computational work (expressed as the CPU time 
in seconds) plotted against the local error tolerance, EPS, for test problems 1 
and 2, respectively. For all codes except EPISODE, EPS is the local relative 
error tolerance. For EPISODE, EPS is a mixed error tolerance — relative for 
species with initially nonzero mole numbers and for the temperature (method B) 
and absolute for species with initially zero mole numbers. Also shown on 
figures 8 and 9 are the CPU times required by the explicit Runge-Kutta method 
for one value of EPS. Note the excessive CPU times required by this technique. 
Its use would make multidimensional modeling of practical combustion devices 
prohibitively expensive. 

For test problem 1, very small values for EPS had to be used for EPISODE 
(figure 8). For values of EPS 2 5 x 10"*, EPISODE predicted little or no change 
in the composition and temperature after an elapsed time of 1 ms. Similar 
remarks apply to test problem 2 (figure 9), for which values of 10" 1 * and 10“° 
had to be used for EPISODE-A and F.PISODE-B, respectively. Although the runs 
with EPISODE-B a.Oi EPS £ 5 x 10" 1 * were successfully completed, the solutions 
(especially for minor species) were significantly different from those given in 
figure 6. With GCKP84 and EPS = 10" 2 , the solution for test problem 1 exhibited 
serious instability and so was terminated. A more detailed discussion of the 
accuracy of the codes tested in this study can be found in reference (2). 

Examination of figure 8 shows that the difference in computational work 
required by methods A and B is small for test problem 1, with method B being 
more efficient. For test problem 2 (figure 9), the difference is small for 
large values of EPS. But for small values of EPS the difference is more marked, 
with method A being significanlty superior to method B. 

Figures 8 and 9 show that LSODE and CREK1D are superior to the other 
codes. EPISODE is an attractive alternative, especially for test problem 2. 
However, in using EPISODE, a word of caution is in order. The computational 
work can bo strongly dependent on the value for the initial steplength (HO) 
selected by the user. A poor guess for HO can make EPISODE pronibitively 
expensive to use. Figure 10 illustrates this behavior for test problem 2. Note 
an order of magnitude increase in the CPU time for a change in HO from 10“ 7 to 
10"* s. Although not shown here, a poor guess for HO also resulted in 
inaccurate and unstable solutions. In addition, as discussed in reference (2), 
the error control performed by EPISODE is unsatisfactory for problems of the 
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type examined in this study. 

A simple method for increasing the efficiency of the algorithms as applied 
to the present problem was explored. This involved updating the rate constants 
kj and k_j (which was calculated from kj and the concentration equilibrium 
constant) only for temperature changes greater than an amount AT. To avoid a 
trial and error search for the optimum value of AT — defined as that value 
which results in minimum computational work — an approximation for it was 
derived and is presented in figure 11. Comparisons of figure 12 with figures 7 
and 8 show the significant reductions in computational work realized by use of 
the above approximation for AT. 
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Adiabatic, Constant-Pressure, Gas-Phase Chemical Reaction 
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In the above equations. 



molar rate of formation of species i per unit mass of mixture, 
kmol -i/kg-mixture s 

forward and reverse rate constants for reaction j 
mole number of species i, kmol-i/kg-mixture 

pre-exponential constants in forward and reverse rate equations for 
reaction j 

activation energy in forward and reverse rate equations for reaction j, 
cal /mole 

number of distinct elementary reactions in mechanism 
number of distinct species in gas mixture 
universal gas constant, 1.987 cal/mol K 

forward and reverse molar reaction rates per unit volume for reaction j, 

kmol/m 3 s 

temperature, K 

mixture mass-density, kg/m 3 

stoichiometric coefficients of species i in reaction j as a reactant, 
and as a product, respectively 


Figure 1 Governing Ordinary Differential Equations 
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Given 


Initial Mixture Composition and Temperature 
Pressure 

Reaction Mechanism 

Find, at the End of a Prescribed Time Interval 
Mixture Composition and Temperature 

Figure 2 Problem Statement 


Me thod 

. 

Description 

CCllPfl'J 

Details not yet available. 

CSEK1D 

Variable-step, predictor-corrector method based on an 
exponentially-fitted trapezoidal rule; includes filtering of 
ill-posed initial conditions and automatic selection of 
Jacobi-Newton iteration or Newton iteration. 

LSODE 

EPISODE 

Variable-step, variable-order backward-differentiation 
method with a generalized Newton iteration*. 

CIEMEQ 

Vai iable-step, second-order predictor-corrector method with 
an asymptotic integration formula for stiff equations. 

DJVSCRU 

Variable-step, fourth-order, explicit Runge-Kutta-Merson 
solver. 


*Othor options are included in these packages. 


Figure 3 Summary of Methods Studied 
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Two Problems Describing Adiabatic, Constant Pressure Chemical Reactions 


Test Problem 1: 


Combustion of a Mixture of 33% CO and 67% H 2 with 100% Theoretical Air 
(taken from reference 11) 

12 Reactions 

11 Species + Temperature 

Pressure = 10 atm. 

Initial Temperature = 1000 K 
Reaction Duration: 1 ms 


Test Problem 2: 


Combustion of a Stoichiometric Mixture of H ^ and Air 
(taken from reference 9) 

30 Reactions 

15 Species + Temperature 

Pressure = 2 atm. 

Initial Mixture Temperature = 1500 K 
Reaction Duration: 1 ms 


Both Problems Include All Three Regimes of Combustion: 
Induction, Heat Release and Equilibration. 


Figure 4 Test Problems 
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Method A 


For adiabatic, constant-pressure combustion reaction, energy conservation gives 


NS 


z 

i = 1 


n. h. 
1 1 


h = constant 

0 


( 1 ) 


where, h ^ = molal-specif ic enthalpy of species i (J/kmol) 
and h Q = mass-spccif ic enthalpy of mixture (J/kg) 

In this method, equation (1) was solved for the temperature using a Nevrton- 
Raphson iteration technique. 


Method B 


Differentiation of eq. (1) with respect to temperature (T) gives 


dT 

dt 


NS 

r 

i =1 


rf.h. 


NS 

I n. c 
1 =1 


where, c n is the constant-pressure 


In this method, the temperature was 


( 2 ) 


specific heat of species i (J/kmol K). 
evaluated by integrating equation (2). 


Figure 7 Evaluation of Temperature 
(for LSODE, EPISODE, and CHEMEQ) 
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GCKP84 

CREK1D 

LSODE-A 

LSODE-B 

EPISODE-A 

EPISODE-B 

CHEMEQ-A 

CHEMEO-B 

DASCRU-A 

OASCRU-B 


r\j 

cr 



a_ 

o 



Figure 8 - Variation of the CPU time 1st with error tolerance, 
EPS, for test problem 1. All runs on IBM 370/3033. 
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— • — LSODE-B 

— -O — EPISODE-A 
» — EPISODE-B 
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HO (s) 

CPU (s) 

io- B 

0.786 

1 

o 

H 

0.703 

io - 7 

0.791 

io-* 

7.91 

io-» 

8.04 

10-*° 

0.772 


Figure 10 Example of Effect of Initial Steplength (HO) 
on Work Required by EPISODE-A (EPS = 10“ B ) 
for Test Problem 2 


An approximate expression for AT — the maximum allowable temperature 
change allowed before the reaction rate constants k- and k . are updated — was 
derived by requiring that the maximum relative error in the resultant reaction 
rates does not exceed the local relative error tolerance (EPS) required 
of the numerical solution. The approximation for AT is given by 


AT 


EPS-T 



E . 

E . 

max 

+ N 

. _lL 

j 

RT J 

’ RT 


(3) 


where, T is 
the maximum 


the current temperature, the bars | 
is taken over all forward and reverse 


denote absolute value, 
reactions. 


and 


Figure 11 Approximation for AT 
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°WGh-,-;^ 

Cf poor 


QJALi j*Y 


Method 

Test Problem 

1 2 

GCKP04 

0.85 

1.73 

CREK1D 

0.23 

1.04 

LSODE-A 

0.31* 

0.52* 

LSODE-B 

0.29* 

0.51* 

EPISODE-A 

0.75* 

0.54* 

EPISODE-B 

0.70* 

0.67 

CHEMEQ-A 

8.-11* 

13.6* 

CUEMEQ-B 

5.69* 

12.3* 


* method incorporated eq. (3) 


Figure 12 Minimum CPU Time (in seconds on IBM 370/3033 computer) 
Required for the Test Problems 
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THE ROLE OF SURFACE GENERATED RADICALS IN CATALYTIC COMBUSTION* 

D.A. Santavicca, Y. Stein and B.S.H. Royce 
Princeton University 


The role of surface generated OH radicals In determining the catalytic ignition 
characteristics for propane oxidation on platinum Is under study. The experiments are 
being conducted In a stacked-plate, catalyst bed. Transient measurements, during 
catalytic ignition, of the catalyst's axial temperature profile have been made and the 
effect of equivalence ratio, inlet temperature and Inlet velocity are being investi- 
gated. These measurements will provide insights which will be useful In planning and 
interpreting to OH measurements. Attempts to measure OH concentration in the catalyst 
bed using resonance absorption spectroscopy have been unsuccessful, indicating that OH 
concentrations are below 10^/cc but still possibly above equilibrium values. 
Measurements are currently underway using forward scatter laser induced fluorescence 
which should extend the OH detection limits several orders of magnitude below the 
equilibrium concentrations. 


INTRODUCTION 


The advantages of catalytic combustion over conventional combustion are lower 
emissions, higher efficiency, increased operational stability, stable operation at 
lower equivalence ratios, improved pattern factors, and wider fuel specifications. 

The performance of a catalytic combustor is determined by the interaction of a number 
of physical and chemical processes, including convection and diffusion of species, 
heat and momentum; conduction and radiation heat transfer in the substrate; and gas 
phase and surface chemical reactions. An important, and as yet unresolved, question 
regarding the interaction of these fundamental physical and chemical processes is the 
effect of catalytic surface reactions on the gas phase chemistry. In particular, it 
is the effect of intermediate and radical species, generated by catalytic wall reac- 
tions, on the gas phase chemistry which has not been established. That radicals can 
be produced by catalytic surface reactions and that at sufficiently high surface tem- 
peratures the radicals are desorbed by the surface into the gas phase has been well 
documented in low pressure flow experiments [1-5]. However, the importance of this 
phenomenon in catalytic combustion and the role of these radicals in determining the 
performance characteristics of catalytic combustors has not been established. 


EXPERIMENTS AND STATUS 


A schematic drawing of the stacked-plate, catalytic combustor in which the 
propane-air catalytic ignition experiments are being conducted is shown in 

* 

Sponsored by NASA-Lewls, Grant NAG 3-353 and AFOSR, Grant 6475. 
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PRECEDING PAGE BLANK NOT FILMED 


Figure 1. The catalyst plates are 50 mm wide, 100 mm long and 1 ram thick. There are 
a total of nine plates in the catalyst bed with an Inter-plate spacing of 6 mm. The 
catalyst plates are made from a steel alloy substrate, coated with aluminum oxide and 
then platinum. Six 0.75 mm diameter holes have been drilled into the side of the 
center plate at six axial locations. Thermocouples have been inserted into these 
holes for measurement of the plate's axial temperature profile during ignition. A 
water cooled gas sampling probe and a thermocouple are used to obtain gas composition 
and temperature measurements at the exit of the catalyst bed. Optical access is 
available through windows in the sides of the test section located at several axial 
positions. The window aperture is 10 mm, therefore the entire distance between two 
plates can be probed. 

Based on a similar experiment by Cattolica and Schefer [6] it was anticipated 
that the OH concentrations would be large enough for detection by resonance absorption 
spectroscopy. Under the conditions of this experiment the OH detection limit using 
this technique is about 10^/cc, as compared to equilibrium OH concentrations of about 
10^/cc (based on the overall equivalence ratio and the catalyst temperature - see 
Figure 2). Measurements were attempted in the stacked-plate combustor using resonance 
absorption and were unsuccessful, implying that the OH concentration is less than 
10^6/cc but still possibly greater than the overall equilibrium concentration. In 
order to extend the OH detection limits it was decided to use laser induced 
fluorescence (LIF). Because of the limited optical access with the stacked-plate con- 
figuration, it has been necessary to use forward scatter collection. The spatial 
resolution with forward scatter collection is comparable to that of the absorption 
measurements, which is adequate for the 2-D flow field between the catalyst plates. 

The forward scatter LIF system is currently being evaluated in a flat flame burner by 
comparing it to simultaneous absorption measurements, after which it will he applied 
to the catalytic combustor. 

Although there have been several [7-9] theoretical studies of catalytic ignition 
there have been no experimental studies reported to date. In order to better 
understand the characteristics of catalytic ignition, especially for the purpose of 
providing insights to guide the OH experiments, it was decided to investigate the 
effects of equivalence ratio, inlet temperature and inlet velocity on catalytic igni- 
tion. Experiments have been conducted where the inlet temperature and velocity are 
fixed, the fuel is turned on and the transient response of the catalyst axial tem- 
perature profile is measured. Preliminary results from the catalytic ignition experi- 
ments are shown in Figures 3-5. 
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AXIAL DISTANCE. X/L 

Normalized temperature versus X/L at different times after Ignition until 
steady state, for 700°K inlet temperature, 0.3 equivalence ratio and 6 
m/sec reference velocity. 
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and 10 m/sec reference velocity. 
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LINER ENVIRONMENT EFFECTS STUDY 


K. S. Venkataramani and E. E. Ekstedt 
General Electric Co. 

Aircraft Engine Business Group 


Estimation of the heat flux to the combustor liner is a key step in the design 
of aircraft engine combustion systems. This forms the basis for determining the 
amount of cooling air and the method of introducing it. Currently, this is largely 
an empirical effort. Future design constraints such as higher pressures and 
temperatures, shorter combustor lengths and tolerance to poorer quality fuels, 
however, accentuate the need for a firmer basis for the heat transfer calculations. 
In particular, it becomes necessary to account for the radiation contributions from 
the flame gases (in spectral bands) and soot particles (continuum) over a wide range 
of combustor operating conditions. Analytical efforts to model the liner heat 
transfer reflecting the above complexities are hampered by the lack of sufficient 
experimental data for model verification. 

The Liner Environment Effects Study Program described here is designed to 
address this need. It is aimed at establishing a broad heat transfer data base 
under controlled experimental conditions by quantifying the effects of the 
combustion system conditions on the combustor liner thermal loading and on the 
flame radiation characteristics. 


Five liner concepts spanning the spectrum of liner design technology from the 
very simple to the most advanced concepts will be investigated. These concepts 
comprise an uncooled liner, a conventional film cooled liner, an impingement /film 
cooled liner, a laser drilled liner approaching the concept of a porous wall and a 
siliconized silicon carbide ceramic liner. The liners will be accommodated in a 
simple test rig housing a three-inch diameter combustor. 


Effect of fuel type will be covered by using fuels containing 11.8, 12.8, and 
14% hydrogen. Tests at 100, 200, and 300 psia will provide a basis for evaluating 
the effect of pressure on the heat transfer. The effects of the atomization quality 
and spray characteristics will be examined by varying the fuel spray Sauter mean 
diameter and the spray angle. Additional parameters to be varied include reference 
velocity, a wide range of equivalence ratio, cooling flow rate, coolant temperature 
and the velocity of the coolant stream on the backside of the liner. 

Both spectral and total radiation measurements will be made in addition to 
obtaining extensive liner metal and film temperature data. 

Reference : 

Claus, RW: Spectral Flame Radiance from a Tubular-Can Combustor, NASA TP-1722, 

1981. 
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Liner Environment Effects Study Objectives 


i 


Establish Broad Data Base on the Effects of Combustion System 
Environment on Combustor Liner Temperatures and Flame Radiation 
Characteristics for a Variety of Liner Concepts and Fuel Properties. 

Data Will Provide Basis for 

• Detailed Combustor Modeling, and 

• Combustor Design 


Liner Cooling Designs 

• Uncooled 

• Film Cooled 

• Impingement & Film Cooled 

• Ceramic Liner 

• Multi-Hole Liner 
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Uncooled & Convectively Cooled Liner 


Cooling Cooling 

Air In Air Out 
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Hot Flow *- 

16 T/C’s 

8 Axial & 2 
Circumferential 
Locations 


Cooling Slot 


Film Cooled Liner 


P 

_1 


16 Metal Skin T/C’s 
8 Axial & 2 
Circumferential 
Locations 


12 Film T/C’s 
6 Axial & 2 
Circumferential 
Locations 
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Impingement/Film Cooled Liner 


Impingement 

Baffle 


1 2 Film T/C’s 1 16 Metal Skin T/C’s 


Ceramic Liner 


Film Cooling Slot 
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3 — « X X X X X 


ID 


Cold Side T/C’s 1 


Embedded 
Hot Side T/C’s 


Ceramic - 
Silicon/Silicon 
Carbide (SILCOMP) 
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Multi-Hole Liner 
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Test Rig Assembly 
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Combustor Section 


»i 




Machined 
Ring Liner 

i 

3" 

r 



Water Outlet 
Ports at 10 
Circumferential 
Locations 


Ceramic Coating 

Base Coat - .002-.005 
NiCrAly 

Top Coat - Magnesia 
Stabilized Zirconia 
Total Coating - .018-.022 


Spectrometer View Port 
2 Locations, 180° Apart 

Wide Band Radiometer 
Port at 1 Location 
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Test Section 


Cooling Air In 



Bleed Air Exit and 
Instrumentation 
Lead Out 


Test Liner 


Exhaust Section 


Port to*- Wide 
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Test Variables 


Parameters 

Liner 

Fuel Hydrogen , Wt. % 

Fuel Nozzle Spray Angle, Deg 

Fuel Nozzle Spray, Sauter Mean 
Diameter, Microns 

Equivalence Ratio 

Cooling Flow Rates (at 2.1 MPa), kg/s 
Cooling Flow Temperature, K 
Internal Reference Velocity, m/s 
Bleed Flow Rates(at 2.1 MPa), kg/s 
Pressure, MPa 


Values 
5 Concepts 
14.0, 12.8, 11.8 
45, 100 

75, 150 

0.3, 0.5, 0.8, 1.2, 1.3 
0.14, 0.23, 0.32 
589, 700, 811, 1000 
9.1, 18.3, 30.0, 41.0 
0.18, 0.32, 0.45 
0.7, 1.4, 2.1 


Test Instrumentation 


Temperatures 

Burner Air Inlet 
Cooling Air Inlet 
Bleed Air 
Air Baffle 

Cooling Air Passage 
Liner Metal ( 1 6) 

Liner Film ( 1 2) 

Radiation 

Total (Wide Band) - 2 

Spectral - Infrared Fourier 
Transform Radiometer 


Flows 

Burner Air 
Cooling Air 
Bleed Air 
Fuel 

Gas Sample 

Exit (7 Element Probe) 
Pressures 

Burner Inlet (Pt and P3) 
Cooling Air Passage 
Inner Impingement Plate 
Exit 


Humidity 

Inlet 
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Droplet Size Measurement 
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Tests will be Conducted with the In-House Droplet 
Measurement Devices at 1 atm 

Air-Assist Fuel Nozzle with Several Combinations 
of Swirlers and Simplex Tips 

Establish Air Flow Requirements to Achieve Desired 
SMD’s - Operational Considerations Likely to Require 
Compromise 


Air-Assist Fuel Nozzle 



Spectral Radiation Measurement Scheme 

P=* Mirror 

Chopper 



Control and Data Aquisition 
Spectral Radiation Measurements 
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SPONTANEOUS IGNITION CHARACTERISTICS OF 

GASEOUS HYDROCARBON-AIR MIXTURES 

G. Freeman and A. H. Lefebvre 
Purdue University 


Experiments are conducted to determine the spontaneous ignition delay times of 
gaseous propane, kerosine vapor, and n-heptane vapor in mixtures with air, and 
oxygen-enriched air, at atmospheric pressure. Over a range of equivalence ratios 
from 0.2 to 0.8 it is found that ignition delay times are sensibly independent of 
fuel concentration. However, the results indicate a strong dependence of delay 
times on oxygen concentration. The experimental data for kerosine and propane 
demonstrate very close agreement with the results obtained previously by Mullins 
and Lezberg respectively. 


INTRODUCTION 

Spontaneous ignition is a process whereby a combustible mixture undergoes 
chemical reaction leading to the rapid evolution of heat, in the absence of any 
concentrated source of ignition such as a flame or spark. In some practical 
combustion devices, such as the lean premix/prevaporize combustor which is now 
being actively developed for advanced aircra c t gas turbine engines, spontaneous 
ignition must be avoided at all costs, since it could damage combustor components 
and produce unacceptably high levels of pollutant emissions (refs. 1 and 2). On 
the other hand spontaneous ignition is sometimes relied upon to provide the main 
source of ignition, as in the compression ignition-engine. 

Spontaneous ignition delay may be defined as the time interval between the 
creation of a combustible mixture, say by injecting fuel into air at high tempera- 
ture, and the onset of flame. In view of their practical irnortance, measurements 
of ignition delay times have been conducted for many fuels, ;i 'er wide ranges of 
ambient conditions and in a variety of test vehicles, including constant volume 
bombs, rapid-compression machines, reciprocating engines, shock tubes, and continu- 
ous flow devices (refs. 3-14). The test methods employed and the results obtained 
are described in reviews by Mullins (ref. 3) and, more recently, by Spadaccini and 
Te Velde (ref. 5). 

Although the subject of spontaneous ignition of liquid fuels has received 
considerable attention in the past, the role of fuel evaporation in the overall 
spontaneous ignition process is still unclear. It was decided, therefore, to 
conduct a series of experiments using only gaseous or vaporized fuels, sometimes 
with the injection of additional oxygen into the flowing air stream, in order to 
obtain more accurate data on chemical delay times under conditions where the effects 
of fuel evaporation are excluded. Having determined the role of chemical kinetics 
in the spontaneous ignition process, the second phase of the research would then 
comprise a number of carefully-designed experiments, using one or more liquid fuels, 
in which fuel evaporation time is varied systematically over a wide range while 
maintaining the chemical reaction time constant. This could be achieved, for 
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example, by conducting a series of measurements of ignition delay time, for several 
different values of mean fuel drop size (between 20 and 200 microns), while main- 
taining constant values of air flow rate, fuel flow rate, air pressure and air 
temperature. The results of these experiments, in conjunction with theoretical 
analyses, should then lead to expressions for the spontaneous ignition times of 
hydrocarbon fuel/air mixtures that take full account of all the physical and 
chemical properties of relevance to the gas turbine. 

The main purpose of this paper is to present the results of some preliminary 
tests that have been carried out at atmospheric pressure using gaseous fuel-air 
mi xtures. 


EXPERIMENTAL 

The apparatus employed in the measurement of spontaneous ignition delay time 
is shown schematically in Fig. 1. Air is supplied from a high pressure blower, the 
air flow rate being controlled by an electrically-activated pintle valve fitted 
into the blower inlet duct. Air flow rates are calculated from the dynamic pressure 
measured in the throat of a bell -mouth orifice through which inlet air enters the 
b 1 owe r . 

Two separate heaters are employed to raise the temperature of the incoming air. 
A preheat section contains two heaters which consist of coiled ni-chrome resistance 
wire would helically around four ceramic support tubes. Each of these preheaters 
operates at 3.5 kW, and can be turned on or off independently of the main heater. 

The main air heater consists of three hairpin-formed, incal loy-sheathed, "cal-rod" 
heating elements. These elements have a maximum sheath temperature of 1150 K and 
a maximum input power of 15 kW. Feedback proportional control is used on the main 
heater. Both heating sections are heavily insulated. 

The fuel injector designed for the experiments on gaseous fuels and fuel vapors 
is illustrated in Fig. 2. For ease of manufacture and assembly, and to ensure 
structural integrity, a square tube matrix was constructed to provide the best 
distribution of holes as dictated by the round cross-section of the pipe. The 
location of each of the twenty-five fuel injection points was determined by dividing 
the cross-sectional area of the pipe into twenty-five sections of equal area and 
locating an injection point as close as possible to the center of each section. 

The objective was to achieve rapid mixing of fuel and air and to attain a uniform 
mixture strength in the fuel/air mixture entering the test section. 

The fuel delivery system for vaporized fuels is illustrated in Fig. 3. Fuel 
is pumped from a reservoir to the fuel vaporizer/heater in the test cell. The 
liquid fuel-flow rate is monitored using a rotameter and is regulated by means of a 
throttle valve. A back pressure valve is installed at the vaporizer exit to 
eliminate the tendency of the vaporized fuel flow to surge due to small pressure 
fluctuations in the vaporizer. Since the vaporizer does not respond well to flow 
transients, the fuel cut-off valve (located just upstream of the injector) is 
equipped with an exhaust line to allow continuous fuel flow, as well as visual 
determination of complete vaporization. 

The test section consists of a single length of drawn stainless steel tubing, 
having an internal diameter of 6.22 cm. Various lengths of tubing are available 
and can be joined together or interchanged to vary the length of the test section. 
Each segment is jacketed by a thick layer of insulating material to minimize the 


axial temperature gradient in the flowing mixture. The test section is instrumented 
with four shielded, stagnation-type thermocouples, located along the axis of the 
tube. The axial locations of the four thermocouples are as follows: 1 - at the 

fuel injection plane, 2 - 18 cm downstream of the fuel injection plane, 3 - halfway 
down the test section, and 4 - at the end of the test section. Thermocouples 3 and 
4 are withdrawn when ignition delay times are being measured. 

Test Procedure 

The experimental procedure developed to determine ignition delay times is as 
follows. Initially flow conditions are set at predetermined values of velocity, 
temperature, and equivalence ratio, and the mixture temperature is then slowly 
increased. The onset of spontaneous ignition is manifested as a clearly visible 
cool flame, accompanied by the odor of aldehydes in the exhaust gases. (Mullins 
reported this same phenomenon (ref. 3)). Spontaneous ignition first appears at a 
distance of thirty to fifty centimeters downstream of the exit of the test section 
as intermittent flashes of flame. As the mixture temperature continues to increase 
the flame becomes stable and its location moves slowly upstream. The flame assumes 
the shape of a cone with its apex pointing upstream, and has an attached flame 
brush extending downstream. The mixture temperature is allowed to increase until 
the tip of the flame is located near to the exit of the test section. At this point 
the inlet air temperature and the mixture temperature at one or more axi 1 locations 
are recorded. Fuel flow is abruptly terminated, and temperature are again recorded 
at various axial locations in the test section. The delay time is then obtained as 
the length of the test section divided by the average velocity in the test section. 

The initial mixture temperature, T cannot be measured directly because mixing 
takes a small, but finite length, so some other method of determining the initial 
temperature must be employed. The method that was chosen is to use a temperature 
measured at a distance downstream of the fuel-injection plane where the results of 
detailed temperature surveys showed that mixing is just complete. This temperature 
is then extrapolated back to obtain the true initial value of T , using as a 
baseline the axial temperature profile as measured at several stations along the 
test section with the fuel turned off. The procedure is illustrated in figure 4, 
which shews that the initial value of T m is determined by lowering the baseline 
temperature profile until it passes through T re f. T re f is then taken as the 
temperature where the shifted profile intersects the temperature axis. If the test 
section were perfectly insulated this procedure would be superfluous because T m 
would be equal to T re f. However, as a finite heat loss occurs, despite the thick 
layer of external insulation, the above procedure is necessary and is considered 
satisfactory. 

One drawback to this method is that the temperature profile of the reacting 
mixture is not exactly the same as that for air only. By plotting the actual 
mixture temperature profile and subtracting from it the air-only profile a graph is 
obtained depicting the temperature rise due to reaction, as shown in figure 5. To 
minimize errors due to this effect it is necessary to keep x * small, although it 
must be large enough to ensure that the temperature measurement carried out at Xf 
is downstream of the mixing region. It would be desirable to measure the axial 
temperature profile for each data point, but this is impractical due to the 
problems that would arise from insertion of the thermocouple into the flowing 
mixture. Insertion of a thermocouple at the upstream end of the working section 
would affect the fuel-air mixing process, while its presence further downstream 
could easily induce premature autoignition of the mixture. 
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In order to determine whether the presence of the flame front can influence 
ignition delay times by supplying heat or active species to the incoming fresh 
mixture, a nitrogen quench system was installed. This system includes a solenoid 
valve which can be opened and closed very rapidly to inject a transient flow of 
cold nitrogen gas into the test section. Injection of nitrogen causes the auto- 
ignition flame to vanish, but when the nitrogen flow is abruptly terminated the 
flame reappears instantly at its original location. This suggests that the 
presence of the flame has no discernible effect on ignition delay times, which is 
consistent with the observations of Spadaccini (ref. 9). 

Although the flame position tends to remain fairly constant it occasionally 
leaps forward a few centimeters and then, almost instantaneously, returns to its 
original position. This phenomenon was also noticed by Mullins (ref. 3) who 
attributed it to inadequate mixing of fuel and air, resulting in pockets of rich 
mixture in the flow approaching the flame. This explanation is considered unsatis- 
factory, at least in its application to the present investigation, partly because 
the results of many detailed temperature surveys have shown that mixing is complete 
well upstream of the autoignition flame, and also because it is generally found that 
ignition delay times are fairly insensitive to variations in fuel/air ratio. No 
attempt has been made to study this phenomenon because it only occurs during the 
phase when the mixture temperature is being slowly raised in order to coax the 
autoignition flame toward the exit of the test section. Close to the test section 
the flame position usually remains quite steady so that measurements of ignition 
delay time can be made with good accuracy. 

RESULTS 


From global reaction rate considerations the chemical ignition delay time can 
be expressed as 

t ct exp (E/RT m ) [Fuel ] m [Oxygen] n (1) 


where 


E 

R 

[Fuel] = 
[Oxygen] = 


ignition delay time, ms 
initial mixture temperature, K 
global activation energy, kcal/kg. mol. 
universal gas constant, kcal/kg. mol -K 
fuel concentration 
oxygen concentration 


The form of equation (1) suggests that a plot of In t versus 1/T should yield 
a straight line with a positive slope, and this is borne out by the results shown 
plotted in figure 6 for propane-air mixtures at a constant equivalence ratio of 0.5. 
The line displayed in this figure is drawn through three sets of data points 
corresponding to three different lengths of test section. The value of activation 
energy, E, calculated from the slope of this line is 38,200 kcal/kg. mol. 


It is of interest to observe that all the data points lie in close proximity 
to the line drawn in figure 6. This is in marked contrast to the results of a previous 
investigation where it was found that ignition delay time was very dependent on the 
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length of the test section employed (ref. 15). This problem was attributed to 
excessive heat loss from the duct walls and was overcome effectively by increasing 
the internal diameter from 2.54 to 6.22 cm, and by applying a thick layer of 
insulating material to the outer pipe wall. 

The results obtained using mixtures of vaporized kerosine with air are shown 
in Fig. 7. The straight line drawn through the data points in this figure corres- 
ponds to a value for E of 40,900 kcal/kg. mol. 

Equivalence ratio was varied over a range from 0.2 to 0.8. For propane the 
results showed that equivalence ratio has an almost negligible effect on x, although 
an increase in equivalence ratio tends to produce a 'stronger' tail flame. 

Equivalence ratios hiqher than 0.8 could not be used due to problems of flashback. 

It was found that as the auto-ignition flame located downstream of the test section 
was brought slowly upstream by gradually increasing the mixture temperature, it would 
suddenly flash upstream and stabilize on the fuel injector, thereby obscuring any 
results and incurring risk of damage to the injector. 

Mixtures of vaporized n-heptane and air also exhibited only a slight dependence 
of ignition delay time on equivalence ratio. Analysis of the results showed that 
x a [fuel ]-0 . 23 _ 

The results obtained for mixtures in air of propane, vaporized n-heptane, 
and vaporized Jet A, are shown in Fig. 8. The equivalence ratio is 0.5 in all 
cases. From inspection of this figure it is apparent that all three lines run 
roughly parallel to each other, thereby indicating that all three fuels have about 
the same value of activation energy. Moreover, it is also clear that ignition delay 
times for propane are appreciably longer than those for Jet A. 

Figures 9 and 10 show for propane and vaporized kerosine respectively the 
effect of replacing some of the nitrogen in the air by oxygen, while maintaining 
the fuel concentration constant. Analysis of the results shows that for propane 
x a [oxygen]"®-^, w h i 1 e for kerosine x a [oxygen]“0-65. 

In figure 11 the experimental data obtained in the present investigation for 
propane-air and Jet A-air mixtures are plotted, for the purpose of comparison, 
alongside the results of previous workers as compiled by Chiappetta and McVey 
(ref. 6). It is of interest to note that the results for Jet A (kerosine) show 
striking agreement with those of Mullins (ref. 3), who also conducted his experi- 
ments at normal atmospheric pressure. Equally noteworthy is the close consistency 
between the propane data obtained by Lezberg (ref. 14) and the results of the 
present experiments. 


CONCLUSIONS 

From the results of experiments conducted at atmospheric pressure on the 
influence of mixture temperature on the spontaneous ignition delay times of mixtures 
with air of propane, n-heptane vapor and kerosine vapor, the following conclusions 
are drawn. 

1. Ignition delay times can be expressed in terms of mixture temperature, fuel 
concentration, and oxygen concentration, by the equation 

x a exp [E/RT m ] [Fuel ] ni [0xygen] n 
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For the fuels employed in this study the value of E is fairly close to 
40,000 kcal/kg mol. For example, for propane E = 38,200 kcal/kg mol., 
while for kerosine E = 40,900 kcal/kg mol. 

The dependence of ignition delay time on fuel concentration is quite small, 
i.e. m = 0. Of the fuels examined the strongest dependence on fuel concentra- 
tion is exhibited by n-heptane, for which m = -0.23. 

Ignition delay times are strongly dependent on oxygen concentration. For 
propane n = -0.59, while for kerosine n = -0.65. 

The experimental data obtained in the present investigation for propane-air 
and kerosine-air mixtures are remarkably consistent with the results of 
previous studies by Lezberg and Mullins respectively. 
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RANDOM VORTEX METHOD FOR COMBUSTING FLOWS 
C. J. Marek 

NASA Lewis Research Center 


The random vortex method RVM of Chorin has been developed by the 
University of California - Berkeley to compute turbulent, reacting* 
recirculating flows, ref. 1. The RVM method models turbulence from first 
principles, tracking the vorticity and obtaining the interaction of vorticity 
with the bull: flow field. A computer program has been produced called MIMOC, 
Modeling the Interface Motion of Combustion, which can be used to calculate the 
reacting flow field behind a rearward facing step, ref. 2 

Several comparisons between experimental data and calculations have been 
niade, refs. 3-4. The RVM method computes qualitatively good results, but the 
quantitative agreement as yet is not completely satisfactory. Much of the 
difficulty may be caused by the treatment of boundary conditions and the 
techniques used for obtaining statistical averages of velocities and turbulence 
quantities. For the rearward facing step the computed reattachment length 
equals the experimental value as shown in figure 1. However the reverse 
velocity in the recirculation zone is over predicted by 300 percent, figure 2. 
In the calculations, a uniform entrance velocity was assumed with no boundary 
layer at the step lip. This high velocity may be overdriving the reverse flow 
region. The profile shown in figure 2 is the worst agreement obtained. Figure 
3 shows that the time steps used to obtain the statistical average also is 
important. As the calculation continues, stationary values should be obtained 
for the mean values, but the averages must bo taken over the time of several 
large scale vortex sheddings. 

In house calculations are being conducted to look at the comparison 
between the turbulence quantities and the experimental values. 

Recent calculations by Hsiao, ref. 4, show much better agreement between 
computations of the mean velocities and experimental data, figure 4. Hsiao has 
included the inlet passage within the calculation domain. The statistical 
quantities are not in good agreement, figure 5, but the averages were computed 
over only 20 time steps. 

The results from the RVM method have been very encouraging and much can be 
learned through continued study of the calculation method. 
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AERODYNAMIC FEATURES OF FLAMES IN PREMIXED GASES 
A.K. Oppenheim 

University of California, Berkeley 


A variety of experimentally established flame phenomena in premixed 
gases are interpreted by relating them to basic aerodynamic properties 
of the flow field. On this basis the essential mechanism of some well 
known characteristic features of flames stabilized in the wake of a 
bluff-body or propagating in ducts are revealed. Elementary components 
of the flame propagation process are shown to be: (1) rotary motion, (2) 
self-advancement, and (3) expansion. 

Their consequences are analyzed under a most strict set of ideali- 
zations that permit the flow field to be treated as potential in charac- 
ter while the flame is modelled as a Stefan-like interface capable of 
exerting a feed-back effect upon the flow field. The results provide an 
insight into the fundamental fluid-mechanical reasons for the experimen- 
tally observed distortions of the flame front, rationalizing in particu- 
lar its ability to sustain relatively high flow velocities at amazingly 
low normal burning speeds. 
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TRANSIENT FLOW COMBUSTION 

Robert R. Tacir.a 
NASA Lewis Research Center 


t.ju steady combustion problems can result from engine sources such as 
accelerations, decelerations, nozzle adjustments, augmentor ignition, and air 
perturbations into and out of the compressor. Also non-steady combustion can be 
generated internally from combustion instability or self-induced oscillations. 
A premixed-prevaporized combustor would be particularly sensitive to flow 
transients because of its susceptability to flashback-autoignition and blowout. 
An experimental program, the "Transient Flow Combustion Study” is in progress to 
study the effects of air and fuel flow transients on a premixed-prevaporized 
combustor. 


This experiment is to be performed at LeRC CE5b, test stand 3. The 
transient capability is provided by high response hydraulic valves (full 
traverse in 100 ms) that are installed upstream and downstream of the test 
section to control airflow and pressure and in the fuel line to control fuel 
flow. Figure 1 shows the salient features of the test rig and the hardware. 
Nonvitiated preheated air flows through the upstream control valve to a 
diffuser section with straightening vanes to provide a uniform velocity profile 
to the test section. The inlet velocity profile is measured with a traversing 
total pressure probe. Test hardware for lean premixed-prevaporized combustor 
tests is shown. The airflow splits at the test section inlet. Half the air goes 
through the center zone for premixing- prevaporizing with the fuel and 
downstream burning. The other half enters an outer annulus and is used for film 
cooling in the downstream burning section. In the center zone, fuel is 
distributed uniformly by a concentric multi-point fuel injector. The fuel 
injector can be located either 15 cm or 38 cm upstream of the flameholder. The 
fuel-air distribution is measured in a plane 4 cm upstream of the flameholder 
with a traversing probe. High response pressure transducers, thermocouples and a 
photodetector are installed in the premixing-prevaporizing section to monitor 
transient response, in particular to detect if autoignition or flashback should 
occur. Doimstream of the premixing-prevaporizing section there is a perforated 
plate flameholder and a combustion section. The walls of the flame tube are film 
cooled. The exit conditions will be measured with a traversing total pressure, 
temperature and emissions probe. Emission measurements will only be taken at 
steady state conditions. The bad: pressure valve is a specially designed water 
cooled valve. This eliminates the need to cool the combustor exhaust products 
with water spray upstream of the back pressure valve. Hater spray upstream of 
the valve can introduce flow transients if there is two-phase flow. 

Preliminary tests have been performed at an inlet air temperature of 600 K, 
a reference velocity of 30 m/s, and a pressure of 700 l:Pa. Ramping the airflow 
down (at constant fuel-air ratio) indicated that this combustor is less 
sensitive to flashback than anticipated. The airflow was reduced to 1/3 of its 
orginal value in a 40 ms ramp before flashback occurred. Ramping the airflow up 
has shown that blowout is more sensitive than flashback to flow transients. 
Blowout occurred with a 25 percent increase in airflow (at a constant fuel-air 
ratio) in a 20 ms ramp. Combustion resonance has been found at some conditions 
and may be important in determining the effects of flow transients. 
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COMBUSTOR FLAME FLASHBACK 
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Flashback, a problem that occurs in premixcd-prevaporized combustors, is the 
upstream propagation of the flame from the combustor into the premixing tube. Not 
only dees flashback change the combustion process from premixed burning to diffusion 
burning, thus creating more pollutants, but it also inflicts considerable damage to 
the fuel injector, premixing tube and other equipment upstream. This study attempts 
to define the conditions at which flashback occurs in steady burning and the mechanism 
that causes flashback in both steady and transient flow. 

Part of the experimental setup is shown in the attached figure. The stainless 
steel test section is a two-dimensional (rectangular cross-section) center dump 
combustor. The 4" x 1" premixing tube, 6"-long, dumps into the 4"x4" combustor, 
10"-long. The combustor empties into the exhaust plenum through a 2"-diameter hole 
that may be partiallv blocked by a conical plunger. Movement of the plunger, powered 
by two solenoids, creates a pressure pulse for transient data. The fuel injector 
stationed IV* upstream of the premixer in the transition section supplies gaseous 
propane to the test section. One or both sides of the test section may be fitted 
with V'-thick quartz windows to view the entire combustion process. 

The equivalence ratio at which flashback occurs is being measured for inlet 
temperatures of 600-910 K, premixer wall temperatures of 450-1050 K and premixer 
velocities of 40-80 ft's. These data will be presented. It is hoped that by the 
time of the conference high speed film will be available to show the mechanism of 
flashback. Future work involves the buildup of a slightly modified Hastelloy test 
section that will withstand inlet temperatures of 1100 K. Transient data will also 
be collected. 
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